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We study periodic structures, such as vortex lattices, moving in a random pinning potential 
under the action of an external driving force. As predicted in [T. Giamarchi, P. Le Doussal Phys. 
Rev. Lett. 76 3408 (1996)] the periodicity in the direction transverse to motion leads to to a 
new class of driven systems: the Moving Glasses. We analyse using several renormalization group 
techniques the physical properties of such systems both at zero and non zero temperature. The 
Moving glass has the following generic properties (in d < 3 for uncorrelated disorder) (i) decay of 
translational long range order (ii) particles flow along static channels (iii) the channel pattern is 
highly correlated along the direction transverse to motion through elastic compression modes (iv) 
there are barriers to transverse motion. We demonstrate the existence of the "transverse critical 
force" at r = and study the transverse depinning. A "static random force" both in longitudinal 
and transverse directions is shown to be generated by motion. Displacements are found to grow 
logarithmically at large scale in d = 3 and as a power law in d = 2. The persistence of quasi long 
range translational order in d = 3 at weak disorder, or large velocity leads to the prediction of the 
topologically ordered "Moving Bragg Glass" . This dynamical phase which is a continuation of the 
static Bragg glass studied previously, is shown to be stable to a non zero temperature. At finite but 
low temperature, the channels broaden but survive and strong non linear effects still exist in the 
transverse response, though the asymptotic behavior is found to be linear. In d = 2, or in d = 3 at 
intermediate disorder, another moving glass state exist, which retains smectic order in the transverse 
direction: the Moving Transverse Glass. It is described by the Moving glass equation introduced 
in our previous work. The existence of channels allows to naturally describe the transition towards 
plastic flow. We propose a phase diagram in temperature, force and disorder for the static and 
moving structures. For correlated disorder we predict a "moving Bose glass" state with anisotropic 
transverse Meissner effect, localization and transverse pinning. We discuss the effect of additional 
linear and non linear terms generated at large scale in the equation of motion. Generalizations 
of the Moving glass equation to a larger class of non potential glassy systems described by zero 
temperature and non zero temperature disordered flxed points (dissipative glasses) are proposed. 
We discuss experimental consequences for several systems, such as anomalous Hall effect in the 
Wigner crystal, transverse critical current in the vortex lattice, and solid friction. 



I. INTRODUCTION 



Interacting systems which tend to form spontaneously periodic structures can exhibit a remarkable variety of 
complex phenomena when they are driven by an external force over a disordered substrate. Many of these phenomena, 
which arise from the interplay between elasticity, periodicity, quenched disorder, non linearities and driving, are still 
poorly understood or even unexplored. For numerous such experimental systems, transport experiments are usually 
a convenient way to probe the physics (and sometimes the only way when more direct methods - e.g. imaging are 
not available) . It is thus an important and challenging problem to obtain a quantitative descriptiipn of their driven 
dynamics. Vortex lattices in type II superconductors are a prominent example of such systemsEI. The motion of 
the lattice under the action of the Lorentz force (associated to a transport supercurrent) in the presence of pinning 
impurities has been studied in many recent experimentgjil3 There are other examples of well studied driven systems 
where quenched disorder is kapjAai to be important, such as the two dimensional electron gas i n M a . magnetic field 
which forms a Wigner crystalljilj moving under an applied voltage, lattices oL magnetic bubblesr tM moving under 
an applied magnetic field gradient. Charge Density Waves (CDW)E3 or colloidai3 submitted to an electric field, driven 
josephson junction arrayst3't2l etc.. This problem may also be important in imderstanding tribology and solid friction 
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pheiiQiaenaEJ, surface growth of crystals with quenched bulk or substrate disorderEil, domain walls in incommensurate 
solidsEj. One striking property exhibited by all these systems is pinning, i.e the fact that at low temperature there 
is no macroscopic motion unless the applied force / is larger than a threshold critical, ^ce fc- Dynamic properties 
have thus-|heen studied for some time, quite extensivelii-|Uear the depinning thresholdL3~Lj, but mostly in the context 
of CDWcH"^3 or for models based on jirijien manifoldscJ^Ey and their relation to growth processed^ described by the 
Kardar Parisi Zhang (KPZ) equationE3'E£l They are however, far from being fully understood. In addition, the full 
problem of a periodic lattice (with additional periodicity transverse to the direction of motion) was not scrutinized 
until very recentlyEj. 

A crucial question in both the dynamics and the statics is whether topological defects in the periodic structure 
are generated by disorder, temperature and the driving force or their combined effect. Another important issue is to 
characterize the degree of order (e.g translational order, or temporal order) in the structure in presence of quenched 
disorder. In the absence of topological defects it is sufficient in the statics to consider only elastic deformations. In 
the dynamics this leads to elastic flow. On the other hand, if these defects exist (e.g unbound dislocation loops) the 
internal periodicity of the structure is lost and one must consider also plastic deformations. In the dynamics the flow 
will then not be elastic but turn into plastic flow with a radically different behaviour. 

The statics of lattices with impurity disorder has been much investigated recently, especially in. the context of 
type II superconductors. It was generally agreed that disorder leads to-a^glass phase (often calledo a vortex glass) 
with many metastable staij^—diverging barriers between these statea^J'tl, pinning and loss of translational order. 
Indeed, general argument£j'Lj, unchallenged until recently, tended to show that disorder would always favot-tlie 
presence of dislocations destroying the Abrikosov lattice beyond some length scale. In a series of recent workaSaEll, 
we have obtained a different picture of the statics of disordered lattices (including vortex lattices) and predicted 
the existence of a new thermodynamic phase, the Bragg glass. The Bragg glass has the following properties: (i) it 
is topologically ordered (ii) relative displacements grow only logarithmically at large scale (iii) translational order 
decays at most algebraically and there are divergent Bragg peaks in the structure function in d = 3 (i.e quasi long 
range order, survives), (iv) it is never|thelfiss a true static glass phase with diverging barriers. There has been several 
analyticanj^LJ and numerical studie^O confirming this theory. Tiie predicted consequences for the phase diagram 
of superconductors compare well with the most recent experimentscil 

While some progress towards a consistent theoretical treatment has been made in the statics, it is still further 
removed in the dynamics. Determining the various phases of driven system is still a widely open question. Evidence 
based mostly on experiments, numerical simulations and qualitative arguments indicates that quite generally one 
expects plastic motion for either strong disorder situations, high temperature, or near the depinning threshold in-ilpsf 
dimensions (for CDW see e.g.c3). Indeed there has been a large number of studies on plastic (defective) flowc3~Ej. 
In the conlext of superconductors a H-T phase diagram with regions of elastic flow and regions otiJaatkL flow was 
observedEllEl. Several experimental effects have been attributed tp-plastic flow, such as the peak effectt3Elc3o, unusual 
broadband noiseE3 and fingerprint phenomena in the I-V curveEaElliiO. Steps in the I — V curve were also observed 
in YBCO near melting irJ3. Close to the threshold aacLin strong disorder situations the depinning is observed to 
proceed through what can be called "plastic channels"E3E3 between pinned regions. This type of filamentary flow has 
been found iiM in simulations of 2D (strong disorder) thin film geometry (with cn S> cge)- Depinning then proceeds 
via filamentary channels which become increasingly denser. Filamentary flow was proposed as an explanation for 
the observed sharp dynamical transition observed in MoGe filmsrIrH characterized by abrupt steps in the-differential 
resistance. Also, interesting effects of synchronisation of the flow in different channels were observed inEj. Despite 
the number of experimental and numerical datacj'E3 a detailed theoretical understanding of plastic motion remains 
quite a challengeEj. 

As in the statics, one is in a better position to describe the elastic flow regime, which is still a difficult problem. 
This is the situation on which we will focus in this paper. Though elastic flow in some cases extends to all velocities, a 
natural idea was to start from the large velocity region and carry perturbation theory in 1/v. At large velocity one may 
think at first that since the sliding system averages enough over disorder one recovers a simple behavior, in fact mucli 
simpler than in the statics. ludeed it was observed experimentally, some time ago in neutron diffraction experimentsEj, 
and in more details recentlyE3, that at large velocity the vortex lattice is more translatiatLaIlii_ordered than at low 
velocity. This tendency to dynamical reordering has also been seen in numy pcal simulationsHc3o. The 1/v expansion 
has been fruitful to compute the corrections to the velocity itself inE3C3'EZl. Recently it was extended by Koshelev 
and Vinokur inLj to compute the vortex displacements u induced by disorder and leads to a description in term of an 
additional effective shaking temperature induced by motion. This description suggests bounded displacements in the 
solid and thus a perfect moving ctystal at large velocity. 

Recently we have investigatedE3 the effects of the periodicity of the moving lattice in the direction transverse to 
motion, in the same spirit which led to the prediction of the Bragg glass in the statics. It was still an open problem 
how much of the glassy properties remain once the lattice is set in motion. We found that, contrarily to the naive 
expectation, some modes of the disorder are not affected by the motion even at large velocity. Thus, the large v 
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expansion oiS breaks down and the effects of non linear static disorder persists at all velocities, leading to new 
physics. As a result the moving lattice is not a perfect crystal but a moving glass. r—. 

The aim of this paper is to provide a detailed description of the moving glass state predicted iii^ and to present 
our approach to the general problem of moving lattices. A brief account of some of the new re s ul ts contained here (e.g 
the r = rcnormalization group equations RG and fixed points) has already appeared inHllJ. We will use several 
RG approach at zero and at non zero temperature. Since several Sections of this paper are rather technical we have 
chosen to expose all the results about the physics of the moving glass in Sections || and III in a s elf con tained manner, 



avoiding all technicalities. The reader can find there t he re sults for the existence of channels ( III A ) the trans verse 



I- V curves at T = and the dynamical Larkin length ( III B ) , the ran dom force and the correlation functions ( III C 
) the various crossover lengths and the the finite temperature results ( III E|) . Decoupling scenarios which distinguish 



between the Moving Bragg glass and th e Moving transverse glass ( p^II D| ) as well as predictions for the dynamical 
phase diagrams are also given in ( [IIIFD Finally we discuss how the moving glass theory stands presently compared 



to numerical simulations ( IIIG) and experiments (IIIH) and present some suggestions of further observables which 
would be interesting to measure. 

The following Sections are devoted to making progress in an an analytic description of the moving state of interacting 
particles in a random potential. Since this is a vastly difficult problem, it is potentially dangerous (and unfruitful) to 
try to attack this problem by treating all the effects at the same time (dislocations, non linearities, thermal effects 
etc.). Already within the simplifying assumption of an elastic flow two main types of phenomena are missed in a 
naive large v approach. The first one is a direct consequence of previous works on driven dynamics of CDW and 



elastic manifoldsEaO. It is expected on symmetry groundsllithat new non linear KPZ terms (Vu)^ will be generated 
by motion, an effect which was studied in the driven liquidO. Another important effect, stud-ied so far only within 
the physics of CDW, is the generation of a static random force convincingly argued by KrugIS and explored iiJU. If 
both effects are assumed to occur simultaneously, they may lead to interesting interplays which have been explored 
only recently and only in simple CDW modeldlj. However there still no explicit RG derivation of those terms even 
in CDW models. In the context of driven lattices, there are not even discussed yet. Our aim in this paper is to 
remedy this situation. We derive these terms explicitly and show that other linear terms, a priori even more relevant 
are generated. Though these additional linear, non linear and random force terms certainly complicate seriously the 
problem, focusing exclusively on these terms only obscurs the physics of the present problem. Indeed the second and 
as we show here more important effect in the moving structure is the crucial role of transverse periodicity to describe 
the dynamics. 

A rigorous study of the problem of moving interacting particules would be to first study the fully elastic flow 
of a lattice. Once the main elastic physics is understood a second step is then to allow for topological excitations 
(vacancies, interstitials, dislocations). In principle the results obtained within the elastic only approach can, as in the 
statics, be used to check self consistently the stability of the elastic flow itself. It is hard to see how one can do that 
in a controlled way without some detailed understanding the elastic flow first ! Here we carry most of the first step 
and propose an effective description of the second. 

Even the purely elastic model turns out to be difficult to treat when all sources of anisotropics, non linear elasticity, 
cutoff effects are included. There are no analogous terms in the statics and thus in that sense the dynamics is more 
difficult. Our strategy has thus been to simplify the problem in several stages and resort to simplified models. The 
simplified models of Moving glasses that we have obtained turn out to exhibit some new physics and become interesting 
in their own. They call for interesting generalizations to other models exhibiting dissipative glassy behaviour, as we 
will propose. 

We will call Model I the full model of an elastic flow of a lattice containing all the above mentionned relevant linear 
and non linear terms. Such models can also be written for oth er elast ic structures with related kind of order (such 



as liquid crystal order) . This model will be discussed in Section VIII B . However its complete study goes beyond the 
present paper. 

Fortunately, a useful and further simplified model can be constructed (Model II) . It corresponds to considering the 
above full elastic model in the continuum limit. It will certainly give a very good approximation of the full model at 
least up to some very large scale. This model was discussed inE3 and is studied in details here. It has both longitudinal 
degrees of freedom (along the direction of motion) and transverse ones. Though it is quite difficult, it can be handled 
by perturbative rcnormalization group studies, as we will show here. It has a new and non trivial fixed point which 
gives a detailed description of the moving Bragg glass phase. 

It turns out that most of the physics of Moving glass is contained in a further simplification of Model II which retains 
only the transverse degrees of freedom (displacements). This model, which here we call Model III, was introduced in 
Rcf. ^ and is described by the equation of motion: 

r]dtu + rjvd^u = cV^u + F"*^* (r, u(r, t)) + C(r, t) (1) 



3 



which we call the moving glass equation. F^*^* is a non linear static pinning force and we have denoted x the direction 
of motion, y the transverse direction(s) and r = {x,y). The model retains only the transverse displaceemnt u = Uy. 

Equation (|l|) was obtained simply by considering the density modes of the moving structure which are uniform in 
the direction of motion. Indeed, the key point of Ref. ^ is that the transverse physics is to a large extent independent 
of the details of the behaviour of the structure along the direction of motion. This is because the transverse density 
modes, which can be termed smectic modes, see an almost static disorder and thus are the most important one to 
describe the physics of moving structures with a periodicity in the direction transverse to motion! Let us emphasize 
that this is explicitly confirmed hece by the detailed RG analysis of the properties of Model II. Note that to obtain 
Model III one sets formally = OIj. 

The hierarchy of models introduced here is represented in Fig. |^. 



Moving Bragg glass 

Elastic flow of lattice 
no topological defects 

Model I full elastic 
V model 



add topological defetcs 



plastic flow 
filamentary flow 



\ add KPZ 
\ linear terms 

Model II elastic limit 
delta correlated disorder 

add 



add some dislocations 
channel decoupling 



Moving Transverse glass 

Smectic transverse order 
of channels 



effective description 
of transverse physics 



loss of order 
of channels 



11,66,44 \ Model III 

moving glass equation u y 
T=0 fixed point 
T>0 fixed point 

FIG. 1. Various models studied here to describe with various levels of approximation the (i) fully elastic flow of a lattice (ii) 
intermediate phase with ordering transverse to motion (iii) plastic flow 



The outline of the paper is as follows. After the Sections || and III where we give a non technical discussion of the 
physical results, we start in Section^ by deriving an equation of motion and, carefully examining its symmetries, we 
introduce the Models I, II, III and explain the approximations leading to them. In Section ^ we perform perturbation 
theory on the full dynamical problem, focusing on Model II. Iir Section VI we use the functional RG to s tudy Model 
III and thus the transverse physics in d = 3 and d = 3 — e. We study T — and T > 0. In Section VII we study a 
two dimensional version of the moving glass equation Model III. This allows to obtain results in d = 2 at T > and 
in d = 2 + e. H aving obtained a good understanding of the trairsverse physics in the previous sections, we attack in 



Section VIII A the RG of Model II. Finally in Section VIII B we examine the full Model I, show that linear terms and 



KPZ terms are generated at large scales and discuss some consequences. 



II. MOVING STRUCTURES AND MOVING GLASSES 



A. Moving structures: general considerations 



All the structures we consider share the same basic features. The static system in the absence of quenched substrate 
disorder consists of a network of interacting objects at equilibrium positions i?^, forming either a perfect lattice 
(periodic case) or elastic manifolds (non periodic case). Depending on the system the objects can be either pointlike 
(e.g. electrons in a Wigner crystal) or lines (vortex lines in superconductors ). Deformations away from equilibrium 
positions are described by displacements Ui or in a coarse grained description u(r, t) where r is the internal coordinate. 
A complete characterization of the structure in motion uses three parameters (i) the internal dimension D (ii) the 
number of components n of the displacement field Ua and (iii) the embedding space dimension d. Two examples 
are shown in Figure ^ and more details are given in the Appendix. Since we are mostly interested here in periodic 
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structures (though not exclusively) we can set D = d. We will consider motion along one direction called x, and we 
parametrize throughout all this paper the space variable r as r = {x, y, z) where x is one dimension, y has a priori 
n — 1 dimensions and z has dz = d — n dimensions, and the displacements along motion as and transverse to 
motion as Uy. Three dimensional triangular flux lines lattices driven along a lattice direction thus have d — S^n — 2, 
r = (x,y, z), u — {ux, Uy) where z denotes the direction of the magnetic fields. Two dimensional triangular lattices of 
point vortices have (i = 2,7i = 2, r — {x,y). 




At finite temperatures or in the presence of quenched substrate disorder the structure is deformed. An important 
issue is then to characterize the degree of order. This can be expressed in terms of displacements correlation functions. 
The simplest one measures the relative displacements of two points (e.g two vortices) separated by a distance r. 

B{r) = l{[u{r)~u{m (2) 
n 

where () denotes an average over thermal fluctuations and is an average over disorder. The growth of B{r) with 

distance is a measure of how fast the lattice is distorted. For thermal fluctuations alone in d > 2, B{r) saturates at 
finite values, indicating that the lattice is preserved. Intuitively it is obvious that in the presence of disorder B{r), 
will grow faster and can become unbounded. B{r) can directly be extracted from direct imaging of the lattice, such 
as performed in decoration experiments of flux lattices. 

Related to B{r) is the structure factor of the lattice, obtained by computing the Fourier transform of the density 
of objects p{r) — J^i^'^i^ ~ ~ ^i)- The square of the modulus \pk\'^ of the Fourier transform of the density is 
measured directly in diffraction (Neutrons, X-rays) experiments. For a perfect lattice the diffraction pattern consists 
of ^-function Bragg peaks at the reciprocal vectors. The shape and width of any single peak around K can be Fourier 
transformed to obtain the translational order correlation function given by 

CK{r) = {e^K.uir)e-^K.u{0)) (3) 

Ck{x) is therefore a direct measure of the degree of translational order that remains in the system. Three possible 
cases are shown in figure ^. 
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r r r 

(a) (b) (c) 

FIG. 3. Various possible decays of Ciiir). (a) For thermal fluctuations alone CK(r) Cste, one keeps perfect 5 function 
Bragg peaks, albeit with a reduced weight (the Debye- Waller factor), (b) CK(r) decays exponentially fast. The structure 
factor has no divergent peak any more, so translational order is destroyed beyond length Ra, although some degree of order 
persists at short distance, (c) Cx (r) decays as a power law. The structure factor still has divergent peaks but not sharp 5 
function ones. One retains quasi-long range translational order. This is for example the case in d = 2 at small temperature 
(Kosterlitz-Thouless) or in the Bragg glass. 

For simple Gaussian fluctuations (and isotropic displacements) Cii-(r) = e"'^^'-'"-' but such a relation holds only 
qualitatively in general (as a lower bound). Depending on how much cristalline order remains in the system the 
structure factor will have extremely different behaviors as depicted in figure H. 



S(q) 



FIG. 4. Depending on the translational order remaining in the lattice the structure factor has different shapes. The thick line 
is the 5 function Bragg peak of a perfect lattice (including thermal fluctuations). The dashed line is the power-law Bragg peak 
of the Bragg glass (which retains quasi long range order and has no topological defects) , the dotted line is the lorentzian-like 
shape of a system loosing its translational order exponentially fast. 



Quite surprisingly, if one takes into account correctly the p eriod icity of the lattice, a thermodynamic phase without 
dislocations was predicted to exist in d = 3 at weak disorder .E3E3. This phase, named the Bragg Glass, posses quasi 
long range order with Bragg peaks diverging as least as (■witiL_43 « 1), similar to dashed line in Fig. ^ At 

the same time displacements B{r) grow logarithmically at large scaleEaO. Similar predictions hold for other elastic 
models such as random field XX.|S¥stems, and a priori also for tauid crystals. The Bragg glass theory has by now 
received considerable numerica]c3cj and analytical confirmationscSci . 

If disorder is increased above a threshold it is predicted that there is a transition at which topological defects 
proliferate. They destroy the translational long range order exponentially fast beyond a length Ro leading to finite 
height diffractions peaks. The height of the peak will be inversely proportional to the scale at which translational 
order is destroyed. This transition is thus characterized by the loss of the divergence iuptiie Bragg peaks. In type II 
superconductors it implies that there is a transition, upon increasing the magnetic fieldEj, from the Bragg glass (at 
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low fields) to another phase. The high field phase is either the putative vortex glasflS or is simply continuously 
related to the high temperature phase. These predictions for the phase diagrams of superconductors has received 
experimental support (see Ref. ^ for a review). 

What happens when an external force is applied to such a structure ? One obviously important quantity to 
determine is the curve of velocity v versus the applied force /. Through this v — f characterics, three main regimes 
can be distinguished and are shown on Figure ^. 




FIG. 5. A typical v — f characteristics at T = (full line) and finite temperatures (dashed line). Three main regimes can be 
distinguished: the creep regime for forces well below threshold, the critical regime around the threshold and the large velocity 
regime well above threshold. 



Far below the depinning threshold fc the system moves through thermal activation. This is the so called creep 
regimej—Since the motion is extremely slow in this regime, it has been analyzed based on the properties of the static 
systemOS. The resulting v — f curve crucially depends on whether the static system is in a glass state (such as the 
Bragg Glass) where the barriers U(f) become very large when / — > 0, or a liquid where barriers remain finite at small 
/, resulting in a linear resistivity. The general form expected in the creep regime is: 



Po/e 



-U{f)/T 



(4) 



Let us emphasize that this "longitudinal" v ~ f characteristics has mainly be used to determine whether the the 
static system (i.e the limit / = f = 0) is or not in a glass state. It may not be enough though, if one wants to probe 
glassiness of the moving system itself. 

The second regime, nearpthe-depinning transition / ~ /c, has been intensely investigated in similarity with usual 
critical phenomena (see e.gE3~E3) where the velocity plays the role of an order parameter. A particularly important 
question is that regime is to determine whether plastic rather than elastic motion occurs. Close to the threshold, in. 
low dimensions and in strong disorder situations the depinning is observpJ to proceed through "plastic channels"E3'E2l 
between pinned regions. This type of filamentary flow has been found inEj in simulations of 2D (strong disorder) thin 
film geometry (with cn ^ cee) where depinning proceeds via filamentary channels which become increasingly denser. 
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FIG. 6. Plastic flow of a network of objects submitted to an external force F, shown for simplicity in two dimensions. The 
motion occurs through plastic channels around pinned regions. Plastic flow might occur close to the depining threshold whereas 
at large velocities one expects to recover elastic flow where the whole lattice moves coherently. 



The third regime is above the depinning threshold f > fc This is situation on which we will focus in this paper 
(though some of our considerations will have consequences in the other regimes as well) . An important phenomenon 
in this regime is that pf dynamical reordering. Indejed, it was observed experimentally, some time ago in neutron 
diffraction experimental, and in more details recently^, that at large velocity the vortex lattice is more translationally 
ordered than at low velocity. Intuitively the idea is that at large velocity v, the pinning force on each vortejt, varies 
rapidly and disorder should produce little effect. This phenomenatLpSiias also known in the context of CDWE3. The 
tendency to reorder has also been seen in numerical simulationafS'OLj. 

Since the effect of disorder were expected taipiaxush at high velocity perturbation theory in 1/v were developped 
mainly to compute the v — f characteristicg£§E30. Recently it was extended by Koshelev and Vinokur inEll to 
compute the vortex displacements u induced by disorder in the moving lattice and in the moving liquid. The effect 
of disorder on the moving liquid was found to be equivalent to heating to an effective temperature T' = T -\- Tgh 
with Tsh ^ l/v. Thus the moving liquid was argued to survive at temperatures lower than the melting temperature 
T < Tm, and 'A-dynamical melting transition to occur below T,„ from a moving liquid to a moving solid upon increase 
of the velocityEll, when T' = Tm- These arguments were then extended to describe the moving solid itself, and it was 
argued that there the effect of pinning could also be describedEZi by some effective shaking temperature Tgh ^ 1 /v^ 
defined by the relation < |m((7)P >= Tshlcm'f' . This suggests bounded displacements in the solid and that at low T 
and above a certain velocity the moving lattice is a perfect crystal. As will be discussed the remainder of this paper 
the picture of the moving lattices emerging from the above bold qualitative argumentstZI goes wrong in several ways. 

There are several other important questions to be answered in addition to the v — f characteristics. The first one 
is the question of the effect of the motion on the spatial correlations and in particular whether translational order 
exist in a moving system. This is related to the question of plastic versus elastic flow. If plastic flow occurs, the 
structure factor should signal some destruction of lattice. However because a moving system is inherently anisotropic 
new effects may appear and the decay of the structure factor will not be as isotropic as in the static system (the 
Lorentzian in Figure-Q). This question thus remains to be investigated. A possibility, suggested by the idea of a 
shaking temperaturel°J, would be that at large velocity one should observe (S-function Bragg peaks characteristic of 
a crystal a finite temperature. Such questions will be discussed in details in section III. Finally determining how 
motion affects the phase diagram of the statics has to be investigated and depends of course on the above issues. In 
particular what remains of the glassy properties of the systems when in motion (slow relaxation, history dependence, 
non-linear behaviors) needs to be addressed. 

For moving periodic systems, an equivalent question can be asked also about "temporal order" and its associated 
effects such as noise spectrum. In particular if one looks at a signal at a fixed position in space but as a function 
of time, one expects a periodic signal with a periodicity of a/v, having S peaks in frequency at the multiples of the 
washboard frequency ujq = 2'Kv/a. If the lattice becomes imperfect one could naively expect the Fourier peaks in 
frequency to broaden in a way that reflects the loss of translationaL order. Quite surprisingly this is not so. Indeed 
it can be shown for a single component displacements field (CDW)Lj that the perfect periodicity in time remains (in 
the absence of topological defects). However this result is not readily applicable to a moving lattice, and it is thus 
crucial to determine whether this remarkable property holds in that case. 



B. The Moving Glass 

To tackle the physics of a structure with a displacement field with more than one component (n > 1), such as 
a trignOTilgr lattice (by contrast with a single Q CDW), two routes seem to be possible—The commonly followed 
oncOQ'CEl is to simply borrow from, or extend, tte physics of single component CDWocj, or of elastic manifolds 
driven perpendicularly to their internal directioncJ. In this case emphasis is put on the displacements along the 
direction of motion and on the proper way to model its dynamics. Such a problem has turned out to be already 
quite complicated in particular due to the generation of KPZ type non linearities in the equation of motion. Evenjf 
degrees of freedom transverse to motion Uy exist as in the cases depicted in Figure ^ they constitute an extensionE2l 
of this "longitudinal" physics. Thus in this "CDW paradigm" it would seem necessary to understand first completely 
the physics of longitudinal modes and then incorporate Uy as an extra-complication. Indeed there were a few 
attempts to describes the physics of driven vortex lattices along those line£3'Efl. 
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FIG. 7. Three types of dynamical systems, (a) a manifold driven perpendicular to itself, (b) a single Q CDW system. Only 
displacements in the directions of motion exist, but periodicity along the direction of motion can play a role, (c) A periodic 
system with transverse degrees of freedom driven along one of its symmetry directions. The correct description of this last class 
of systems is the moving glass fixed point where the transverse degrees of freedom are the important one as represented here. 

The second approach is based on thc_realization that the physics of periodic structures driven along one of their 
internal direction is radically difFerentO from the above descriptions. This stems from the fact that due to the 
periodicity in the transverse direction Uy a static non linear pinning force F^*'^'^ persists even in a fast moving system. 
We want to stress that this is a very general property of a any moving structure which contains uaiform density modes 
Kx = in the direction of motion (as can be seen on the Fourier decomposition of the densityl£3). As illustrated in 
Fig. 1^ the substructure formed by these modes can deform elastically in the Uy direction and sees essentially a static 
disorder. As is obvious from the picture |^ (c), this substructure has generically a liquid crystal type of (topological) 
order and can be termed a "smectic" (though when dz — 0, e.g for d = 3 and n = 3, it is rather a "line crystal" - see 
below) . In all cases the basic starting point thus involves the transverse degrees of freedom as shown on figure 0, 
and is quite different from the "CDW description" . The equation which capture the main ingredients of such moving 
systems was derived in Ref. It leads to a new and interesting model for transverse components u = Uy, which has 
the general form in the laboratory frame: 

j]dtu + rivd^u = cV^M + i^"*^*(r, u(r, t)) + C(r, t) (5) 

Since this equation captures glassy features of moving systems we call it the moving glass equation. Although it 
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looks like a standard pinning equation the convection term ijvdxU dissipates even in the static limit (a reminder that 
we are looking at a moving system) and does not derive from a potential. Thus we will consider this problem and 
its generalizations as a prototype for a new class of physical phenomena which are glassy and do not derive from a 
potential (or from a hamiltonian) . The first example is to choose F^^'^^{r,u) periodic in the u direction: 

F^'-\r,Uy) = V{r)po ^ KySmKy{uy - y) (6) 

Ky^Q 

and corresponds to lattices (or to liquid crystals) driven in a random potential with a short range correlator 
V{r)V{r') = g{r — r') of range r/. The study of this case in Ref. ^ gave the first hint that non potential dy- 
namics can indeed exhibit glassy properties and lead to dissipative glasses. This is a rather delicate notion because 
the the constant dissipation rate in the system would naturally tend to generate or increase the effective temperature 
and kill the glassy properties. However this type of competition between glassy behaviour and dissipation arises in 
other systems which are a generalization of the above equation. Let us briefly indicate some of the generalizations 
that we are proposing which are being studied here or in related works. l— , 

An interesting generalization is the case of a periodic manifold with correlated disorden^. This is relevant to 
describe the moving Bose glass state of driven vortex lattices in the presence of correlated disorder. 

Another generalization is to extend the equation (||) to a component vector Ua- It is easy to see in that case 
that a non potential non linear disorder is generated if v > (which reduces to the "static random force" for A^ = 1). 
Thus in that model it is natural to look at a generic non potential disorder F^*'^*(r, u) from the start. The mean field 
dynamical equations for large A^ and the FRG equations at any A for a large class of such models are derived in 
Section VI and in Appendix A subclass of these models is non periodic models (manifold) . They are relevant to 



describe the random manifold crossover regime in the moving glass (see below). A further subclass is then obtained 
by setting w = 0. Interestingly the resulting modeL describes polymers (and manifolds) in random flows and can be 
studied both in the large A limited and using RGeII for any A^. 

Finally, there are other simpler but interesting situations such as disorder correlated along the direction of motion 
or lattices moving in a periodic tin roof potential. These potentials which are independent of x have the interesting 
property that the steady state measure P[u{r)] is identical (at any T > 0) to the one with v = 0. 

Thus we see that the moving glass equation hides a whole class of new interesting dissipative models whith glassy 
properties. 



III. PHYSICAL RESULTS 

In this Section we present all the physical results on the Moving Glass that we have obtained in Ref. and 
in the present paper. We deliberately avoid technicalities and refer to the proper Sections for details. 



A. Channels 

One of the most striking property of moving structures described by (|^) is that the non linear static force f 
results in the pinning of the transverse displacements Uy{r, t) into preferred static configurations Uy{r) in the laboratory 
frame. Thus the resulting flow can be described in terms of static channels where the particles follow each others like 
beads on a string. In the laboratory frame these channels are determined by the static disorder and do not fluctuate 
in time. They can be visualized in simulations or experiments by simply superposing images at different times. What 
makes the problem radically new compared to conventional systems which exhibit pinning is that despite the static 
nature of these channels there is constant dissipation in the steady state. This can be seen in the moving frame where 
each particle, being tied to a given channel (which is then moving) must wiggle along y and dissipate. In fact the 
existence of the channels shows in a transparent way that the wiggling of different particles in the moving frame is 
highly correlated in space and time, thus leading to a xadically different image as the one embodied in the "shaking 
temperature" based on thermal like incoherent motiontZI. 

The channels are thus the easiest paths followed by the particles. One can see that the "cost" of deforming a 
channel along y is that dissipation is increased. Thus the channels are determined by a subtle and novel competition 
between elastic energy, disorder and dissipation. As a consequence these channels are rough. This is a crucial difference 
between what would be observed for a perfect lattice as illustrated in Figure 
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(a) v=0 

t;'-y--K--V--V--*- 

♦ ♦ ♦ » ♦ 



(c) T>0 d=3 




(b) 



T=0 



(d) T>0 d=2 



(e) T=0 Fy>0 




FIG. 8. (a) a snapshot of a perfect (non-disordered) lattice moving along the x direction, (b) upon superposing images at 
different times one would see that at T = the particles follow perfectly straight lines, (c) at < T < Tm in d = 3 the 
channels remain perfectly straight with a finite width due to the thermal fluctuation of the particles, (d) in d = 2 since thermal 
fluctuations are unbounded channels are completely blurred and cannot be defined even for T < Tm- (e) Even in situation (b), 
(c) applying a additional small force along y immediately results in tilted channels with angle Fy/Fx. 



By contrast the channels which are predicted in the Moving Glass axe illustrated in Figure || and in Figure It is 
important to stress that the Moving Glass equation (|^) introduced inE3 does not assume anything about the coupling 
of the particle in different channels but only implies that the channels themselves are elastically coupled along y, and 
thus through compression modes. Indeed on specific models such as model II one can verify explicitely that although 
coupling between longitudinal and transverse degrees of freedom exi sts a p riori, the longitudinal degrees of freedom 
Ux do not feed back at all in the moving glass equation (see section VIII A). 

The existence of channels naturally leads to several a priori possible regimes for the coupling between particles in 
different channels. The first case, represented in Fig. ^ (b), is a topologically ordered moving structure corresponding 
to full elastic coupling between particles in different channels. Since, remarkably, this structure retains perfect 
topological order despite the roughness of the channels, it is reminiscent of the properties of the static Bragg glass, 
and thus we call it a Moving Bragg glass. A second case of a Moving Glass corresponds to decoupling between the 
channels, by injections of dislocations beyond a certain leng thscal e Rd and is called the Moving Transverse Glass. 
These two regimes will be discussed in more details in section [II D . Finally note that in d = 3 channels can be either 
"sheets" (for line lattices) or linear (for point lattices) as represented in Fig. 

It is important to note .that the channels in the Moving glass are fundamentally different in nature from the one 
introduced previouslyOEj to describe slow plastic motion between pinned islands, as illustrated in Fig. (|^). In the 
Moving glass they form a manifold of almost parallel lines (or sheets for vortex lines in d = 3), elastically coupled along 
y. For that reason we call them generically "elastic channels" (whether or not they are fully coupled or decoupled) 
to distinguish them from the "plastic channels" (even though some plastic flow may occur when elastic channels 
decouple) . 

Note that in the above discussion we have concentrated on elastic channels which can spatially decouple. It is 
possible a priori that they may still remain temporally coupled, i.e synchronized. Indeed, examples of synchronization 
where observed even in extremely plastic filamentary flow. 
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FIG. 9. (a) The motion in the moving glass occurs through rough static channels. Relative deformations grow with distance 
and become of order a at distances ~ Rt^ ■ Only the channels themselves are elastically coupled along the y direction. 
Depending on the dimension, velocity and disorder strength two main cases can occur: (b) an elastic flow where the particles 
positions are elastically coupled between channels (in d = 3 and weak disorder or large velocities). In this regime the lattice is 
topologically ordered (no dislocations) and the rows of the lattice follow the channels. This is a Moving Bragg Glass, (c) In 
d = 2 or at stronger disorder in d = 3 the positions of particles in different channels may decouple. Dislocations with Burgers 
vectors along x (indicated by the square) are then injected between some channels beyond the length Ra. This situation 
describes a Moving transverse glass (with a smectic or a line crystal type of topological order). 



(a) (b) 




FIG. 10. Different types of topological order for the manifold of channels in d = 3 (a) for line lattices in motion the channels 
are "sheets" and thus form an anisotropic type of smectic layering (b) for 3d lattices of pointlike objects or for triply periodic 
structures (triple Q CDW) they have instead the topology of a line crystal 

B. Dynamical Larkin length and transverse critical force 

Another important property of the moving glass intimately related to the existence of stable channels is the existence 
of "transverse barriers" . Indeed it is natural physically that once the pattern of channels is established the system 
does not respond in the transverse direction along which it is pinned. Thus we have predicted in Ref. |6^ that the 
response to an additional small external force Fy in the direction transverse to motion vanishes at T = 0. A true 
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transverse critical force exists (and thus a transverse critical current in superconductors) for a lattice driven 
along a principal lattice direction. 

The transverse critical force is a rather subtle effect, more so than the usual longitudinal critical force. It does not 
exist for a single particule at T = moving in a short range correlated random potential. By contrast even a single 
particule experiences a non zero longitudinal critical force. It does not exist cither for a single driven vortex line or 
any manifold driven perpendicular to itself in a pointlikc disordered environment. It would exist however, even for 
a single particle if the disorder is sufhciently correlated along the direction of motion (such as a tin roof potential 
constant along x and periodic along y). Such disorders break the rotational symmetry in a drastic way. Still, in the 
case of a lattice driven in an uncorrelated potential it does seem to break the rotational symmetry of the problem. 
In some sense in the Moving glass the transverse topological order which persists (and the elasticity of the manifold) 
provide the necessary correlations (through a spontaneous breaking of rotational symmetry). Thus the transverse 
critical force is a dynamical effect due to barriers preventing the channels to reorient. 

We have investigated the equation ^numerically in d = 2 and found that indeed starting from a random configuration 
and at zero temperature the field u{r, t) relaxes towards a static configuration Ustat{f) solution of the static equation. 
Applying a small force in the y direction (i.e adding fy in P) yields no response. The manifold is indeed pinned. 

Thus we have proposed the Moving Glass as a new dynamical phase (a new RG fixed point) and the transverse 
critical force as its order parameter at T = 0. The upper critical dimension of this phase is d = 3 instead oi d — A 
for the static Bragg glass. Above d = 3 weak disorder is irrelevant and the moving glass is a mpjjing crystal. For 
d < 3 disorder is relevant in the moving crystal and leads to a breakdown of the 1/v expansion ofo. Divergences in 
perturbation theory can be treated using a renormalization group RG procedure (Section 
new fixed point which confirms the prediction that the moving glass is a new dynamical pt 



VI). One indeed finds a 
Using RG and the 



properties of this new fixed point one can compute various physical quantities (Section |VI B 3 ) . We find that the 
transverse critical force is given by: 

with c = cii in c? = 2, c = y/cuc^i in d = 3 and A is a nonuniversal constant. The length scale i?^ is the dynamical 
Larkin length. It is defined as the length scale along y at which perturbation theory breaks down, non analyticity 
appears in the FRG and the (scale dependent) mobility vanishes. Before we proceed further let us define now disorder 
strength parameters. For uncorrelated disorder the random potential V{r) which couples to the density of the structure 
has short range correlations of range r/, V(r , z)V{r' , z) — g(r — r')(5'^^(z — z') (see Section ^V] ). As in ^ we denote 
by A (also denoted A(u = 0), see Section pv|) the bare static pinning force correlator A = PoSif k =o ^y9K where 
p is the average density and gx the Fourier transform of g{r) at the reciprocal lattice vectors. Throughout A2 will 



denote the second derivative of the non linear pinning force correlator A2 = — A"(0) « A/r^ (see Section VI). Our 
result is that in d = 3 the dynamical Larkin length is given by 

Ry-ae.^'-^ (8) 



while in d < 2 it reads 



Ry^^l + 4^W3-rf))V(3-.) 



with again c — cn in d — 2 and c = y^ciiciJ in d = 3. These results are valid for large enough velocities v ^ v* (see 
below for the definition of v* and results for all velocities). Note that for v > v*c the dynamical Larkin length depends 
only on cn (and of C44 in d = 3) as it should since the physics of the moving glass is controlled by the compression 
modes and thus largely independent of the detailed behaviour along x. 

Another way to estimate the Larkin length is to compute the displacements in perturbation theory of the disorder. 
At very short distance one can treat the pinning force in (|^) to lowest order in u. This gives a model where disorder 
is described by a random force F^^'^^lx) independent of u whose correlator is (^F^^^^(^r)F^^^^{r')) = A6'^{r — r'). This 
regime is the equivalent of the short distance Larkin regime for the statics. In the Moving glass at very large velocity 
V ^ V* the displacements along y grow as B{r) — Bjip{r) (at T = 0) with: 

/• dq^dqyd'^^q^ A(l - cos(g^y)) 
""^^-^^ J i^nY {rivq.? + {cl,ql + cl,qlf ^ ^ 

The scale along y at which Uy becomes of order r/ defines the dynamical Larkin length R^, i.e Biip{y = R^, x = 0) ^ 
r J . The resulting expression coincides with the one obtained within the RG approach (up to non universal prefactors) . 
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Similarly one can define a Larkin length for transverse pinning along the x direction by the condition that Uy{x = 
Rc,y = 0) — Uy{x = 0,y = 0) Tf. Since what determines this length is only Uy (and not Ux) it is independent of 
the detailed behaviour along x. It is important to note that the Moving glass is a very anisotropic object at large 
scale with a scaling a; ~ j/^ of the internal coordinates. This implies that at large velocity [v > v*) the Larkin 
length along x is very large (much larger than Kj!), with — v{Ry)'^ /cu (One has also the more conventional 
behaviour R^r^ \JC^aJc\\R\ d — Z ). Estimating the random force acting on a Larkin volume for the transverse 
displacements^ one recovers the above estimate for . 

The resulting transverse TV characteristics at T = is depicted in Fig. ^ The transverse depinning is studied in 
Section VI and we find the behaviour near the threshold Vy ~ \F'^ — F"^^ for F"^ > F^ with ^? = 1 to lowest order in 
e. A reasonable conjecture which would be interesting to verify is that it remains 9 = 1 to all orders. Thus it starts 
linearly with a slope which depends on the velocity v. It is very large for v v* and diverges in the limit w — > 0. 




FIG. 11. Transverse v-f characteristics at T = 0: transverse velocity Vy as a function of the applied transverse force Fy at a 
fixed longitudinal velocity. The behaviour near threshold is found to be linear with a large slope for v < v^- 

The existence of a transverse critical force in a moving state raises interesting issues about history dependence. 
These issues are largely open and should be explored in further numerical, experimental and theoretical work. 

Let us for instance consider consider two experiments. In the first one a force fxBx + fy^y is applied to the lattice at 
time t — and then wait until a steady state is reached. The velocity is then {vx{f),Vy{f). In the second experiment 
one first apply a force fx along the direction x, wait for a steady state and then apply fy along y. One then measures 
the velocities {Vx if) , Vy {f) . The question is should one find the same result in the two experiments or not. Of course 
there are subtle issues which complicate the problem and needs to be further investigated (such as (i) the order of the 
limit system size versus waiting time before deciding a steady state is reached (ii) whether the lattice will globally 
rotate or break into crystallites, (iii) some non universality of T = dynamics) but one should still be able to find 
an operational answer. If it is found that there are such history dependence effects then that would be a strong 
characteristic of a glassy state (it should not happen in the liquid where one expects both answers to be the same, 
but in the same trivial sense as for a single particle). On the other hand, if no clear history dependence is found 
it has interesting consequences. We assume in the following that the global orientation of the lattice is unchanged. 
Then the first consequence is that there is a well defined history-independent global v-f function. This function 
however is non analytic in a large region of the fx, fy plane. Vy{fx,fy) should remain zero at least in the region 
fy < F^ (vxifx, fy)) and similar regions near each of the principal symmetry axis of the crystal. This is clearly the 
result of the FRG calculation presented here. But then one may also guess that it may be non analytic too along other 
lattice directions (though it is possible that some of the higher symmetry directions be screened by lower ones). The 
transverse mobility as a function of the angle and the force should exhibit a complex (and rather strange) behaviour 
which would be interesting to investigate further. A second interesting consequence would be that if in the above 
described first experiment one chooses a /j, > smaller than F^^Vx), the lattice would first glide in the direction of the 
applied force (as small time perturbation theory would indicate) but would soon change its velocity to lock it along 
a symmetry axis. It is quite possible that this locking effect exists and be a possible explanation for the behaviour 
ubiquitously observed in experiments, namely that lattices tend to flow along their principle axis directions. Such a 
behaviour near depinning was observed in recent decoration experimentsE3. 

Another important question for experiments is to determine the transverse critical force as a function of the 
longitudinal velocity v. As v decreases F^ increases but it is intuitively clear that F^ cannot become larger than the 
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longitudinal critical current (strictly speaking in the same direction y). We will neglect for now the dependence of 
the longitudinal critical current in the orientation with respect to the lattice (which gives a numerical factor which 
can be incorporated). We will call the critical current for v — 0. As d is decreased below v* the transverse 
critical force saturates at This is depicted in Fig. ^ (the large v behaviour was given in|^,^,^). There is thus a 

crossover towards the static isotropic behaviour (e.g in the Bragg glass) - assuming no dynamical phase transition as 
V decreases which would complicate the analysis. 




FIG. 12. Transverse critical force as a function of the longitudinal velocity. For a relation between tjv* and the longitudinal 
critical force see the text. 



This crossover can be explicitly estimated using the FRG in Sections Vl, VIII A and physical arguments. It is 
convenient to discuss it using the Figure ^ (also useful for studying the crossover in the correlations- see next 
Section). Let us first discuss it for simplicity with isotropic elasticity cn = cge — C44 = c. There is a crossover length 
scale Rcr = c/v below which the Moving glass looks very isotropic and very similar to the Bragg glass. This length 
scale is represented in Figure [T^ as a dashed line. Increasing the length scale R starting from a, at fixed w, one is first 
controlled by the static behaviour until reaching that line (i? < Rcr) and then one is controlled by the dynamical MG 
regime for R > Rcr- Similarly one can also represent the Larkin length at w = Rlf° = (c^r^/A)^/^''"''' (in d < 4). 
The crossover velocity v* corresponds to the velocity at which R^ = Rcr when one has also = Rl^°. One finds 
that: 



(A/cV^)ln(c2r^/aA) (d = 3) 
c(A/cV2)i/(4-<i) < 3) 



(11) 

(12) 



These results are valid when (-R*^° > a), i.e. in the collective pinning regime for the statics. We denote rj = 
min(r/, a), i.e if r/ > a one can simply replace rf by a in all the above formulae. 

Thus for V < v* the transverse critical current becomes of order the longitudinal one F^"" = cr f / {R]^°Y . It is 
useful for purpose of comparison with experiments to compare 7]v* with Fl^° . One finds the general relation: 
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(13) 



with logarithmic corrections in d = 3, i]v* = Flf°R]?° /rf\ii{R]^° /a). This result is remarkable. Since for weak 
disorder one has usually that R]^° ^ it shows that for a system with isotropic elasticity, the transverse critical 
force should remain of order the longitudinal one up until very far above the longitudinal threshold {Fx S> Fc) (very 
high up in the Vx-F^ curve in Fig. H). As we see now this is different when cee cn. 

Incorporating cgg, cn and C44 one finds that v* is indeed smaller, with: 
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(14) 
(15) 



up to a logarithmic factor ln(i?**°/a) in d = 3. we have used i?*"" = t^Cqq C44 /A in d = 3 and Rl^° = r/C66/(A) 
d = 2 and FJ:''° = cqqt f / {Rl''°)^ (we have assumed cge <C cn and neglected the contribution of compression modes). 
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Thus the value of v* is then much smaller. One sees from that a measure of the transverse critical current may 
lead to interesting information about the elasticity of the lattice. 

Finally note that one can make a simple minded argument showing directly on the equation (||) that the new 
convective term should not change pinning much at small v. Indeed starting from the case v — 0, where one has a 
pinned state u^^{r) and treating the convection term as a perturbation (which should be OK at small scales) one 
sees that this terms acts on the v — pinned state as an additional quenched random force. Since there is a critical 
force fc{v = 0) in that case, it is intuitively clear that this term will not destroy completely the state Mst^?(r) until 
VTf/Rc ~ /c- This argument gives back the correct value for w*. 




FIG. 13. Crossover as a function of the length scale R and longitudinal velocity v from the static Bragg glass behaviour 
(at small v) to the Moving glass behaviour (at large v). The dashed line represents the crossover between these two regimes 
R — Rcross = c/v. The dynamical Larkin length R^iv) as a function of v and the transverse translational order length Ra{v) 
are indicated as plain curves. This is valid in the collective pinning regime R^"" > a where _R"° is the static Larkin length. 



C. Displacements and correlation functions 



Due to the presence of the static disorder one expects unbounded growth of displacements in the Moving glass. The 
relative displacements induced by disorder in the moving system can be first computed in naive perturbation theory 
using (^. One finds: 

B{x,y)^A^H{ -) (16) 

crjv jjvy^ 

where H{0) = est and H{z) ^ z'^^^'^^l'^ at large z. Thus x scales as and the displacements are very anisotropic. 

The above formula, if taken seriously, leads to displacements growing unboundedly for (i < 3. This is similar to the 
Larkin calculation for the static problem. As in the statics it indicates that the crystal is unstable to weak disorder 
in d < 3 and that perfect TLRO is destroyed. Note that due to motion the upper critical dimension is now d = 3 
instead of d = 4 for the statics. As in the statics, the above formula and perturbation theory breaks down above W}. 
and an RG approach is absolutely necessary to compute the displacements. 

Using that the RG calculation one finds that the behaviour of displacements is controlled by a new fixed point 
characteristic of the Moving glass phase. One finds that the correlation function of displacements averaged over 
disorder can be rigorously separated into two parts: 



B(r) = {[uir) - u(0)]2) ^ B^Fir) + B^i^{r) (17) 

where -Bnl('") comes from the nonlinear part of the pinning force. While this part is dominant in the Bragg glass (and 
was computed inEj) in the Moving glass this contribution is subdominant and we will neglect it for now. The main 
contribution comes from the static random force which is generated both along y and x direction. The generation of 
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such a random force, forbidden in a static system, occurs here because of th e non p otentiaUty induced by the motion 



The complete expression of the generated random force is given in Section VIII A 



(see also Section |VI 



This random force gives a contribution to the displacement which at large scale has the same spatial dependence 
than the one naively extrapolated from Larkin regime formula (^) and thus (|l^). One thus finds: 
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(18) 
(19) 



At large scales the random force contribution to B{r) dominates. Although the formula resembles the perturbative 
one, the amplitude of the random force is given by the renormalized A^^ which has been to be extracted from the RG 
analysis and is determined by the non linear pinning force. In general A^j can be different from the perturbative A. 
In particular A^^ must vanish when u — > 0. 

A/j is a non universal quantity (contrarily to the behaviour in the Bragg glass) but one can still obtain a reliable 
estimate for Ai? by studying the crossover depicted in Fig. (|l^). If the velocity is smaller than the crossover velocity 
V* the random force will be renormalized downwards according to the behaviour in the Bragg glass phase. Thus A^ 
will be smaller than the bare A. This is illustrated in figure O 
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FIG. 14. Renormalized random force strength as a function of the velocity (left) and resulting amplitude in displacement 
correlations (see text). 



The amplitude of the displacements (e.g the prefactor of the logarithmic growth in (^6|) ) generated by the renor- 
malized random force is maximum around the velocity f * . Even at this velocity the displacements can be estimated 
as -BRF(f) ^ ^/^"-(^/^c) ^ — 3 and B-RFir) ^ r^{y/ K-i) in d = 2. At all other velocities the amplitude is much 
smaller. 

Given the form of the displacement correlation function the Moving glass will have QLRO in d = 3. One finds for 
transverse translational order correlations: 
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(20) 



In particular Akq — Tra^ Afl/(4uc). The dependence in the coordinate x is: 

X 

and thus one finds an anistropic divergence of the Bragg peaks corresponding to Kx = of the form: 

S{q) ^ ^ 



(21) 



(22) 



The question of the divergences of peaks associated to Kx > is discussed in the next section. 

In d = 2 algebraic growth of displacements imply a stretched exponential decay of CK{f) and thus that the peaks 
in the structure factor are rounded (as the dotted line in figure U) 

Note that in each configuration of the disorder the random force along y and the transverse critical force will 
compete. The physics of the moving glass will be determined by this competition. 

The roughness of the channels define an additional lengthcale at which the wandering becomes similar to the lattice 
spacing. As in the statics (Bragg glass) it is possible to estimate these lengths. At large velocity these lengths are 
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large and at these scale the system is very anisotropic. A simple argument a la Fukuyama-Lee similar to the one in 
Ref. ^ give: 



/c 



(23) 



At large v one can also obtain these lengths by looking at the displacements generated by the random force. For small 
V < V* there is a long crossover since at small scales the system looks more like the Bragg glass. As a consequence 
the estimates for Ra change. This illustrated in Fig. is determined roughly by Rcr — Ra- 

Let us summarize the main regimes as a funtion of the velocity of the moving glass, as can be seen on Fig. At 
large velocity v > v* the system is already anisotropic at the scale Rc and pinning and correlations are determined 
directly by the asymptotic moving glass behaviour. For v* < v < v* the system is isotropic at the Larkin length and 
pinning is similar to the static, but the system is still very anisotropic at scales i?^. Finally for v > v* the system is 
almost static like up to R^ and isotropic. The randon force is enormously reduced, and transverse barriers are very 
large 



Most of the novel properties of a moving structure discussed in the previous Section were obtained from the moving 
glass equation (H), which contains only the transverse displacements Uy. They thus rest only on the channel structure 
itself and not on the precise motion of the individual particles along these channels. Let us now discuss the problem 
of the coupling of particles between different channels which is important for the issue of topological, translational 
order and structure factor. 

An oustanding problem in the statics is whether or not topological defects will be generated by disorder in an elastic 
structure. Using energy arguments it was predicted that due to the periodicity a lattice is stable to dislocations at weak 
disorder in c? = 3 giving rise to the Bragg glass. The similar question of whether disorder will generate dislocations 
arise also for moving structures. At first sight the situation looks even more complicated to tackle analytically and 
furthermore precise energy arguments cannot be used because the system is out of equilibrium. However, as is 
becoming clear from the discussion in the previous Section the problem of dislocation can now be discussed here in 
term of decoupling of channels. Even in the presence of dislocations our picture of pinned channels should remain 
valid as long as periodicity along y is maintained. Before the channel structure was identified in |6^ it was unclear how 
dislocations could affect a moving structure. The existence of channels then naturally suggests a scenario by which 
dislocations will appear. In fact Ref. ^ naturally suggests that transitions from elastic to plastic flow may now be 
studied as ordering transitions in the structure of channels. Wc thus discuss the problem of the coupling of particles 
between different channels which is important for the issue of topological, translational order and the structure factor. 

The peaks at vectors with a non zero Kx will distinguish between the Moving Bragg glass and the Moving Transverse 
glass. Indeed Pj^ now determines the large scale behaviour, and thus cgg. 

Thus we call this system the Moving Bragg glass. The question of the decoupling is examined in the next section. 

Let us examine first whether dislocations will appear in the Moving Bragg glass in d = 3. The relative deformations 
due to disorder grow only logarithmically with distance, resulting in quasi long range order. At weak disorder or 
large velocity (since the relevant parameter is A/w) the prefactor of the logarithmic displacements is very small. This 
suggests, by analogy with the statics, that dislocations will not appear, leading to a stable Moving Bragg glass at 
weak disorder or large velocity. In that case the structure factor will exhibit Bragg glass type peaks (at all the small 
reciprocal lattice vectors). Note however that due to the anisotropy inherent to the motion the shape of each peak 
will be highly anisotropic the length R^ being much larger than R^. 

Upon increase of disorder the first likely transition corresponds to a decoupling of the channels, while the periodicity 
along y is maintained. This corresponds to the loss of divergent Bragg peaks at reciprocal lattice vectors with non zero 
component along the direction of motion. The peaks at reciprocal lattice vectors along y will still exlubit divergences 
(computed in the last Section). This particular case of a Moving glass was observed numericalljS and called the 
Moving Transverse Glass (see next Section). This phase has also a smectic type of order. One question is whether 
particles can hop between the channels in this phase. This however seems unlikely at zero temperature provided the 
channels are well defined. In the absence of such hops this decoupled phase can still be described by equation (|^) 
it has non zero transverse critical current. Increasing further the disorder should destroy even the channel structure 
leading to a fully plastic flow. 

An estimate of the locus of the transition between the Moving Bragg Glass and the Moving Transverse Glass is 
given by a Lindemann criterion. 



D. Decoupling of channels and dislocations 



(K(y = a) 



«.(0)]2) 




(24) 
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This indicates that edge dislocations appear first and decouphng of channels come from displacements along the 
direction of motion. 

In (i = 2 displacements grow algebraically. It is thus much more likely that dislocations will appear at large scale. 
Presumably this scale will be when the displacements will be of order a. One then easily see that dislocations will 
appear between the layers. Indeed 



(K(L?J)-«,(0)]2)^a^ (25) 

is controlled by the random force along x and by cgg (the displacements Uy are down by a factor R^/R^ - see 
below). In this regime blocks of channels of variable transverse size (depending on the strength of the disorder) will 
be separated by dislocations. 

The peaks at vectors with a non zero thus allow to distinguish between the Moving Bragg glass and the Moving 
Transverse glass. In systems with a small ratio cee/cn and stronger disorder the peaks with Kx > have a tendency 
to be smaller (and decoupling becomes easier as Ux becomes larger than Uy). Indeed the displacements at large 
scales (and thus the de cay of T LRO) are controlled by the random forces along y, Ayy and along x, Axx (they remain 



uncoupled -see Section VIII A). They act differently on Uy and Ux- Indeed, only (shear, cge) and Pyy (compression 



- , , , vv , 

cii) lead to unbounded displacements (e.g. in d — 3). The y random force will thus act mainly on Uy via compression, 

and the x random force mainly on Ux via shear. Though generally one has Axx < Ayy at weak disorder if the ratio 

cee/cii is small this could strongly favor the weakening of the Kx > peaks and channel decoupling. 

Thus, the problem of the behavior of dislocations in the moving glass system is of course still open, and constitutes 

as for the statics one of the most important issues to understand. It is noteworthy, however, that although these 

issues are of course important to obtain the structure factors and such, they do not affect the physics of the moving 

glass. 



1. the Moving Bose glass 

In presence of correlated disorder along the z direction we predict the existence of a "moving Bose glass" in d = 3. 
Indeed the same calculation as above: 



= A / '-^i^Si,.) . l!-;7^rj\ ... (26) 
(27r) ivvqxV + [ciiq^. + ci^^qiV 



now yields a fast growing displacement. Thus the disorder effects are stronger and one can expect thermal effects to 
be weaker for correlated disorder. The situation resembles the d = 2 case at T = 0. One can still predict a transverse 
critical current. Full translational LRO is unlikely but a Moving Transverse glass type of order along y is likely. This 
should be enough to guarantee, a localization effect of the layers and thus a transverse Meissner effect along y. A 
detailed study will be given irJljj-^rhe effect of correlated disorder on another dissipative glass system (nonpotential) 
were found to be quite strong inEil. 



E. Moving glass at finite temperature 

Thus the moving glasses, in its different forms, described by equation (^) is a new disordered fixed point at T = 0. 
An important question is to understand what is the effect of thermal fluctuations. 

Indeed in moving systems, as can be seen by perturbation theory, the fluctuation dissipation theorem is violated 
and a generation of temperature by motion will occur. This corresponds to the^^ysical effect of heating by motion. 
Note however that this effect is different from the "shaking temperature" effectl£3 which would occur even at T = 0. 
Indeed at T = a system in the absence of thermal fluctuations retains perfect time order which implies that no 
temperature can be generated ! 

Although the temperature will grow due to motion, this effect competes with the fact that naive power counting 
in glassy systems suggests that the temperature is an irrelevant variable flowing to zero. The competition between 
these two effects is highly non trivial and leads to new physics which needs to be investigated for the whole class of 



moving glasses of section [I B 



Remarkably, for this class of systems, new finite temperature fixed points exists. In the case of driven lattices, we 
find a fixed point at finite temperature in a d = 3 — e_expansion and d = 2 + e expansion. Similarly for randomly 
driven polymers very similar fixed points are obtainedEiJ. Thus a large class of moving glasses seems to exist at non 
zero temperature. 
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In (i = 3, the fixed point is slightly peculiar since both temperature and disorder flow to zero but can be analyzed 
along the same lines. The properties of the finite temperature phase are continously related to the T = one. 
In particular the finite temperature moving glass exhibits the same type of rough channel structure. Channels are 
slightly broadened due to bounded thermal displacements around the average channel position. Thus the asymptotic 
behavior the displacements, and structure factor, still remains similar to the ones at T = discussed in previous 
Sections. There is in addition a contribution of thermal displacements. In rf = 3 they are small, and one sees that 
the RG methods developped here allow to estimate more precisely the thermal heating effect, and to distinguish it 
clearly from the disorder effects. This will be important for determining the dynamical phase diagrams. In d ~ 2 the 
thermal displacements are large (see Section VIl). 

The main effect of temperature is to modify the v — f characteristics. One finds (Section VI) that the asymptotic 
mobility f^n is non zero. However at low temperatures or at velocities not too large the v — f characteristics remains 
highly non linear. There is still an "effective" (or apparent) transverse critical force Fy{T) as shown on figure |l5|. At 
low temperature the mobility /i^ is very small. If the velocity is not too large v < v* there are several regions in the 
v-f curve. Below the transverse pinning force slow motion due to effective barriers exist. They are a growing function 
of 1// until one reaches the finite temperature moving glass fixed point. Indeed reducing the transverse force probes 
larger and larger lengthscales. As depicted in figure O one is dominated until the scale i?cr by the Bragg glass fixed 
point for which temperature is strongly renormalized downwards. In this regime the v — f curve are nearly similar to 
the one in the static Bragg glass and thus highly non-linear. This corresponds to a creep regime. For smaller forces 
(i.e. when probing scales larger than i?cr) one crosses over to the moving bragg glass fixed point. At that point the 
FRG calculation shows that the barriers saturate. Thus below the scale /* one recovers a linear v — f charactersitics, 
with an extremely small mobility. Note that the scale /* which corresponds to the crossover scale i?cr can be much 
smaller than the critical transverse pinning force ii v < v*. 

These properties shows that even at finite temperature the moving Bragg glass remains different from a perfect 
crystal. The definition of what is "glassiness" in a moving structure is a new concept which has to be defined. In 
that respect too close analogies with the statics can be misleading. A first obvious glassy characteristics is the loss of 
translational order, contrarily to the crystal. Note however that a similar effect could be obtained by adding a random 
force by hand to a perfect crystal. However the response of such a structure to an externaly applied force would be 
identical to the one of a perfect crystal. Thus the glassy properties of the moving Bragg glass are necessarily stronger 
than such a state. The same question of history dependence as discussed above at T = can be asked. If these effects 
exist the question of the finiteness of the barriers might not be as an important issue as in Xhe static case. Note 
however that in some other examples of non potential dynamics barriers can indeed be infinitetll 

Since the finite temperature Moving Bragg glass is described by a new fixed point which still contains non-linear 
disorder A(u) the system remains obviously in a glassy regime. Some correlation functions of the system will necessary 
depend on the existence of this finite A(m). Ti4s is reminiscent of what happens in the statics where the order 
parameter is the fluctuation of the susceptibilityE3 whereas the averaged susceptibility itself remain inocuous. 

Note finally that for driven lattice (and for experimetal purposes) , the predicted existence of high (even if asymp- 
totically finite) barriers in some regime of velocities is a totally un-anticipated property of disordered moving systems. 
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F. Phase diagrams 



Having established the existence of the moving Bragg Glass in d = 3 and of the moving Transverse Glass in d = 3 
and d = 2 and discussed their properties we now indicate in which region of the phase diagram these phases are 
expected to exist. We study the phase diagram as a function of disordei\ temperature and applied force (or velocity). 

Let us first discuss the case d = 3. We have represented in Figure a schematic expected phase diagram in the 
three axis. For clarity we have not represented intermediate phases. The static phase diagram at _F = was dicussed 
in . There is a transition at finite disorder strenght between the Bragg glass to an amorphous glass where dislocations 
proliferate. Upon applying a force the bragg glass phase becomes the Moving Bragg glass in the low velocity regime 
(creep regime) and continuously extend to the Moving Bragg glass at higher drives. At weak enough disorder the 
continuity between the two phases suggests that depinning should be elastic without an intermediate plastic region. 
Upon raising the temperature the Moving Bragg glass melts to a liquid, presumably through a first order dynamical 
melting transition. 

The T = plane contains a pinned region for F < Fc{A) it is natural to expect the Bragg glass to still exist even 
for a finite force F < Fc until the depinning transition. At higher disorder dislocations appear and the Bragg glass is 
replaced by an amorphous glass. The depinning of this amorphous glass should be via a highly disordered filamentary 
plastic flow. Upon increasing the force and thus the velocity, the system should reverse back to the moving bragg 
glass. At strong disorder and finite drive the liquid extends to zero temperature. 

This different behaviors are also represented along each plane in Figure |l^,|l^,|l^, where wc have also indicated 
intermediate phases such as the Moving transverse glass. 

Determining the exact shape of the various boundaries is still an open and challenging problem, in particular in 
the square regions in Fig. 
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FIG. 16. Phase diagram in the temperature T, disorder A, applied force / variables. In d = 3 for very weak disorder since 
the moving glass is likely to be topologically ordered, the possibility of a depining without a plastic regime exists. Note that the 
Moving Bragg glass then should extend all the way down to small /. We have not represented intermediate phases (Hexatic, 
Moving Transverse Glass) for clarity. Note also that the lower plane corresponds to a small but non zero A. 
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FIG. 17. Phase diagram in the force /, disorder A in d = 3 at T = 0. The behaviour in the square region is unclear. 
An interesting possibility would be a direct depinning of an hexatic into the Moving Transverse glass, but other scenario are 
possible. 
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FIG. 18. Phase diagram in the temperature T, disorder A in d = 3 for a fixed velocity (not too small). This phase diagram 
is the dynamical version to the static one (containing the Bragg glass, the vortex glass, and the field driven transition). The 
MBG can either thermally melt (via a first order transition) or decouple because of disorder. 
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FIG. 19. Phase diagram in the force /, disorder T in d = 3. The Bragg glass phase also exist at T = for / < Fc. 



One of the strong features that emerges from this phase diagrams is the fact that the Bragg glass is able to 
survive motion by turning into the moving bragg glass. On the other hand other, more disordered phases such as the 
amorphous glass (vortex glass) are likely to be immediately destroyed at finite drive (and finite temperature) and to 
be continuously related to the liquid. 

In d—2, most of the transitions reduce to simple crossover. At F = and finite disorder dislocations are expected 
to be present. The resulting phase should thus be continuously connected to the liquid, although it can retain good 
short distance translational order. At T = there is a pinned phase until Fc, which should depin by a plastic flow. 
At larger drive disorder effects become smaller and one expects the system to revert to a moving glass state. As 
discussed earlier, due to the presence of disorder induced dislocations, this state is a Moving transvers e gla ss (if d = 2 
is above its lower critical dimension). At any finite temperature, one can use the RG flow of section VII. Since the 



temperature renormalizes above the melting temperature and disorder flows to zero the resulting phase should be a 
driven liquid. 
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FIG. 20. Phase diagram in the force /, disorder T in d = 3. There will be very long crossover not represented. 



G. comparison with numerical simulations 

Some of the predictions of the Moving Glass theory contained in the short account of our workEl have been later 
verified in several numerical simulations in d = 2 and d = 3. j— . 

The static channels were clearly observed at T = in d -j T'l They, were also-abserved in Ref . |5[ The transverse 
critical force at T = in d = 2 was observed very clearly irO and inE3 (see alsoQ) and found to be a fraction of the 
longitudinal critical force which is a reasonable order of-|magnitude. The effect of a non zero temperature T > was 
observed to weaken the effect of transverse barriers inE3 in d = 2. Note that some non linear effects were observed 
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to persist for low enough transverse force and temperature. Such an observation is in agreement in the discussion of 
Section () and can be interpreted as a long crossover. ■— ■ 

Sharp Bragg peaks were observed in the direction transverse to motion inE^ at T = 0. However the order along 
the X direction was found to have fast decay. This is consistent with a decoupling of the channels, and the resulting 
state was termed the "moving transverse glass" . This is consistent with the expectations from the theory presented 
here, as illustrated on the phase diagram (|2^). This phase is presumably the T = moving transverse glass fixed 
point analyzedJn Section () which does have a non zero transverse critical force. This is in fact confirmed by the 
observation oiEj of a smectic type of order where well separated dislocations exist between the chapaels consistent 
with the expectations discussed in Section (). The absence of long range order was also observed inEZI in a stronger 
disorder situation. r— . 

In d = 3, a simulation of a driven discrete superconducting XY gauge modeO finds not perfect but still well defined 
Bragg peaks at T > (near the melting) , a result which indicates that the driven lattice is in a quasi ordered moving 
Bragg glass state. In another study on the simpler d = 3 driven XY model at T = 0e3, it was found that indeed there 
is a phase without topological defects at large enough drive. If it carries to the lattice problem it would indicate that 
indeed there is a d = 3 moving Bragg glass state. rnn 

Finally note that there are also very recent simulations of a lattice with a periodic substrateE^El. This is a simpler 
case where a transverse critical current exist (it does exist for a single particule). It may be worthile to investigate 
this case in all details. i— . 

All the above numerically observed effects seem to be in qualitative agreement with the predictions inO. However, 
it would be very useful to be able to make a more quantitative comparison. This should now be possible, as we will 
give here more detailed predictions than the short accountEj. Among the various interetsing topics to check are the 
algebraic decay of translational order in d = 3, a detailed study of the dependence of the transverse critical force on the 
velocity, the exponent 9 of the transverse depinning, a measure of the barriers at low temperatures, a characterization 
of the history dependence, and zero and at low temperature. 



H. comparison with experiments 



The moving glass picture has also been confronted with experiments. Since these experiments need the character- 
ization of a moving structure they are challenging. The transverse critical current can in principle be observed in 
transport experiments (and will show up as hysteretic effects). These are difficult though because of dissipation in 
the longitudinal direction. i— . 

Decoration experiments on the moving vortex lattice have been performed recently by Marchevsky et alEj. In 
these experiments the external field is slowly varied and vortices are decorated while they move. The decoration 
particles thus accumulate on the regions where vortices are flowing preferentially. The lattice is observed to move 
in the symmetry axis direction and relatively large regions of highly correlated static channels are observed. These 
channels do not look like "plastic channels" but rather like the channels predicted inl£3. Note however that some 
dislocations along y appear (defects in the layered structure). This may be due to strong disorder effects or since the 
advancing front geometry is iru-the shape of a droplet some dislocations are unavoidable. Another set of experiments 
in NBSe2 was also reported inEa. Note that there has also been several decoration experiments performed just after 
the current is turned off. These can in principle probe the defect structure of the flowing lattice (though one may 
worof about transient effects) but cannot show the channel structure. 

InE3 we have suggested that the transverse bacriers may explain the anomalies recently observed in the Hall effect 
in a Wigner crystal in a constant magnetic fieldtjO. 

The qualitative analysis suggested by the moving glass theory is as follows. An electric field is applied in along 
the X direction. The Wigner crystal starts moving along x when the applied field is larger than a "longitudinal" 
threshold Ex > Ec- It produces a current along cc, Ix = qvx which is directly measured. Below the longitudinal 
threshold a highly non linear regime is observed where activated motion dominate. Since it is moving in a high 
magnetic field, the moving Wigner crystal is submitted to a transverse Lorentz force = qVxB. The geometry 
of the experiment is such, however, that no transverse motion is permitted in the stationary state (because of zero 
current boundary conditions), and thus Vy = 0. Thus the transverse Lorentz force must be balanced by a transverse 
electric field, which is thus generated, and is measured as the Hall voltage Vy. In the absence of transverse pinning 
the Hall voltage is Vy = LBIx- Remarkably, it is found in the experiment that the actual measured Hall voltage is 
indeed Vy = LBIx for small Ix, then experiences a plateau, and finally starts again growing linearly with a slope 
dVy/dIx ~ LB. We have interpreted the different behaviours upon increase of Ix as follows. For small Ix one is near 
the longitudinal depinning and it is probably a plastic flow regime with little transverse barriers. Then upon motion, 
a transition to a moving solid occurs, which is presumably a moving glass. The existence of a non zero transverse 



24 



critical force > then immediately implies that there are sliding states with Vy = as long as Fy < Fy and no 
Hall voltage necessary. r— ■ 

Finally note a recent experimentt^l on superconducting multilayers where it was iaund that the flux flow resistivity 
exhibit quasi periodic oscillations as a function of the field. This was interpretedtj in terms of dynamics matching 
of the moving vortex lattice with the periodic substrate. This is compatible with the presence of a quasi ordered 
structure in motion. 

It would be interesting to probe further the channel structure by direct imaging techniques. In particular one 
may investigate the degree of reproducibility of the channel pattern. It is predicted that upon sudden reversal of 
the velocity the channels should be different. The question of order and quasi order can be probed in experiments 
such as neutron scattering, flux lattice imaging magnetic noise experiments, NMR experiments and more indirectly in 
transport measurements. Other imaging techniques such as /i-SR NMR electron holography can also be used. Finally 
it would be interesting to check for similar effects in the presence of columnar defects since, as discussed in this paper 
we predict the formation of a moving Bose glass. 



IV. THE MODEL AND PHYSICAL CONTENT 



A. Derivation of the equations of motion 



Let us first derive the equation of motion for a lattice submitted to external force /. We work in the laboratory 
frame. This offers several advantages that will become obvious later. We denote by Ri{t) the true position of an 
individual vortex in the laboratory frame. The lattice as a whole moves with a velocity v. We thus introduce the 
displacements Ri{t) ~ + vt + Ui{t) where the i?° denote the equilibrium positions in the perfect lattice with no 
disorder. Ui represent the displacements compared to a moving perfect lattice (and corresponds to the position of the 
i-th particle in the moving frame). The definition of v imposes '^^Ui{t) = at all times. We furthermore assume 
that the motion is overdamped. The exact equation of motion can then be obtained from the Hamiltonian H by 

v^^^-^ + f-nv + Q:{t) (27) 

at ouj 



where rj is the friction coefficient and the thermal noise satisfies 0(00 (^) = 2TrjSijS{t — t'). The Hamiltonian is 
the standard Hamiltonian for periodic structure in a random potential H = Hd + -ffdis • ^oi is the standard elastic 
Hamiltonian, and i/dis describes the interaction with the random potential 

ffdis = / V{r)p{r) = Y,[ V{r)5{r - i?^ + vt + udt)) (28) 

Jr J Jr 

where the random potential has correlations {V{r)V{r')) — g[r — r') of range rj. 

In order to use the standard field description of the displacement u instead of focussing on the equation for one 
particle, one rewrites ( p7| ) as 

^dMd^^JH^^ /ay(r)<5(r-i?o + t,i + u,(i)) + /-r;^; + C(i?»(i),t) (29) 
dt oui 

In doing so one would get the same thermal noise for two particles being at the same place at the same time, instead 
of the two independent noises of equation (|2^). Since such a configuration cannot happen, going from ( p7| ) to (|2^) is 
essentially exact. ■— 1|— ■ 

As for the static caseEaES the difficulty is to take the continuum limit of (^9|) since the disorder can vary at a much 
shorter scale than the lattice spacing a. To proceed one follows the same steps than for the static case, suitably 
modified to take into account the time dependence of the displacements. One first introduces a smooth interpolating 
displacement field u{r,t) such that u{R^ + vt,t) = Ui{t) (see formula A2 of |38| , |39| ) . The field u{r,t) is the smoothest 
field interpolating between the actual positions Ui{t). All coordinates r are expressed in the laboratory frame. The 
field u{r,t), whose components we denote by Ua{r,t) thus expresses the displacement in the moving frame, as a 
function of the coordinates of the laboratory frame. 

As for the static, if one assumes the absence of dislocations at all times the particles can be labeled in a unique 
way. One then introduces the continuous labelling field: 

(l){r,t) =r-vt~u[(l>{r,t)+vt,t] (30) 



25 



Thus 4>{Ri{t),t) = 4>{R^ + + Ui{t),t) = i?° by definition, and (j) numbers the particles by their initial positions. In 
the absence of dislocations the field (p{r, t) is uniquely defined. 

To obtain the continuum limit of the equation ( p9| ) one first performs the continuum limit in the Hamiltonian as in 
|39| , to obtain for the disorder term 

iJpi„ = j dxV{x)p{x) = ~poJ dxV{x)dc,Ua+ po J dx ^ i/(a;)e*^ ('-^*-"('-'*)) (31) 



where K spans the reciprocal lattice and po is the average density. 
In (|3^) we have made the approximation 

M((/)(r,t) - M(r,i) (32) 

Such an approximation is exact up to higher powers of du, negligeable in the elastic limit, as for the static caseil. 
However the dynamic case is more subtle since such terms could generate when combined with a non-zero velocity 
relevant terms. This is the case for example of the so-called KPZ terms generated through cutoff effects. Since it 
is hopeless to try to tackle from first principle all such additional terms the only safe procedure is to assume that 
every t erm allo wed by symetry will be generated, and has to be examined. We will proceed with such a program in 
section VIII B. For the moment we only retain the dominant terms of (^l|). If one then takes the derivative with 



respect to the smooth field u{r,t) one obtains for the equation of motion in the laboratory frame 

7?9t<i -t- 771; • V<, = / $a/3(r-r')u^t, +Fp"i„(r,t)-f -77t;„+Ca (33) 

J r' 

where $Q^(r — r') is the elastic matrix. The term rjv ■ Vuq comes from the standard Euler representation when 
expressing the displacement field in the laboratory frame, —"qva is the average friction and in the continuum v is 
determined by the condition that the average of u is zero. The thermal noise satisfies in the continuum limit 



(:a{r,t)Ci5{r'.t')^2Tj^5o.p5\r-r')5{t-t') (34) 

and 

Fr{r, t) = -(5i?pin/5u,(r, t) = V{r)po ^ tK^ exp(zif ■ {r - vt - u{r, <))) - poVc.V{r) (35) 

K 

is the pinning force. Note the difference between our equation ( |3^ ) and the one derived in|6^, which does not contain 
the convective term. This difference comes simply from a different definition of the displacement fields. They consider 
displacements fields labelled by the original position of the particle (i.e. the actual position of the particle is r -t- w) 
whereas for us r denotes the actual position of the vortex considered (i.e. in the presence of an external potential V 
the potential acting on the vortex at point r is V{r) instead of V{r + u) for |65|) . 

In Appendix ^ we give a more general derivation of ( ^5[ ) valid even for cases where the equation of motion is not 
the derivative of a potential. 



B. Models and Symmetries 



Before we even attempt to solve (33), let us examine the various symmetries of the problem and define several 



models which approximate the physical problem at various levels. 

The physical symmetry of the original equation of motion (|2^) is the global inversion symmetry (r — > — r, u — > —u, 
V — + —V f — *■ — /). When the force (and thus v) is along a principal lattice direction, one has then two independent 



inversion symmetries Ix ^ {x 



X, Ux 



-V f -/) and ly ^ (y -i 



iy —> —Uy). These symmetries 
are exact and hold in all cases. They are the only symmetries of the original model (|27D. The proper continuum limit 
of equation ( p7|) must thus include all terms which are relevant and consistent with these exact symmetries. We define 
such a model as Model I, which will be studied in more details in Section ( VIII B). The additional terms can orig inate 



from e.g anharmonic elasticity, cutoff effects or higher order terms in Vw, as will be dicussed in Section VIII B 



If one drop in Model I the terms which are small in the elastic limit Vm <C 1, one obtains another model that we 
call Model II. It corresponds to the continuum limit of the equation of motion to obtain equation ( ^3| ) i.e. ( |33| ) in 
the elastic limit. Although model II is slightly simpler than model I, it only misses terms which are small in the bare 
equation but would be allowed by the above symmetries. Even if some of them are relevant, the would only be able 
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to change the physics compared to model II at very large length scales. One thus expects model II to give in practice 
an extremely accurate description of the physics. 

The Model II posesses a higher symmetry than Model I: let us examine the symmetries of the pinning force (|3^). 
Using the correlator of the random potential V, the correlator of the pinning force is: 



Since u is a smooth field it has no rapidly oscillating components and thus in ( |3^ ) the terms that are rapidly oscillating 
in r + r' can be discarded. Setting K' = —K in (pq), one is left with 

= pIY^ K^Kpg{r - r') exp(iX • (r - /) - iK ■ {urt - Ur't' + v{t - t')) (37) 

The symmetries of ( p^ ) thus a priori depends on the precise form of the correlator g{r). However in the elastic limit 
it is legitimate to replace Ur't' by Urt' in the above expression. Integrating then over r' one obtains: 

Aq/3 = Po X! KaKpgK exp{~iK ■ {urt - Urt' + v{t - t')) (38) 

where is the Fourier coefhcient of the correlator g{r). Since gx is essentially zero for K 3> l/?"/, the error made 
in the above approximation is itself of order Vu and thus consistent with the elastic limit approximation. This 
justifies the choice of ( p8| ) as the random force correlator in Model II. The disorder term then posseses the statistical 
tilt symmetry (STS) Urt — > Urt + /(?") where /(r) is an arbitrary function. In this case one can absorb any static 
change in u without affecting the correlations of the pinning force. Furthermore in the case of isotropic elasticity, the 
additional inversion symmetry y —y holds. i— . 
Though we will study the complete Model II in Sections ^ and VIII B its main physics can be understoodl£3 by 



noticing that the pinning force F^™{r, t) in ( |35| ) naturally splits into a static and a time-dependent part: 

i^^'^'(r,u) = y(r)po *A'„exp(iX-(r-w))-poVal/(r) (39) 

K.v=Q 

F^^''{r,t,u)^V{r)pQ ^ iK^ei^^{iK ■ {r - vt - u)) 

The static part of the random force comes from the modes such that K.v — which exist for any direction of 
the velocity commensurate with the lattice. The maximum effect is obtained for v parallel to one principal lattice 
direction, the situation we study now. This force originates only from the periodicity along y and the uniform density 
modes along x, i.e the smectic-like modes. Since this static pinning force F^^^(r, u) is along the y direction, it is useful 
consider only the transverse part (along y) of the equation of motion ( p3| ) dropping F^y^. This leads to introduce 
Model III, defined by the following equation of motion in the laboratory frame: 

ridtUy + rjvd^Uy = cV^Uy + F"*^'(r, Uy{r, t)) + Cy{r, t) 

F''-\x,y,Uy) = V{x,y)po ^ KySinKy{uy - y) - podyV[T) (40) 

Although for non isotropic elasticity Ux also appears in the equation of motion for Uy, and can in principle lead to 
additional static effects, it is not included in Model III for reasons which will be explained below. Thus Model III 
only involves the transverse displacements Uy . It posesses the same symmetries as Model II with the three additional 
independent symmetries y — > — y, Uy — > —Uy and [x — + ~x, v — > ~v) and is also defined in the elastic limit. 

It is to be emphasized that although the derivation of model III was given here systematically starting fcom an 
elastic description, the only serious hypothesis behind Model III is the existence of transverse periodicityO. As 



discussed in section flB\ ( [1 01) will be the correct starting point to describe any kind of structure have such transverse 
periodicity properties. One thus expects model III to be the generic equation containing the physics necessary to 
describe these structures. 



V. PERTURBATION THEORY FOR THE COMPLETE TIME-DEPENDENT EQUATION 



Let us start by a simple perturbation analysis of the equation of motion Moc 
disorder expansion has a long history in various contexts such as vortex latticesE^'^ 



. Such a large velocity or 
and charge density waves 



27 



The natural idea is that at large velocity the disorder term oscillates rapidly and averages to a small value and that \/v 
is a good expansion parameter. As we will see such an idea is in fact incorrect, since previously unnoticed divergences 
appear in the perturbation theory. 



A. Analysis to first order 

We start from the initial equation ( ^ ) defining Model II that we rewrite as: 

iR~^)rtr't''^r't' = f» ^ Va0Vp + fa{r,t,Urt) (41) 

where from now on we drop the pin subscript on /. The response kernel R is defined in Fourier space: 

{R-XL-f = ^tr5,'^-,{Vc.pdt + iri^v^q^So^p + CT{q)Plp{q) + (42) 

where the elastic matrix is : Criq) = c^^q^, CL{q) = Cnq^ for a two dimensional problem and Cxiq) — c^^q^ + cuAq^, 
Criq) = c\\q^ + 044^2 for a three dimensional problem. The bare value of the friction coefficient 770 is defined as 
Tiap = SaprjQ. In (|4^) the velocity is fixed by the constraint that Ui^) = to all orders in perturbation theory. This 
is equivalent to enforce that the linear term in the effective actiorO is exactly zero. 

Instead|-Qf working directly with the equation of motion it is more convenient to use the Martin-Siggia-Rose 
formalismEZl. The generic de Dominicis-Janssen MSR functional is given by 

Z[h, h]= [ i:)uL>ue"'^["'"l+''"+*''" (43) 



where h, h are source fields. The MSR action corresponding to the equation of motion ( p] ) and the disorder correlator 
(§) is 



with 



S[u,u] = So[u,u] + Sint[u,u\ (44) 



So[u,u] = / iKt(.R-'):t'X't' - - - r^T / («<,)(*<,) (45) 

J rtr't' J r,t 

S,nt[u,u] = j drdtdt'{iu'^t){lU%)^.'^^{Urt-Urt' +V{t-t')) (46) 



Note that ( ^4| ) corresponds to the action derived in Appendix 



The fundamental functions to compute are the disorder averaged displacements correlation function C"^'^,^, 



and the response function R"^^^,^, — 5 {u'^^) / Sh^^, which measures the linear response to a perturbation 

applied at a previous time. They are obtained from the above functional as C^f^,^, — {u"fU^n,)s and R'^f^,^, — 

{u'^^ivf^ni) s respectively. Causality imposes that Rrt,r't' = for t' > t and we use the Ito prescription for time 
discretization which imposes that Rrt,r't — 0. We assume here time and space translational invariance and denote 
indifferently Crt,r't' — Cr-r',t-t' and Rrt,r't' — Rr-r',t-t' by the same symbol, as well as their Fourier transforms 
when no confusion is possible. Note that in this problem C_,-^t ^ Cr.t when v is non zero. In the absence of disorder 
the action is simply quadratic S = Sq. The response and correlation fucntion in the absence of disorder are thus (for 
t > 0) and introducing the mobility /i = l/ry: 



Cq,t = Pi)3(9)--^e-('=^('?)''l*l+™'?="'^*) + Pj^(g)— y-^e-('=^(*)^l*l+™'==^*) (48) 



R-^^l = Pi'^(g)/ie-("^(9)+™«-)^*0(t) + Pj^((7)/ie-("^(«'+*"'?.)Mt0(i) (47) 
T T 
'cUq) ^ '^-^'^'^^ CT{qy 

Note that the fluctuation dissipation theorem (FDT) 

TP^^^ = -d{t)dtC^^^ (49) 

does not hold here (it holds only for v — 0) since we are studying a moving system which does not derive from a 
Hamiltonian. 
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It is easy to show that the disorder does not produce any correction to the part iut{cq^ + ivqx)ut of the action, and 
thus that the parameters c (cn and cge ) and tjqv are not renormalized (we consider here for simpHcity the isotropic 
version c — cn — cgg but this property holds in general). This is similar to the property of non renormalization of 
connected correlations in the statics (for v — 0) due to the statistical tilt symmetry. Indeed here one has the exact 
relation: 

\ogZ[ht,h] = - j dt ht{cq'^ +ivq^y^h + \ogZ[ht,0] (50) 

where h is an arbitrary static field. This relation implies that the static response function J dt' Rq t.t' — (^9^ + *i"?x)^^ 
is not renormalized. The property (M) is demonstrated by performing the change of variable ut^r = ,. + fr with 



fq = —(cq'^ + ivqx)~^hq in the action and noticing that the change in the random force F{r, u) — s- F'{r, u) = F{r, u + f) 
is thus that all averages over F and F' are identical since they have the same correlator from (|3^). 

Let us now study the perturbation theory in the disorder and compute the effective action r[u,u] to lowest order 
in the interacting part Smti using a standard cumulant expansion 

T[u,u] = So[u,u] + {S.m[u + Su,u + Su])su,Su (51) 



where the averages in (Bl) over Su, Su are taken with respect to Sq. 

The calculations are performed in Appendix One finds that the effective action has the same form as the bare 
action, up to irrelevant higher order derivative terms, with the following modifications. First the full non linear form 
of the correlator of the pinning force is corrected by thermal fluctuations A'^ ^"if ■ In d > 2 it reads: 

A'^^ = e-*^'^~ (52) 

or equivalently gK — gxe"^^^^'^ where Boo = {u'^)th- Wc have defined i3"f — 2(C^g — C"f ). This amounts to a 
smoothing out of the disorder by thermal fluctuations. Secondly, the friction coefficient matrix is corrected by Srjap, 
and the temperature by ST. Finally, the driving force is corrected by Sf (we are working at fixed velocity, enforcing 
order by order that f + Sf ^ rjv). 

Let us start by the corrections to the driving force. We find: 



K I=L,T 



To derive this formula from (??) we have used the symmetry K —K. This formula gives the lowest correction to 
the driving force at fixed velpeity or, equivalently to the velocity at fixed driving force. It is identical to the formula 
(22) of Schmidt and IIaugerE3. There are small differences, unimportant in the elasti c limit, which come from the 



different definitions of the continuum limit of the model (see discussion in section |VA|) . 

A salient feature of the above formula was noticed by Schmidt and Hauger, i.e the velocity and the force are not, 
in general, aligned. There are aligned however when the velocity is along one of the principal lattice directions, i.e 
Kq.v = where Kq is one of the principal reciprocal lattice vector (note that this is also the case for the median 
direction 7r/6). Such a feature is reasonable on physical grounds and can be confirmed by higher order analysis of the 



perturbation theory (see section VI) 



Furthermore, using the approximation v{K + 9) ~ vK Schmidt and Hauger found that, in d = 2, the transverse 
pinning force versus the angle a between the velocity and one principal direction of the lattice has a discontinuity 
at a = 0. One could naively think that such a discontinuity could be interpreted as the existence of a transverse 
critical current. Indeed a natural interpretation of figure 1 of ^ would be that one need to apply a finite force to the 
lattice (opposite to the transverse pinning force) to tilt sighthLits velocity from the principal axis direction. Notable 
confusion on this subject existed at that time in the litteratura£3^L3. Such an interpretation is in fact incorrect. First, as 
Schmidt and Hauger correctly pointed out such a discontinuity is an artefact of the approximation v{K + q) ^ vK, 
and diseappears if the correct expression (^3|) is used. Furthermore it is easy to check that even with the above 
approximation the discontinuity exists only in d = 2 and the function is continuous for d > 2. Thus the first order 
perturbation does not exhibit any divergence and does not give rise to a transverse critical current. 

In order to have divergences in the perturbation theory (and the associated effects) it is thus necessary to examine 



the perturbation theory to second. We will perform such a calculation in section VI 



Before we do so it is interesting to examine the first order corrections to the friction coefhcient and the temperature. 
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At zero temperature it reads: 



(5ry"'3(t;) = _ y drdrrW^r = 0, t)A°'^'^^{vt) (55) 
Still at T = and using the bare form of the disorder one finds: 

'5'?"/3 = E E / dqK^Kp{K.P\q).K)gK , (56) 

When the velocity is along a principal lattice direction one finds that Srjxy = and thus the friction matrix remains 
diagonal. However the corrections to the friction is clearly not the same along x and y. 
Next we give the corrections to temperature. They are: 

5{^T)^p = ^ E / A^^e-'^-"*(e-^^-^°.-^ - e^^^-^^-^) (57) 
2 Jt 

Contrarily ot the velocity corrections, corrections to the temperature (|5^) exhibit divergences for any v already at the 
first order in the disorder. However these divergences are weilJiiddcn and can only appear if one looks at the non 
zero temperature perturbation theory, which was not done inE3E3. At T = one finds trivially that 5T — 0, showing 
that disorder alone cannot generate a finite temperature. Such corrections are thus non trivial and there is in general 
no simple relation between S(riT)ap and 6{vi)ai3, due to the absence of FDT theorem. Only in the particular case 
where v = and of potential random forces A^* = KaKpgK the FDT theorem enforces 5T = Q (see, e.g ( ^ ) andEj). 
The way to treat these divergences will be examined together with the second order in perturbation in section VI. 



B. Correlation functions 



The last physical information that can be extracted from the perturbation theory is about the correlation functions. 
The calculation is performed in appendix 1^, and this gives (see (??)) in the static limit (see (??) for the full expresion) 



E 

K.K.v=a I=L.T.I'=L,T 



E 



dq{2TT)gK—r^——. . . 

BZ ci[q)+ir]oV.qCj{q)-i'noV.q 



D 
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where the disorder is renormahzed by the temperature 



(58) 



(59) 



Let us first focuss on the T = limit of (^8|). The length i?^'^ can be estimated by looking at ^ = {u^^Ug + u^gUq) 
This gives 



E 

K,K.v=0 



BZ 



dq{2^)Kl 



ql 



qliv^ql + (C66(g2 + gg) + c^iqlf) qliv^ql + {cii{ql + (?2) + a^qiy 



(60) 



where q\ — q^ + qy. Since qx q^ one immediately sees from ( pCj ) that the compression modes are the one responsible 
for the divergences. Other modes are not divergent around d = 3. Lengthscqles qre theus controled by cn. 
On the other hand in order to obtain a decoupling of the channels one can use a simple Lindcrmann criterion. 



{[ux{y = a) - UxiO)Y^ 



Using 



one gets 



Bxxiy 



A 



vv 



q.BZ 



q\{v'^ql + (C66(g^ + ql) + Ci^qD) 



(61) 



(62) 



we have neglected all terms containing cn. One can thus neglect the compression mode. The decoupling between 
the channels is thus controled by cee whereas the roughness of the channels and the characteristics lengthcales of the 
moving glass will directly depend on the compression mode cn. Estimating the integral one finds in rf = 3: 
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BUV = «) - ^ min(( — {a/vf/') (63) 
One recovers the Bragg glass estimate (in the simpler case a — rj). 

The above perturbation expansion for the Lindemann criterion implicitely supposes that the random force is directed 



along the y direction. In fact, under renormalization a random force along x is also generated as discussed in section II]. 
One can use again the Lindemann criterion with a random force along x. In that case only the transverse modes 
remain are important (fo small enough qx). 
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VI. RENORMALIZATION GROUP STUDY OF THE TRANSVERSE PHYSICS (MOVING GLASS) 



Up know we have studied the perturbation expansion of the full continuous model II, keeping both the x and y 
directions. Doing the second order perturbation on that full model is tedious. Since one knows on physipal grounds 
that the singularities in the perturbation theory will come from the static components of the disorderES , which as 
was discussed in Section 11 B originate from the transverse degrees of freedom, we now study the perturbation theory 
of the simplified transverse equation of motion^ Model III. If, as we indeed find, this perturbation theory is singular, 
this implies divergences in the full model II as well. Model III will thus g ive a g ood description of the physics. We 



thus study it here and come back to the full Model II in the next Section VIII A 



A. zero temperature perturbation theory to second order 



To avoid cumbersome expressions, and since in this whole section we only discuss transverse degrees of freedom we 
skip the index y for Uy . We also discuss here for simplicity a, n = 1 component model for the transverse displacement 
u (which is appropriate for flux lines in d = 3 and point vortices in d 2). Generalizations to n > 1 will be briefly 



mentionned in Section VI D 



We thus study the dynamical equation for Model IIIE3. Since we will be dealing with an anisotropic fixed point it 
is useful to distinguish Cx and Cy. 

rjdturt = (c^V^ + cy^l - vdx)urt + F{r, Urt) + Cir, t) (64) 

The bare value of the friction coefficient along y is denoted by ?7o, and for simplicity we denote by v the quantity t^qv 
(which remains uncorrected to all orders in this model). The correlator of the static transverse pinning force is (EO): 



F(r, u)F{r', u') = A(w - u')5'^{r - r') (65) 
Averages over solutions of (^4|) can be performed using the Martin-Siggia-Rose (MSR) action: 



5*0 [u, m] = / iUrtivdt + Vdx - ^~ Cy^ y)Urt ~ r]T{iUrt){iUrt) (67) 

Jrt 

Sint^-\( {iUrt){iUrt')^{Urt- Urt') (68) 
^ JrW 

In this section we will restrict ourselves to T = 0. The corrections coming from correlation functions {Sutdut') ^ T 
then vanish, which simplify the analysis. This can be used to show that to all orders the temperature remains zero. 
The first order corrections where computed in the Section ??. At T = there is no correction to the disorder term (to 
this order). There is a non trivial correction to the kinetic term, which gives the following correction to the friction 
coefficient ry: 

/< i<+oo 

Sr] = -A"{0) / / dttR{q,t) (69) 

Jq Jo 

This leads to: 

^ ^ _A"(0) / '^'"'g^ (70) 

r, ^ V (27r)<'-i(4c,Cj,92 + ^2)3/2 ^ > 

(check) This is not a divergent integral, except when v ~ 0. To find divergences in the perturbation theory at T = 
one has to go to second order. 

The second order corrections to the effective MSR action, and thus to the coarse grained equation of motion, are 
computed in Appendix |^. To second order a correction to the full nonlinear disorder correlator appears and reads: 

dA{u) = A"(ii)(A(0) - A(w)) f G{r)G{r) - A'{uf j G{r)G{-r) (71) 

where G(r) is the static response function 
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CW^/ , Cl,)^-^—^^ (72) 



At zero velocity both terms in ( ]71| ) are infrared divergent for d < 4, as is well known leading to the glassy effects in 
the statics. The key novelty with respect to the problem at f = is that due to the assymetry introduced by motion, 
G(r) is different from G(— r). As a consequence the second term in ( [7l| ) is now convergent for v > 0. Indeed the 
integral 

G{r)G{-r) ^ f G{qf = f ^f^— r (73) 

is convergent in all dimensions for v > 0. On the other hand, one divergence remains from the first term: 

G{r)G{r) = / G{q)G{-q) = / / - / ^Jw-^. (74) 



Q 



i^^r-'2cyql^Ac.Cyql+v^ J C^nr-'^vcyq^ 



The integral (^4|) is divergent for d < 3, even for v > 0. Thus, contrarily to general belieS~il originating mainly form 
the study of the first order perturbation in 1/v, analysis to second ord^K. confirms the surprising conclusion that even 
at large velocity infrared divergences occur in the perturbation theoryE30. Such divergences indicates the instability 
of the zero disorder fixed point and the breakdown of the large v expansion. They lead the system to a novel fixed 
point where the disorder is plays a crucial role. The above divergence is the key to the physics of the moving glass. 



B. Renormalization group study at zero temperature 

In order to handle these new divergences, and to find the new fixed point which describes the large scale physics, we 
use a dynamical functional renormalization group (DFRG) procedure on the effective action using a Wilson scheme. 
This allows to keep track of the full function A(u), which is necessary since the full function is marginal at the 
upper critical dimension. This is equivalent to decompose the fields into fast and slow components u u + Su and 
it —> u + Su and to integrate the fast fields Su and Su over a momentum shell. This method is very similar to the 
method introduced in Ref. p5|p4| for the v — case, though differs in details. 



1. derivation of tfie RG equations 



The first task is to perform a dimensional analysis of the MSR action and to determine the appropriate rescaling 
transformation. Since we want to describe both the v — and v > fixed points, we perform the following redefinition 
of space, time and the fields (keeping also an arbitrary T) : 



y = y e x = x e t = t e (75) 



We now impose that the action S in ( pq ) is unchanged, which yields redefinitions of the coefficients. Since Cy 
since this quantity remains uncorrected to all orders, this fixes a — i — d— a — z — C,. One finds the rescaling: 
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T^-^n' ^ r;e(2-^)' v^v' ^ t>e(2-'^)' (76) 

A ^ A' = Ae(5-'i-'T-2C); (73) 

In the case w = the natural choice is cr = 1, which yields A' = Ae*^"'^''^^^^'. Power counting at the gaussian fixed 
point {z — 2^ C, — Q) yields iiie-,upper critical dimension due — 4 below which disorder is relevant (and a c? = 4 — e 
expansion can be performedE3^E3. 

For D > since v is uncorrected to all orders a natural choice is cr = 2. Power counting near the Gaussian fixed 
point (z = 2, C = 0) indicates that now the upper critical dimension is thus d^c = 3 ( with A — > A' = Ae^^"''"^^'' ). 
As a consequence disorder terms are relevant for dimensions c? < 3, whereas the temperature appears to be formally 



irrelevant (see however Section VIC ). The elasticity term along x {cx) now corresponds to an irrelevant operator 



at the anisotropic fixed point. Note that the above rescaling (iTq) indicates that the proper dimensionless disorder 
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parameter is A/wA'^"'^. Here we are mostly interested in periodic systems for which, as in the staticsE^, one must set 
C = 0. For completeness we give however the equations for non periodic systems (C > 0). 

The standard RG method consist in two steps. Firstly, one integrates the modes between a < y < ae' or equivalently 
Aq > Qy > Aoe^' with Ao ^ 7r/a, which yields corrections to the bare quantities. The cutoff procedure we choose here 
for convenience is to integrate over the following momentum shell: 

^ d'i-^qy [+^ dq.. 



dci= 7^ ^ (79) 

This results in the same theory but with a different cutoff and corrected parameter. Secondly,one perform the length, 
time, and field rescaling (|75|), as well as the corresponding change of quantities ([75|), so as to leave the effective action 
invariant. The cutoff has thus been brought back to its original value. Scale invariant theories will thus correspond 
to fixed points of this combined transformation. 

The RG equation for the disorder can be also established. The shell contribution of the integral (jTJ) is asymptoti- 
cally: 



«GW.-GW~//%4^~5^-4,_,A-3 (80) 



where = and in d = 3 we need A2 = l/(27r). Using ( |7l| , |7^ , p0| ) we obtain after rescaling we obtain the 

following FRG equation for the disorder correlator: 

^ = (3 - d - 2C)A(.) + C.A'(.) - A'iuf^ + (81) 

^ A"(u)(A(0)- A(u)) (82) 



'Invcy v^l + 4c^(?)cyA^/i;2 
where Cx{l) — CxC^'^K In the large scale limit it reduces toEBS 
dA{u) 



{3 ~d~2C)A{u) + CuA'iu) + - A"(m)(A(0) - A(w)) (83) 

dl AiTVCy 

The equation |^ allows in principle to examine the intermediate scales crossover when v is not very large. Indeed 
there is a characteristic crossover length scale: 

Lcross = 2^CxCy/[riov) (84) 

such that ( ^ ) becomes valid for e' ^ LcrossAjj-^ote that setting w = in the above equation ( ^l|) leads back the 
FRG equation for the usual manifold depinningSj^a (up to numerical factors originating from choice of short distance 
cutoff, and different choices for rescalings). 

We also give the RG equation for the friction coefhcient for the case of the periodic problem (C =0). Before rescaling 
it reads: 

^ = -^"'°'-^-- (4e.i:AC)./^ <'*^' 

Taking into account that A — Aoe^' and Cx{l) — Cxe~^\ and AJ'(O) = A"(0)e'^"'''' one find the RG equation, after 
rescaling: 

. A"mU 2c..(/)Ag-^ 
thus except for u = 0, 77 is corrected only by a finite amount (as long as A"(0) is finite, see below). 
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2. study of the FRG equation 

We now study the FRG equation (|8^) for the perfodic problem (C = 0). Thus we impose A(u) to be periodic of 
period 1 and study the interval [0, 1]. Let us look for a perturbative fixed point in d = 3 — e. Absorbing the factor 
^^J^ ^ in A(m) and redefining temporarily el — > I. the FRG equation reads: 

^ = A(u) + A"(u)(A(0) - ^{u)) (87) 

No continuous solutions such that dA{u)/dl — exisJlS^. This is due to the fact that the average value of A(m) on 
the interval [0,1] must increase unboundedly. Indeed integrating inside the interval one finds: 



A(^.) = / A{u) + / A'iuf + [A'H(A(0) - A{u))]i; = / A(u) + / A'^ 



d 
Jl 

It is thus natural to define A{u) — A(0) — A{u) which satisfies: 

^=A{u){l + A{u)) (89) 

Note that physically one expects A{u) > 0. This equation has a fixed point A (u) — u{l — u)/2. It is shown in 
Appendix ^that this fixed point is stable (locally attractive). Since the flow equation for A(0)(Z) is simply: 

^ = A(o)(o m 

the value of A(0)(/) thus grows unboundedly as A(0)(/) = A(Q)^p-|(restoring the e factor). 
Thus the full fixed point solution in a d = 3 — e expansion isCjS: 

Ai{u) ^ A*{u)^ A*{0) + Ce'^ A*{u)" = l (91) 

where C is an arbitrary constant and: 

A* {u) = C* + ieATTvCy){^u^ - ^u) (0 < u < 1) (92) 



and the solution repeats periodically as shown in Fig. 21 We have restored the factor l/{4:T:vCy) and e = 3 — d. 

In X-space the fixed point solution can be written Aa- = ioi K {K ~ 27rfc with k integers) and 

Ak=o{1) = Ai{u = 0) + (1/12) = Ao(u = 0)e' + (1/12). 



A(u) 



t) T ^ ^ " 

FIG. 21. Solution of the FRG equation. Note the non analyticity at all integers 
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Thus there is an ever growing average to the correlator. Remarkably, does not spoil the above fixed point, since one 
can always separates A(0) and A(m) — A(0) in the starting MSG action. In perturbation theory one sees that A(0) 
has no feedback at all into the non linear part. It simply means that there is an unrenormalized random force which 
simply adds to a non linear pinning force, which is described by A*{u). 

Note that this solution has cusp non analyticity at all integer it. This was to be expected since one has exactly: 

= A"(0)- A"(0)2 (93) 

dl 

with C = 1 in the asymptotic regime. At the start A"(0) is negative (since A(w) is an analytic function with a 
maximum at u = 0). One easily sees that A"(0) becomes infinite at a finite length scale (interpreted as the Larkin 
length, see next section), and the function becomes non analytic (also A"(0+) becomes positive). 

Once the solution is known in d = 3 — e it is straightforward to obtain it in the physically relevant dimension d — 3. 
In d = 3: 

^ = A"(u)(A(0)-A(tt)) (94) 
Defining A/(u) = jA/(?i) and introducing /' = In^ one finds that: 

l'^^ = A(^) + A"{u){AiO) - Hu)) (95) 

which is identical to (|7|). Thus the physics will be controlled by the slow decrease to zero of disorder at large scale 
with the following stable fixed point behaviour: 

AKw)~ A(0) + ^p(u2-ii) (96) 

( note that the random force term does not grow by rescaling in d = 3). 

Finally, note that the equation (^3|) presents several differences and some remarkable similarities with the one 
describing the statics FRG and the dynamical FRG for i; = 0ina(i = 4 — e expansion. The statics FRG equation is: 

^ = R{u) + hi^uf - R'\u)R"{0) (97) 

where R(u) is the correlator of the random potential and the dynamic FRG equation for (i; = 0) is 

- A{u) + A"(u)(A(0) - A{u)) - A'(m)2 (98) 

Since one has A{u) = —R"{u) it yields the solution periodic in [0,l]ii: 

A*(m) = — (1 -6u + 6u2) (99) 
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Remarkably both the solution for v ~ and for w > are non analytic at integer u, though the detailed form of 
these solutions is different. Since this non analyticity is related to glassiness and pinning one can expect a certain 
continuity of properties between the moving and non moving case. 

The main difference however is that in (9^) A(0) = A(0''") starts growing at initial RG stages as for d > but is 
stopped to its fixed point value ^ beyond the scale at which a nonanalycity develops (Larkin length). This effect 
is due to the term — A'(0+)^ and physically means that in the case v = displacements grow much more slowly at 
larger scales. The system remembeiii^that it is a potential system and thus J A{u) remains zero if it is zero at the 
start (at least formally, see howeverca). By contrast for the moving system one has asymptotically A/(0) ~ Aooc"^'. 
As will be discussed later this corresponds to the generation of a random force (which cannot exist in the statics). 

Note that for v not very large one can see on equation |l]) that there will be a long crossover during which the term 
— A'(0+)^ will act (random manifold regime, see below). Thus the actual value of Aqo should be decreased compared 
to the value naively suggested by perturbation theory, an effect studied in the next Section. 
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3. Physical results at T — 



We now extracts some of the physics of the moving glass from the FRG analysis. ^From the equation for the second 
derivative of the force correlator: 

dA"(0) 



dl 



with l/C{l) = Anvcy^yi + Ac^e 
infinite. We first estimate it in the large velocity regime L 
one has 



(3-d)A"(0)-C(0A"(0)2 (100) 
^'cyAp/iJ^ it is possible to extract the length scale Kj! at which A"(0) becomes 



<C a where one can set C{1) = l/{AiTVCy). In d = 3 



a;'(o) = - 



(101) 



where A"(0) = — A2 is the bare value. Thus: 



(102) 



This length scale, introduced irj£3'E3, is the analogous of the Larkin length for the statics. Indeed Kj! coincide with 
the scale at which the scale dependent mobility fJ.{L) vanishes as depicted in Fig. ^ This can be seen from (|8^). 
The divergence of A;'(0) drives at L = eJ = Ry drives fi{L) ^ l/-q{L) to zero for all larger scales. Beyond that scale 
pinning starts to play a role. 









\ ^ 




► 



FIG. 22. Scale dependent mobility. It vanishes beyond the dynamical Larkin length it vanishes (at T = 0). 



In c? < 3 one has: 



Thus the dynamical Larkin length is: 



'(0) 



Aae^ 



A2(e 



47rf Cyg 



Rl = {l 



(103) 



(104) 



Note that it is the second derivative A2 of the force correlator which appears in the Larkin length. For a realistic 
disorder with a correlation length r/ one has A2 = A(0)/r^. Using this relation, one checks that (104) is the one 
obtained inil by estimating the length scale at which Udis ^ ff. 

One can also determine the dynamical Larkin length when the velocity is not very large. Restoring the proper 
dependence of C{1) in I in (100) one gets that I/A2 = /q" e'^^C{l) with Ic = ln(i?^/a). This yields after some algebra 
in d = 3: 
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47r77ovcy 1 -S7r?7Qi;cy 

RV = e-^^-{l + e ^2 + (1 - e ^2 ' )^/ 1 + 7"7;;" ) (105) 



and in d = 2 one finds: 



A2 A2 y (ryot^)^ 

where we recall Aq ~ tt/o. These formulae int erpo late smoothly between the Bragg glass and Moving glass results. 

Finally, note that since the above equations 103| are exact (within the 3 — e FRG approach) the calculations of 
the Larkin lengths are independent of whether there is an intermediate random manifold regime, i.e it holds both for 
a > rf and a < Vf. Non universal irrelevant operators will of course change the numerical values of the prefactors 
but the above expressions should be correct when all the Larkin lengths are large. j— . 

One of the remarkable properties of the moving state is the existence of transverse pinning23. This is demonstrated 
from the FRG fixed point, due to the non analyticity of the fixed point (p^). Adding an exter nal fo rce fy along y (i.e 



in the l.h.s of (|64|)) generates a velocity Vy. The naive perturbation theory result , formula ( |B37| ), for Sfy{vy) (the 
correction to the applied force at fixed Vy) reads: 

SfiVy)^ J^Rr=0,t^'{Vyt) (107) 

In the limit of vanishingly small Vy one gets a non zero limit 6fy{0^) ~ -^F^, i.e a transverse critical force only if 
the function is nonanalytic with A'(0+) < 0. The critical force is thus given by summing up the contributions of all 
the successive shells 

« - / d/A;(0+)A,_i y ' e-^^--')' (108) 

where quite logically only scales beyond the Larkin length give a non vanishing contribution. Using the asymptotic 
value for A'(0+)* = e'in'qQVCy one finds: 



where C = Arf_i(3 - d)'KK'^-^/a. 
In d = 3 one finds: 



« CcyaiRl)-"^ (109) 



with C" = Ad-iTiP^ V«- 

Note that the prefactors are not universal (and in d = 3 affected by the highr orders). In the above formulae we 
have assumed a direct passage from the Larkin scale regime to the asymptotic periodic regime, thus r/~a. Ifr'/<Ca 
an intermediate random manifold regime will be first reached were the typ ical value for A;(0^) will be rather a/rf. 
This will yield to the replacement of a by r/ in the nuHierators of ( |l09| , [ll0| ). 

Remarkably this coincide with the estimate given inl£a obtained by balancing forces. Note that this result can be 
obtained here without any reference to a Larkin length along the x direction. This illustrate that the physics of the 
moving glass depends only on the periodicity along the y direction. 

Note that a more general expression of Sfy{vy){l) can be given and reads: 

SUvym= f dl I i(-^)i/2e-Mt(c„e--A^+.V(4c.))Ad-ie-('i-i)'A;(«,„i)e-^' (111) 

2 Crirt 



t>0 



It is interesting also to study the behaviour near the transverse depinning threshold. It is easy to see that because of 
the absence of IR divergence in the integral for 77 the exponent at the threshold remains uncorrected, i.e Vy ~ {fy — F-iY 
with = \ (to first order in e). The slope can easily be estimated from A"(0) the above results and becomes large at 
small velocities. 

We now study the displacement correlations. The growth of the average of A(u) implies that there is a static 
random force generated. However, unlike in the v = case, the critical force does not kill the random force in the FRG 
equation. In fact the moving glass will be dominated by the competition between the random force and the critical 
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force. Although the existence of such a random force has no effect on pinning it will affect strongly the positional 
correlation functions. In particular the relative displacements correlation function 



B(r) = ([u(r) - u(0)]2) (112) 
becomes, in the presence of the random force at large scale: 

^rf{X,y) ~ Aren(Uj / — — ^ — j-^— — - Aren(Uj U { -) (lid) 

J (2i:)'^ (V^qxj+c{qi+qyj c-qv rjvy^ 

where H{0) — est and H{z) ~ ^'^■^"''^/^ at large z. Thus x scales as and the displacements are very anisotropic. 
This behaviour results in an algebraic decay of translational order in d = 3 {B{x, y) log \y\) and exponential decay 
in d = 2. 

It is thus important to estimate Arcn(0)/w which sets the growth of displacements. At large velocity one finds that 
\-cn{0)/v — A{0)/v the bare disorder, which at large v is very small. When decreasing v the displacements due to the 
random force at a given scale increase but eventually saturate and decrease again. This is because in the intermediate 
velocity regime the random force will be reduced strongly, from its original value. 

An estimate, from ^is that: 

A,,„(0) = A(0) - / d/e-'A;(0)^^ ''f J"' (114) 
Jo 27r (4c^e~2'cj^A5 + t;2)2 

and one can simply set e = to get the result in d = 3. 

We have given here mainly the results at large v. These calculations can extended to any v by studying the crossover 
from the Bragg glass to the Moving glass, as well as the intermediate random manifold regimes. This crossover is 



discussed in the Section III 



C. RG study at finite temperature T > 



We now extend the analysis to finite temperature. In principle the FRG equations can also be written for any 
temperature. We will study both the case v = and v > (since no such derivation exist in the litterature) . 

In the V = case the temperature is formally irrelevant. In fact it is dangerously so (see below) as it will cutoff the 
properties of the fixed point (the non analyticity) and thus modify some observables leading to barrier determination. 

Here as we will see temperature for w > is even more so - very - dangerously irrelevant by power counting. The 
dimensional rescaling ( [zsl) yields T T' = Te'^^^'^^^''^' and thus that T is irrelevant. This turns out to be incorrect: 
if one add a small T > onto the T — moving glass fixed point, it will flow upwards very fast (while if T = to 
start with it remains so). This is at least what the FRG seems to indicate. 



1. derivation of the RG equations 



Letting a T > leads to several important modifications in the perturbation theory. The general idea is that the 
disorder will be changed everywhere roughly as A/<-e~ a^'^"-*'^ (with t — > oo). Of course this has to be checked 
carefully which is done in detail in the appendices and we give here only the main results. Since near the upper 
critical dimension (d^ — 4 for v — and du = 3 for u > 0) thermal displacements are bounded: 

r 1 2T 

lim Bo,t - Boo = 2T / ^ = —AdK''-^ {v = 0) (115) 
t^oo cq-" c 

^2T I : [v>Q) (116) 

Jq Cxqi + Cyqfi + ivq^ 

We will denote rescale T by a non universal quantity and define T the "dimensionless temperature" " 

\Bo. = f (117) 
One will need to replace everywhere A/^ by the smoothed disorder: 
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Ak = Akc-^^" K(u) ^ A^e^^^e-^^' (118) 



K 



The divergent part of the correction to the friction coefficient is now (B23): 

5r,^ j tRr-^o.tK^AK (119) 



thus the same as before, except one must use the smoother disorder. ■— ■ 

Let us now compute the renormahzation of the temperature by disorder (for a related calculation see alscO) . As 
was mentionned earlier there is a new and non trivial divergence in the correction to the temperature. Using (54 5^ ) 
one finds: 



^ Jt>o 



) (120) 



When V = this integral can be simplified using the FDT relation 2TRr=o,t — 0{t)d/dtBo^t which gives 6T — 0. 
This yields the RG equation, after rescaling: 

^ = (2-d-2C)T (« = 0) (121) 
For V > Q the second part does not diverge anymore, and the divergence of the first one can be extracted as follows: 

riST^Y. A^(e-^^(^"'-^-)-^^'-l)^ / cHAkK^B^ - Bo,t) (122) 

Jt>0 if "'0 

with: 

2T 



Thus: 



Boo - Bo,t = / _f,-{c.q^+cyqy)^^\t\+^vq.^.t (^23) 

q CxQ^ Cyqy 

6T ^Y^K^Kk I ^r-^^ , . ^ tYk^AkI-^^ , ^ (124) 

Extracting the IR divergent part this yields, after rescaling the following RG equation near d — i: 

^^2-d-2C--^A"(0) {v>Q) (125) 

Sin ce we will also discuss the crossover from u = to w > we also give the following more precise estimate of 
( |120| ) (obtained from the large time behaviour): 

Note that it vanishes, as it should when w ^ 0. This yields the following RG equation for the temperature: 

^ = 2-.-2C , (127) 

™' (47rc,z;)(l + i2i^^^)3/2 

We now look at the corrections to disorder. The calculation is done in the Appendix The divergent contributions 
to the disorder correlator are, adding first and second order in perturbation: 

AP = e-^P" \Ap + e^P" J2 A^e-^^' Ax'e-^^" [K^ f Gir)Gir) + K'K f G{r)G{~r)) (128) 
~P'ApJ2^K'e'^'''" I G{r)G{r)] (129) 



40 



The key point is that using K.K' = (P^ — — K''^)/2 all exponential factors rearrange and at the end everything 
can be written only using the smoothed function Aa'- Using this smoothed function can write: 



5^(u) = A"(m)(A(0) - A(u)) / G{r)G{r) - A'(u)^ / G(r)G(-r) (130) 
This yields the RG equation for the disorder: 

dKiu) 



di 



TA"(u) + {3-d- 2C)A(w) + CuA'iu) + /i(/)A"(w)(A(0) - A{u)) - f2{l)A'{uY (131) 



where /i and /2 are the same coefhcients as in |8^. We have used that the smoothed function A(u) itself has an 
explicit cutoff dependence. Note that this equation is correct for any T and to second order in the disorder. Note 
that it can be obtained also by a small T expansion, assuming T small (and expanding the first order correction to 
A in T). This equation adds to the renormalization equation for the friction coefficient: 

— - 2 - z + A (0)^.-1 (4^^(;)^^A.^^.)3/2 (132) 

with here z = 2. 

2. analysis of FRG equations at T > 

Let us now analyze the FRG equations at T > for the periodic case in an e = 3 — d expansion and in o? = 3. We 
write A instead of A and T instead of T everywhere for convenience. One has: 

= eA{u) + TA"{u) + A"{u){A{0) - A{u)) (133) 

al 

|^^-l + .-A"(0) (134) 

We have absorbed the factor ^J^^ in A{u). Let us first search for a fixed point. We thus assume that dT/dl = 
with T = T* which implies that A"(0) = -1 + e. We search for a fixed point for A(u) - A(0) (as we did for the T = 
fixed point). Let us set A{u) — A(0) — T*g{u) with g{u) > and periodic. One gets: 

(135, 

This is the motion in the potential ^(5). It always has a periodic solution starting from g = 0. Thus the solution is: 

Vig)^±,g-'-^Logil + g) (136) 



This yield a condition since we have fixed the period to be m = 1 

1 _ dg 
2~ Jo ^^Wii) 

The other condition being: 



Vigma.) = (137) 



A"(0) = -l + e (138) 
Both conditions determine T* and A"(0). From this we get the fixed point temperature: 



e2 £2 



^(-p^ax ^^^^^p 81n(l/e) 



(139) 



Thus we find that there is a finite temperature fixed point. This is the moving glass T > fixed point. Though we 
have not investigated in details the stability of this fixed point it is likely to be attractive. Indeed one sees clearly on 



41 



the equation of renormalization of temperature that at high T one expects A"(0) to be small (since A is smoothed 
by temperature), while at low T A"(0) grows very fast thereby T increases. Note that similar finite T fixed points 
were found for other non potential systemS. 
Note that one has also: 

^=em-T*{l^e) (140) 

and thus the random force is still generated (though it grows slower). 
Form this we get the temperature: 

T* = ~ (141) 

l(f^'— , Y 81n(l/e) ^ ' 

Thus we find that there is a finite temperature fixed point. Though we have not investigated in details the stability 
of this fixed point it is likely to be attractive. Indeed one sees clearly on the equation of renormalization of temperature 
that at high T one expects A"(0) to be small (since A is smoothed by temperature), while at low T A"(0) grows very 
fast thereby T increases. 

Note that there will be an interesting crossover at low T where A"(0) will first start to increase violently before it 
finally decreases again to its fixed point value. This will be discussed in following section. 

The case of d = 3 can be studied similarly. As in the case T = one looks again for a solution decaying as 1/Z (as 
we noticed before e and 1 /I play the same role) . 

Here the solution in d = 3 is as follows. One has: 

A(0) - K{u) ^ -^g{u) (142) 

This ansatz gives the same equation as in d = 3 — e with e l/l. Thus g — g*{u,e — l/l) is still the asymptotic 
solution. Also: 

Thus at large scale temperature decays back to T = 0. Indeed the fixed point function is very similar to the T = 
fixed point function except in small layer around integer u. Near the origin the term TA"(0) is of same order as 
eA(0). Thus the main effect of temperature is to round the non analyticity. 



3. Physical results at finite temperature T > 

We now discuss the behaviour of the mobility fiu = 1/7]r and of the I-V characteristics. 

One can compute the mobility from the RG equation by integrating the I dependent solution over all scales: 

In(^) = /' Ar (0)^.-1 .^ ^-fltl\w2 (144) 

Since asymptotically there are only finite corrections to 77 as soon as v > this integral converges. Thus there is 
a non zero asymptotic mobility fiR in the T > moving glass (by contrast with what one would have in the static 
Bragg glass at T > 0). 

However the renormalized mobility fxji will be very small (for experimental purposes) in several cases (i) at low 
temperature (ii) for velocities not very large v < v*. 

A complete calculation of fiR in all regimes can be made by examining the RG equations derived above. Here we 
give an estimate for the low temperature behaviour. A key point is that at low temperature will be determined 
by the short scale contributions. Indeed there must be some continuity with the T — flow, where — A5'(0) diverges 
after a finite length scale, the Larkin length (as discussed above). Thus at low temperature — A"(0) will first 
shoot up near the Larkin length, strongly renormalizing the mobility downwards, before the temperature catches on 
and reduces it back to its fixed point value — A*"(0) ^ 1. Note that this fixed point value corresponds to values 
of disorder much larger than the original disorder. Indeed restoring the factors one sees that (in d = 3) that the 
original disorder dimensionless parameter is A2/(47ruc) ~ ln(a/i?^) <C 1 at weak disorder while asymptotically one 
has A^/(47r'yc) = 1. The global behaviour with length scale is illustrated in the Fig. (p^). 
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l=log(L/a) 



FIG. 23. Behaviour of the second derivative — A"(0) of the disorder correlator as a function of the scale around the T = 
dynamical Larkin length i?^. At T = there is a divergence at which is rounded at T > 0. However — A"(0) still passes 
through very large values before eventually decaying slowly towards its fixed point value. This results in high barriers at low 
temperature as discussed in the text. 

The small T behaviour in d — 3 can be estimated as follows, let us denote Tq the bare value. One has the exact 
equations: 

= TA""(0) - A" (Of (145) 

dl 

= -7A"(0)A""(0) + rA(f^)(0) (146) 

|^ = -1-A"(0) (147) 

We will roughly estimate the scale Rc{Tq) = ae' at which the term rA""(0) starts slowing down the growth of 
A"(0). For this we will drop the term (thus assume that for R < rXTq) the flow is identical to T = 0). 

Integrating the equations then yields (1 — A2l*)^ = Toe~' A4/(A2)^ denoting the bare values of the derivatives of the 
disorder correlator by A2 — — A"(0) and A4 = A""(0). This length scale is very close to the Larkin length and the 
end result is that: 

A9 At Ry 

Thus this argument indicates that as T ^ the renormalized mobility will vanish as: 

m - Moe"^'/"^)" (149) 



D. extensions 



As was mentionned in the Section II B there are many possible generalization of the moving glass equation to a 
larger class of non potential model. Let us indicate here the general FRG equation for these models. We consider 
a dynamical equation with a N component field Ua and static disorder described by some force correlator Aq^(u). 
Then, as shown in the Appendix, the corrections to disorder will be: 

6Aa(i(u) = AaS3-^f,{u){Aai ~ Aa' f3' {u)) / G~^a' (r)G sp' {r) 

J r 

-Aaa.'Au)^l3f}'A'^) / G-ya'(r)G5/3'(-r) (150) 



where we have defined the static response G{r) = dTR{T,r). Note that this formula is valid for a large class 
of models with arbitrary response function. It does not suppose for instance that the random force correlator is the 
second derivative of a random potential. 

^From this a generalized FRG equation can be derived, which depends on the divergences contained in the response 
function. A special case is when v > and isotropy is assumed (the simplest N component generalization of the 
moving glass equation (|^). Then one finds: 
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— = eA^fsiu) + Cu^A^p^^iu) + CA^^,^s{u){A^s{0) - A^s{u)) (151) 
where C is a numerical constant. 

The temperature can be added. In the isotropic case it simply produces an extra term —TAap_-yy{u) in the above 
equation. 

Finally let us close this section by a remark on the nature of the MG fixed point. The role played by Cx (apart 
from a regulator) is a bit mysterious. Indeed one can wonder whether these results could be obtained directly from 
the following static equation for u^y '■ 

vdxU = cld^ + F{x,y,u) (152) 

which contains only the operators which are relevant at the fixed point. Indeed the final results for the fixed point 
of the moving glass involve only (and presumably to all orders in perturbation theory) zero frequency (i.e static) 
propagators (at zer o temperature and at zero applied external force). If one is too naive and uses an Ito scheme in 
X this equation |152| becomes trivial and leads to only a random force. It is probable that when the full non linearity 
of F[x,y,u) is kept, and the x discretization treated properly (which may be subtle) the proper result is found. It 
probably then lies entirely in a proper regularization of the corresponding functional determinant (in a sense the 
dynamics in t that we use here can serve as such a regularization). We leave this for future investigations. 



VII. MOVING GLASS EQUATION IN Z? = 2 AND D = 2 + e 



As stressed in the Introduction, it is important to first study the elastic theory as a function of the dimension d, 
before attempting to include topological defects. Up to now we have studied the Moving Glass equation in an c? = 3 — e 
expansion. This study is of course mostly relevant for the physical dimension d = 3. To study the other physically 
interesting dimension d — 2 another RGLcalcluation can be performed. For the statics v — th(j_RG approach was 



constructed by Cardy and Ostlund CO 



and, with some 



It was later extended to study equilibrium dynamicsL22 
additional assumptions, to study the problem in d = 2 + e. In the case ?; = it yields a marginal glass phase in d = 2 
for T < Tg described, by a line of perturbative fixed points. —Extensions to models with n > 1 components necessary 
to describe a latticaHa and to far from equilibrium dynamicst23 were also studied, ((((((((rsb quelque part )))))))))))) 
In this Section we first show that in d = 2 the CO fixed line is unstable to a finite v on the simplest case of n = 1 
component Moving Glass equation. We derive the RG equations for the case v > 0. We stress that this is a toy 
model since it is clear that in d = 2 additional instabilities to dislocations will occur at the temperatures where we 
can control the behaviour of the model. However it is instructive, as we will see, and provide the first necessary step 
to introduce the other instabilities. 



A. d = 2 



1. RG equations in d = 2 

The equation that we study is: 

{r]dt - Cxdl - Cydl + vdx)urt = /i(r, Urt) + f2{x) (153) 



where fi{r,Urt) is the random nonlinear pinning force with fi{r,Urt)fi(r' ,Ur't') = A{uj.t — Ur'f) and f2{x) is the 
disorder originating from long wavelength disorder (|^) with: 

iMhFq) = + Ao (154) 

In addition such terms are generated in perturbation theory (at least for v = 0) and should thus be added from the 
start. 

We use everywhere the shorthand notation v = rjQV where is the bare value of the friction coefficient. Since we 
are looking at a periodic system one has A{u) = A^e*^". However in d = 2 where temperature is marginal the 
harmonics are relevant at different temperatures. This remains true at u > 0. It is thus enough to consider the lowest 
harmonic 

A{u — u') = g cos(u — u') (155) 
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Perturbation theory is carried in Appendix ^ using the MSR formaUsm. Note that the random forces /2 can be 
ehminated by a shift and do not feedback in the RG (see ^) . In addition due to the tilt symmetry (gaUlean invariance) 
Cx, Cy and v have no corrections. 

One finds to first order in g the following corrections to the friction coefficient, the temperature, and the disorder: 

r+oo 

Srj^g / firTi?(0,T)e~*^("'^) 

f+oo 

S{vT) = g / dT(e-^^("^^) - e-^s(o,r=oo)) 

Sg = ge-^B{0.r=co) (^gg) 

where B{r,t) = {[{urt — uoo]^)o ^-nd R{r,t) — (Suj,t/ShQo)o are respectively the correlation and response functions in 
the theory without disorder. In the case v = the fluctuation-dissipation theorem (FDT) holds: 2TR{t) — dB/dr. 
This ensures that the correction of the friction coefficient and of the thermal noise are exactly proportional and can 
be both absorbed into a single correction of the fric tion coefficient (leaving T unchanged ST = 0). This can be seen 



immediately by integrating by part (156) using (D3). This property does not hold any more when v ^ and T will 
renormalize upward in d = 2. As we will discuss in more details later, the, generation of a temperature due ot disorder 
occuring here is quite different from the notion of shaking temperatur£3 

Similarly than for the statics, disorder will be relevant below a certain temperature Tg. To determine Tg one 
computes the mean square displacements in the absence of disorder 



d'^q 2T 
(27r)2^' 



B{r,t,a) = / 7^^(l-e-=« ^.\t\^^qr+^vq^^.t^^~a q- ^^^^^ 



Bir = 0, t, a) = —iLogi^^] +C + Log[^^-^^^^] - E^[^^-^^^]) (158) 

where the mobility /i = l/j] has been introduced. In (??) we have used the following regularizations: (i) an infrared 
regulator by defining a large time tj^ax but no infrared regulator in momentum q (ii) an ultraviolet cutoff is enforced via 
a gaussian cutoff in momentum. Such a regularization procedure turns out to be extremely conveninient to establish 
the RG equations. Thus according to whether w = or i; > one has the two different large time behaviours: 

4T 

B{0,t,a) = —{lIl[v^^t/a] + C/2) {v > 0) (159) 
2T 

BiO,t,a)^—iHcfit/a'] + C/2) [v = 0) 

J- r 



where — 47rc is the transition temperature of the static system. Remarkably, as can be seen from (159), Tg is half 
of the Cardy Ostlund glass temperature Tc of the statics ! 

Thus the CO line is unstable and both disorder and temperature will be generated. To obtain the RG equations we 
restrict ourselves to the case when the starting cutoff is large enough a^v^ /{4:C^) ^ 1 (or the velocity large enough) 
so that one is already in the asymptotic regime. Of course at small velocity there will be a complicated crossover 
where the short distance properties are dominated by the static solution, but the large distance properti es w ill be 



again given by the present RG equations. Introducing the dimensionless coupling constant g = ga/{vTc) ( |156[ ) and 
( |159D allow to obtain the correction to rfT as: 

which yield the RG equations upon a change of cutoff a' — ae'": 

''^'^ (l-^).g(0 + O(,^) (161) 



v{l)Tc dl 



dl ' T, 
~9il> 



1 dmni)) _ ,,„^_c/2 



the last legality being valid near the transition at T — Tc/2. (161) can be compared to the Cardy-Ostlund RG 
equationiy for v — 0: 
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f =2(1-^)5 + 0(3^) [v^O] (162) 



Finally, because of the exponential decay of the response function 



R{0,t) = , , ^ e^^<»^+<=^-) (163) 



(as in Section ??) which cuts all divergences, the correction to rj in (156) is only a finite number. There is thus no 
renormalization of "q at large scales. The final RG equations in d = 2 read: 

- = (l--).g + 0(.9) 
^ = 2013 + 0(5^) (164) 



with Ci = e '^/^ is a nonuniversal constant. 



B. Analysis of RG equations in d = 2 

Let us analyze the RG flows. We introduce the reduced temperature r — (2T — Tc)/Tc and g = 2Cig. The RG 
equations are: 

The trajectories are the arches of parabolaes represented in Figure ^ centered around Tc/2, of equation: 

5-50-^(^0 -t') (167) 

(note that close to T = Tc/2 these trajectories are not modified by the higher order terms). 

As can be seen from Figure ^ if one starts at small disorder with temperature Tc/2 — AT, both disorder and 
temperature first increase pushing the system in a region where the disorder is irrelevant, ending up with a disorder 
free system at about Tc/2 + AT. This has several physical consequences: 

(i) at finite velocity the effect of disorder is weaker than in the statics, which manifest itself in the RG equation 
since weak disorder becomes irrelevant for T > Tc/2 a region which is already deep in the glass phase in the statics. 
This effect is analogous to the dimensional shift from due = 4 in the statics to due — 3. 

(ii) However we still find a transition at T = Tc/2 below which disorder is relevant and grows under RG. That such 
a region where disorder is relevant exist in c? = 2 is compatible with the FRG findings in d = 3 — e and clearly shows 
that even in motion one still has to consider the effect of the random potential. However due to the importance of 
the thermal effects in d = 2, at large enough scales the disorder will stop being relevant since the temperature also 
increases. The lenght scale ^ at which disorder becomes again negligible can be estimated from the RG and reads: 



^^(^)l/(4ro) (168) 

9 

^ cna become extremely large when the disorder is weak or when one starts at low enough temperatures. 

(iii) Finally, we find that disorder generates an additional temperature. This renormalization of temperature is 
physically very different that the "shaking' temperature" of In particular the value of the generated temperature 
in our case does not depend on the strength of the disorder but on the temperature itself and the distance to Tc. In 
particular, and in similar way than for the FRG, if one had started at T = no temperature is generated as can be 



seen from (172) 



Using the RG flow one can compute the displacements. For the connected correlations one finds: 
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{[u{x) ~ u(0)]2) - {u{x) - m(0))2 - T*{T,g) \nx 



(169) 



We have used the exponent at the fixed point, which is a correct procedure because the fixed point is approached 
fast enough as tq — t(0 - roe-"'o'. These correlations are nonmonotonic as a function of T with an almost cusp 
(rounded by g) at 7^/2, and increase below Tc/2. 

Given the present result near Tc/2 possible topologies of the fiow in d = 2 can be proposed as shown in Fig. In 
any case the relevance of disorder at low temperature is a very good confirmation of the presence of the moving glass 
phase at least at T = for which temperature renormalization effects are absent. 

However in d = 2 thermal RG effect are obviously important and the question of whether the MG phase exist at 
finite T arises. Two fiow diagrams are in principle possible. In one of them a low temperature phase exist (MG). In 
that case an additional fixed point which controls the transition is necessary. In the absence of such a fixed point the 
moving glass will always be unstable due to temperature renormalization. This last scenario is supported by the FRG 
results in d = 3 — e and by the next Section. 



g 



? 




FIG. 24. RG flow diagram in d — 2. The flow is circular around a new instability temperature at T = Tg = Tcol'i.. The 
Cardy Ostlund line of fixed point of the statics which start at T = Too in unstable when i; > 0. For d = 2 + e there is a new 
finite T moving glass fixed point (presumably attractive) at g ~ e (the resulting spiraling flow is not shown). Continuity with 
the FRG result suggests that this fixed point moves upward to the T = axis as d goes from d = 2 \,q d = Ind = 2a 
zero temperature moving glass fixed point is expected at infinite g — — A"(0) (if the lower critical dimension die for the T = 
moving glass is die < 2). 



C. Moving glass equation in d = 2 + e 

We now follow Goldschmidt et al. and try to continue the above RG equations to d = 2 + e. The e simply shifts the 
dimensions of the operators. We stress that this is based on the assumption that the RG functions are well behaved 
around d = 2. It has been used though in other cases such as the 0{n) model. 

Taking into account the dimensions, one readily obtains using the same reduced variables, to lowest order: 



dg 

U 
dr 

m 

dr] 



= -{^ + r)g-bt 

^-2e + g (170) 
= 



with b = B / [ACD being a universal (regularization independent) number (by analogy with CO). 
These equations now have a new fixed point at 5 = 2e and r = — (1 + 2b)e (note that B /Ci is universal). To lowest 
order the eigenvalues are: 

\±=-be±iV2k (171) 
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Thus without needing to compute the coefficient b we know that there is a fixed point, and we know the leading 
behaviour of the eigenvalues X± ~ iiy/e. Such spiralling fixed points have been obtained in other problems (e.g 
Ref.t2il). However to know whether the fixed point is attractive or repulsive one needs to know the real part, which 
is controlled by the 0{g^ , gr) terms in the RG equation. For instance, the RG equation will contain at least: 

J = (-2e + 5 + cg2)(l + r) (172) 

Inspection shows that b actually controls the leading behaviour of the real part. So this is the only non linear term 
we need to compute. Results from FRG (see') and static CO lead to expect that 6 > 0, but to settle the question 
and obtain the value of b an actual calculation along the lines oM is needed. 

It is tempting to associate this fixed point to a finite temperature moving glass fixed point (analogous to the new 
finite T fixed points found in recent manifold studieS). 

Finally we not that at this fixed point one will have z — 2 since we find that drj/dl = 0. (thus there will be large 
but finite barriers). 



VIII. TOWARDS A COMPLETE DESCRIPTION OF ELASTIC FLOWS 



In this Section we investigate Models II and Model III. 



A. Study of the complete dynamical equation in elastic limit (Model II) 

In this Section we come back to the Model II which contains both degrees of freedom transverse Uy and along 
motion Ux and the full time dependence of the pinning force. If one describes the elastic flow of a solid, i.e in a regime, 
or a range of scales where there are no topological defects, this model improves on Model III (see Section Intro for a 
discussion of the regimes where it will be useful). As discussed in Section |^ it is still an approximation, but a rather 
good one, of the full (but intractable) model of the elastic flow (Model I). 



Since Model II is still quite difficult, our aim in this Section is more limited than in Section |VI[ We show that the 
main features of the transverse physics of the Moving Glass embedded in Model III are also present when the degrees 
of freedom along the motion are added. In fact we show explicitly that in Model II the RG equations for along y 
remain identical to the one of Model III ! The two important issues we discuss are the one of the existence or not of an 
extra temperature generated by motion, and of the generation of a static "random force" . We perform perturbation 
theory up to second order in the disorder and examine the new terms generated as well as the divergences. We develop 
an approach which allows to treat all harmonics A'^ of the disorder correlator. First, we find that a static "random 
force" is generated in the direction of motion. This may seem surprising at first, because first order perturbation 
theory gives a pinning force along x which rapidly oscillates. However, as our calculation shows, to second order the 
various washboard frequency harmonics interfere to produce a static random force. Second, we identify the divergences 
in perturbation theory and follow the evolution of the full disorder correlator under renormalization. We show that 
up to some details the resulting picture is close - if not identical - to the one given by Model HI. We work at T = 
but the approach can be extended to T > along the lines of Section [v| . 



1. general properties 
The equation of motion ( |4l| , ^ ) keeping the time dependent terms reads: 

(i?"^)"f.,i,uf,j, = fa -rjVa + Fa{r,t,Urt) (173) 

More specifically we will be interested in the model for a triangular lattice {n ^ 2 component model) with the force 
applied along a symmetry direction. The equation of motion reads for lines in d = 3: 

VxxdtUx + m'vdxUx + (ceeV^ + Ciidl)ux + (cn - CfiQ)dx{dxUx + dyUy) = / - rj^xV + FP^'^{r,t,u) (174) 
rjyydtUy + TjovdxUy + (ceeV^ + cnd'^)uy + (cn - cee)dyidxUx + dyUy) = FP''"{r, t, u) (175) 

(176) 
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and setting C44 = for points in d = 2. Note that rjxy = from the symmetry {uy —Uy, y —y)- We have 
aUowed for different rjxx and rjyy since, even if they start identical (and equal to 770) they will not remain so under 
renormalization. The statistical tilt symmetry ensures that the elastic coefficients and rj^v rem ains uncorrected. Note 
that in later calculations it will be convenient to rescale the z direction, setting z — ca/ cqqz' . 

The correlator of the pinning force can be written as: 

F„(r,t, u)Fp{r',t',u') = S'^ir - r')A"^(w ~ u' + v{t - t')) = 5'^{r - r') ^ ^^js ^-^K.{u-u■+v{t-t■)) ^^^^^ 

K 

It contains all lattice harmonics K . Due to the modes with Kx 7^ it is an explicit periodic function of time 
with frequencies all integer multiples of the washboard frequency 27ru;o = v/a. In addition it contains non linear 



static components Kx = 0, Ky ^ (which lead to the Model III treated in Section VI). Finally it contains a static 
u-independent component, Aif=o which we call - by definition - the static "random force" . We believe this is the 
"correct" definition of the random force. Note however that one must distinguish it from A(m = 0). If one was to 
compute the displacements correlator to order one would need A(m = 0) to second order, which is a different 
calculation. The important point is that the static "random force" is strictly zero in the bare model Ax=o(0) = 0, 
but that it is generated in perturbation theory, as we show below. 

The idea behind the method presented here is that strictly at T = all the time dependence remains strictly 
periodic to all orders in perturbation theory in the disorder (plus a static part). Only frequencies multiple of ojq 
can be generated. Indeed to lowest order FP^"'{x,y,t,u = 0) is periodic in time, and yields a u periodic. Iterating 
perturbation theory thus leads only to periodic u and F^^"'{x, y, t, u) with integer multiples of loq. This property allows 
us to construct a closed RG scheme of the above model with a renormalized disorder which remains of the form 

An immediate consequence is that no temperature is generated when T = at the start. Temperature is defined 
as the zero frequency limit of the incoherent noise (or using the MSR formalism the vertex function r/T = Tuu{q — 
0,uj ^ 0+)). Thus here one has: 

rjST = fdrY, {Ken)f e-'""-^^ = (178) 

The static random force K = Q mode leads to a 5{lj) part in u and thus is distinct from the temperature. The 
fact that no temperature is generated strictly at T = is a rather strong property of the elastic flow. In physical 
terms in the T — Q elastic laminar periodic flow all particles strictly replace each others after time tq — a/v, i.e 
^i^,iyi^ + r) = Ri^+i^i^{t + t) where i = ix,iy are integer labels for the particles. It is clearly possible that this 
laminar periodic flow becomes unstable to chaotic motion. It is still very interesting to investigate this laminar 
periodic flow. We will thus proceed here assuming here that instabilities to chaos happen only at finite large enough 
disorder, or at large enough scale. We reserve the study of the stability of this flow (chaos, non perturbative effects, 
etc..) to future study. Finally note that this periodic flow is allowed by the assumed absence of topological defects 
in the system. Dislocations, if present would probably ruin periodicity and introduce a small additional temperature 
(though this remains to be investigated). Finally, it is. interesting to note that the present considerations are in stark 
contradiction with the Koshelev-Vinokur argumentsEZi about the generation of a "shaking temperature" . 

We start by establishing the possible form for the disorder correlator (at any order in perturbation theory) based 
on the symmetries of the problem (Model II). This is highly necessary here because the bare disorder is potential 
A^'^(O) = gKKaKjj, but it does not remain of this form in perturbation theory ! We are interested here in the case 
when the velocity is along a the lattice direction. Our analysis here is very general, and we will specify when we 
apply it to the case of a triangular lattice. More details are contained in the Appendix using the more rigorous MSR 
formalism. 



The symmetries are as follows. First one can exchange t with t' and u with u' in (VIII A 1) and relabel the disorder 
term. This gives A"'^(w) = A^{v). This is by construction and from the speciflc dependence in t and t' of the 
disorder. Second, the action must be real and thus A"'^(?;) — A'^{v)*. Third, the symmetry Ty { Uy ^ ~'^yj 
y - ^y - -Uy) yields that A^£ _^J^) = A^'^ ,,J^), A- __^J«) = A- ^^^(^;), A^^^k^^) - ^^Z^kS^). 

= Similarly, because of Tx one finds A^^^^^^^(-w) = A|^^^^^^(w), A'L^^j^Ji-v) = 

^k^^kM. ^"i\^,KS~y) = ~^Z.kM^ = -^k^^kM- Note that the global symmetry T.T, 

implies that A^^'^(w) = A^'^(-w). Thus one finds that one can split the disorder correlator into: 

Al^{v) = A%%{v) + Af^{v) (179) 

where A^'^(w) is real, symmetric in a/3, even in K, and even in v and /S.°^j^{v) is imaginary, antisymmetric in a/3, 
odd in K and odd in v. This naturally leads to the decomposition 
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A^^^(i;) = A^v^S^p + A^K^Kp + Afv^vp (180) 
^A^if(«) = ^^fivo^Kp ~ vpK^) (181) 

where all Af are even in K and v and real. The bare disorder has only non zero and thus posses the extra 
symmetry A"'^ = A'^ or equivalently {u —u, v — > ^v). Because of the convection term vdxU in the equation of 
motion, this additional symmetry will not hold to higher orders in perturbation theory. It is natural to suppose that 
to any fixed order in perturbation theory the Af are regular when w — > 0. Thus in the limit v = one recovers a 
strictly potential problem (all terms except Af vanish). 

Finally note that Model 111 is a particular case of Model II which corresponds to the foUwoing choice of bare 
parameters: (i) isotropic response cn — cge (ii) A^ = 0. Then clearly the equations along x and y decouple. 

Another particular case that we note "en passant" is to start from A'^ — gxKaKfj and elastic elasticity. Then the 
equations are only coupled through the time dependent part of the non linear pinning force along y which depends 
on Ux (one has A|^ 7^ 0). It would be interesting to check whether this is enough to change the behaviour. 

At r = the lowest order corrections to the disorder come from second order perturbation theory. The calculation 
of the effective action to second order for the most general model is performed in the Appendix We first give 
the full, awsome looking expression for it (|l82| ), and then will extract some of its main features. A full consistent 
treatment is left for future work. The result of Appendix |^ is that the correction to the disorder is: 

5Af = ^K^KsAf ^ / AP^,R^P{-ivK\ -q)R^\ivK', q) (182) 

A" "'9 

+ \Y.I^^- ^')^(^ - K')sAf^j,,[AP^,WP{-ivK', ~q)R'\ivK', q) + Ap\,R'^p{ivK' , ~q)R'\~ivK' , q)] (183) 

^ K' ''1 

-Y^l K's{K' - KyAZ^K'-RR'^i^^K', q)R'\iv{K' - K), q) (184) 

K' •'1 

+KsA%pY. I K'^^K-R'^iHK + K'),q)R'\ivK' ,q) - K^Af!_\Y. j K'.A^PR^PitvK' , q)R'\iviK' - K), q) (185) 



2,. generation of the static random force 
Setting K = in the above formula one gets: 

dAf = ^ / K.KsiliA'^Jl, + Al^)AP^R'^P{^ivK, -q)R'\tvK, q) - A^''A^^i?^''(zT;i^, q)R'\ivK, q)] (186) 



K 



This is the general expression for the static random force correlator, i.e a u independent gaussian random term 
F(r) in the original equation of motion. 

The first remark is that the random force cannot have crossed correlations, i.e 8^^^(v) = 0. Thus are two 
independent random forces one along x of strength Ag^ and one along y of strength Aq^. This is consistent with 



formula (18C) which gives l^^(v) = iA^ KyV and vanishes for K — 0. It can also be seen explicitly on the above 
expression which is found to be symmetric in a/? (using all the above symmetries). Since we know that Aq^{v) must 
be antisymmetric in xy (from the above Ky —Ky transformation) it must vanish. 

Since the above expression is still complicated as it does involve all disorder harmonics, one must carefully distin- 
guish between: 

(i) the static random force generated to lowest order in the bare disorder O(A^) (i.e at the initial stage of the RG). 
To this order one can use the hare disorder in ( |186D and the resulting perturbative expression of the random force 
(evaluated below) is found to be well behaved and without IR divergences. Of course once even a small finite random 
force is generated it is relevant by power counting and must be taken into account (though it does not feedback in 
the RG for the non linear disorder). 

(ii) the static random force generated to higher orders in perturbation theory (or at the next stages of the RG). 
There we find IR divergences. This means that non trivial corrections originating from the non linear disorder will 
have to be also taken into account to estimate the random force generated. This is done in the next subsection. 

Thus we start by giving the expression of the random force generated from the bare disorder (i). 
Setting A°^ = QKKaKp in (186) one finds: 
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^J2f KaK0gj^{{K.R{~ivK,~q).K){K.R{ivK,q).K) - {K.R{ivK,q).Kf ) (187) 

Note that it does vanish for u = as it should since the problem becomes potential in that limit. 
One can first specify this formula for the case cn = cge- We also set rjxx = Vxy (which is consistent since we are 
just looking at the lowest order contribution in perturbation theory). This yields: 



= V / X^j^.i^.g^ ,, + ,,,, (188) 



One thus finds that a random force is indeed generated, i,e there is a positive Ag^ and Aq'^. One checks that indeed 
Aq^ = 0. The integral for Aq^ is infrared divergent, as discussed above (which is natural from the analysis in Section 
0) and will be examined in the next Section. 

Let us estimate the magnitude of the static random force generated along x. From the above expression ( |188| ) one 
finds in the case = relevant for point lattices (in d = 1, 2, 3): 

AS-»^ E A-,|- + ^L_ J2 (188, 

' K.\K^\>K,„ ' K.\K^\<K^ 

where Km = \Ko\ max(l, Vcr/v). There is a crossover velocity rjoVcr ^ cir/a. One has defined Cd — ^Sd/ ((i(27r)'^) and 
= Sd{4: — (i)7r/(8(27r)'' sin(7r(i/4)) and n — 2 for triangular lattices. This yields for ry ^ a: 

« min(l, i^f/^) (190) 

(»7ow)V/ " Vcrrf 

In the case cge *C cn one can simply retain only the transverse mode (thus setting c = cge in the above formulae. 
In the case dz — 1 (relevant for line lattices) we will only give the large v estimate (valid for v ^ Vcr = ceen/a). It 
reads: 

Let us stress again that we have defined the random force as A^'Lq (which we believe is the "correct" definition). 
If one was to compute the displacements correlator uu to order A^, in order to show that the displacements "feel" 
the random force (e.g and grow unboundedly in the a; or y direction) one would need A{u = 0) to second order (and 
to be consistent the response function corrected to order A) which is a different calculation. While it will make little 
difference in an order of magnitude estimates for the effect of the x random force, it is drastically different along y 



(one quantity IR diverges and the other does not - see Section VI, A(ii = 0) has no divergence) !! 



3. RG study of Model 11 

We will now look for the divergences in perturbation theory which appear in Model II. We will address only the 
case V > 0. Let us look again at (|182|). It contains infrared divergences of the same type that was discussed in Section 



VI. These divergences occur only for momentum integrals which are both (i) zero frequency integrals {K.v — 



terms) (ii) involve q and —q. These are of the form: 

D^p.M = / G^P{^q)G'\q) (192) 
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where G°''^{q) = R"'f^{uj = 0,g) is the static response function. Here it has the form: 

G"'^(9) = E-7%^ (193) 
^ ci[q) + ivqx 

where I — T,L index the transverse and longitudinal projectors and elastic eigenenergies as defined in (). The key 
point (for n — 2, e.g triangular lattices) is that one has here qx ^ Qy ^ ql- Thus all projector elements are subdominant 
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for small q (i.e kill the IR dive rgences) except the two elements Pj^ ~ + qI) ^- Thus the only IR 

divergent elements among (|l92| ) are Dyyyy, D^xyy, Dyyxx-, D^xxx- Explicit calculation of the divergent parts gives: 

[I f 1 

Jqy,q. 2^(0119^ + CAAqi) Jq^^q^ 2v(C66qy + caqi) 

Dxxyy = Dyyxx - .((cn + Cee)?^ + 2c44g.2) 

Let us now analyze the consequences of these divergences. ^From ( p_82[ ) the relevant corrections to disorder are: 

J A^'^ = i?^p,5A {~K^KsAfA''p^ + {K-PUK-P)sAf_p^[A''p^ + A'_^p] (196) 

P=(0,P„) 

Note that there is no component Ap^(-Q p -j in the bare action but that it is generated in perturbation theory. This 
expression simplifies because of the symmetries discussed above and the fact that only the terms ( |194| ) appear in {19t ) 
(the term Ap^ always occurs in sums symmetrized over pX and one can use that A^ p =o ^ —A^p p to cancel all 
crossed Dxxyy tcrms). From what remains one finally obtains the following RG equations: 

-^ = eAf+ J2 {AnAypl{-KyKyAf + {K-P)y{K-P)yA%%)+AeeA^/KxKx{Af_p-Af) (197) 



with All = l/(47ru^ciiC44) and Aqq — 1/ [ATiv.JcQ^tm) if one uses the same regularisation as in Section VI It turns 
out that these complicated looking equations have some simplifying features. 

The RG equations ( |199| ) are thus the generalization to Model II of the RG equations of Model III and contain both 
the physics of Ux and Uy. They show that in the Moving Bragg glass, non trivial non Ikiear effects also occur in the 
direction of motion. A more det ailed study of these equation will be given elsewhereEHI. Here we give their salient 
features. The RG equations (199) exhibit remarkable features. First, the subset of these equations for Aj'p ^ ^ closes 
onto itself !! Indeed one finds: 

dA^^ 

Py 

This is exactly the RG equation of the Moving glass Model III !! Thus we have shown that the Model III describes 
correctly the transverse physics, as announced, within Model II. 

One can also write the above RG equations as coupled differential equations for three periodic functions of two 
variables A"^ {ux,Uy) with a,/? = xx,yy,xy. We will temporarily use the shorthand notation u = Uy and v — Ux- 
Let us denote A^ =o(") ^i(")- Let us absorb the factor e/An in the A. These coupled RG equations become: 

'^^"^['''''^ = A^0{u,v) + (Ai(0) - Ai{u))dlA^p{u,v) + ^{A2{Q) ~ A2{u))dlA^p{u,v) (199) 

with 7 = Aqq/Ah and Ai(u) = J dvAyy(u,v) and A2{u) = J duA2.j.(u, w) 
Or equivalently: 

^^1^ = AkM + ((Ai(0) - Ai{u))dl - (A2(0) - A2{u))^KI)AkM (200) 

(dropping temporarily the a(3 index). 

The function Ai(u) obeys the closed RG equation: 

= Ai{u) + A'/(u)(Ai(0) - Ai(^.)) (201) 

di 

and thus Ai(u) converges towards the Moving glass fixed point A*i[u) — AJ(0) = \u{u — 1). Let us now examine 
the behaviour of the other components of the disorder with Kx ~ (and a/3 — xx, xy). They obey the equation: 

= A^p{u) + A%{u){Am - Ai{u)) (202) 
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One easily sees that Aap(u) also becomes non analytic beyond the dynamical Larkin length since: 

= A::,^(0)(1-A'/(0)) (203) 



dl 

and thus the divergence of A"(0) — oo at the Larkin length 1^ = Ini?^ (see Section ^ ) implies that A^^(O) also 
diverges. Thus Aq^(m) does become non analytic. Note that the solution of this equation is simply A2(/) = (^)Ai(Z). 
It is then easy to show that A"'^(m) = CA^(it) is a stable fixed point solution (up to the usual growing constant). 



Indeed, inserting the fixed point value Ai(m) in ( 203 ) one finds exactly the stability operator of the original fixed 
point. This was discussed in Appendix ^ and the results there shows the stability with in the space of (non analytic) 



periodic functions. The constant C can be determined. Indeed one can also use (203) to study g — A'^^(0+). One 
gets g ^ /;^,^ dl{l - A'/(0+,0)- The exponential convergence of Ai towards its fixed point implies that 5 is a finite 
constant. This determines C. Thus at the fixed point one can replace A2(u) = CAi(u) in ^00| . 

Thus we have shown that the RG equations in a Moving Bragg Glass decouple compl etely along y, giving back the 



one of the generic moving glass studied in Section VI . A complete study of the system [199| will be given elsewhere. 



B. Full model for the elastic flow (model I) 

We now come back to the problem of establishing the correct long wavelength hydrodynamic description of a moving 
structure with some internal order described by a displacement field Ua{r,t). The first step is to write an equation of 
motion which contains all terms which are (i) allowed by symmetry and (ii) a priori relevant in the long wavelength 
limit by power counting. We carry this step here, check that all these terms are indeed generated in perturbation 



theory from the original equation of motion 27, and estimate their magnitude. The second step, which is to solve the 
universal large distance physics of such an equation turns out here to be a formidable task, which goes beyond this 
paper. 



As we have discussed in Section IV the problem of driven lattices posses some additional "almost exact" symmetries 



which allow to si mplify th e hydrodynamic description and to extract some of the physics. This has lead us to study 



Model II (Section [VIIIAD which posses the statistical tilt symmetry forbidding many terms, and the simpler Model 



III (Section VI) (the Moving glass equation) which we believe contains most of the physics of moving structures, i.e 



the physics of the transverse degrees of freedom. 



As discussed in Section IV the only exact symmetries of the problem, for motion along a symmetry direction of the 
moving structure {x axis) are the spatial inversions along the directions transverse to the velocity. Power counting 
shows that the general form for the equation of motion in d < 3 is (Model I): 

■qapdtup + LlpdjUp - Clpd^dsup - K^^p^dsupdeU^ = F^^^r, u, t) + Ca{r, t) + fa - iiafivp + 5fa (204) 

where the velocity v is fixed by the convention that d/dt j\,ujf = and / is the applied force. The KPZ terms 
are allowed because Ua is dimensionless at the upper critical dimension duc=3 in a power counting at T = 0. This 
can also be seen by writing a MSG formulation of ( ^04[ ). The above equation of motion is not fully complete 
unless one specifies the relevant disorder and thermal noise correlators. The thermal noise has a gaussian correlator 
rjarjfs = 2{'qT)afj5{t — t')5'^{r — r'), in general anisotropic. The correlator of the pinning force F^^^^{r,u) (which has 
zero average) is gaussian and of the form: 



F^i'^ir, u, t)FI^'^{r', u' , t') = A^piu - u' + v{t - t'))S''{r - r') (205) 

where Aa[3{u) is a periodic function in the case of periodic structure. Note however that in general, Aq^(u) is not a 
potential disorder. 

Compared to Model II, cutoff effects which break the exact statistical tilt symmetry allow new terms to be generated, 
such as linear terms which correct the original convection term vd^u and non linear KPZ type terms. The linear 
terms are obviously relevant and the non linear KPZ terms, in presence of the disorder, are relevant for d < 3. Thus 
once they are generated, even if their bare value is very small it may grow under RG and become important at large 



scale (a full solution of the RG equations for (204) would be needed in order to conclude). However one may guess 
that since the statistical tilt symmetry is "almost exact" , the scale at which these new term are able to change the 
physics compared to Model II and Model III (if they do) may be very large here. Finally note that there are also 
small corrections to the elastic matrix. 

The above approach consists in writing a model independent equation ( ^04] ) based on symmetry arguments. It may 
be useful in proving the universality of the behaviours of various structures. However, in many cases it is much more 
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instructive to start from a given simple model without disorder, such as ^ and to estimate the bare values of the 
new terms to first order in a perturbation theory in disorder. Indeed, in the absence of disorder the above equation 
of motion reduces to: 

€p{dtup + vd^up) = {CyJijd^dsup + U- itpvp + Co(r, t) (206) 



disorder to all terms of the equation (204) 



We have thus computed in the Appendix B 3 the corrections to first order in perturbation theory with respect to 



Though the above equation (210) looks formidable many terms are zero from the exact inversion symmetry. We 
thus now explicity specify the terms allowed in the equation of motion, for the case of an elastic structure described 
by a n = 2 component displacement field [ux, Uy). In d = 2 and d = 3 the equation along y should be odd under 
the inversion {uy — s- —Uy, y —y) and also under (z —z) while the equation along x must be even under these 
transformations. This yields in d = 'i: 

rjyydtUy + vidxUy + V2dyUx = {cidl + C2dy + csdfjuy + ad^dyUx 

{aidxUx + a2dyUy)dxUy + [a^d^Ux + aidyUy)dyUx + a^dzUxdzUy + Fy'-^r, u, t) + Cy(r, t) 

rixxdiUx + vsdxUx + VidyUy = {cr^dl + c^dy + C'jd'^)ux + csdxdyUy + aeidxUx)^ + aridyUx)^ + a^idxUyY + ag{dyUy)'^ 
+aiodxUxdyUy + aiidxUydyUx + ai2{dzUxf + aiz{dzUyY + +FP''"'{r,u,t) + fx - VxxV + dfx + Cx{r,t) (207) 

and the same in d = 2 with C3 — ~ ci — ai2 — 013 — 0. 



While the full analysis of (21C) goes beyond the present we now present some arguments showing that the linear 
terms, by themselves are unlikely to alter drastically the physics of the moving glass. The physical interpretation of 
the linear terms is that now the local velocity explicitly depend on the local strain rates of the structure. Though 
this may generate instabilities (see below), unless an instability occurs, this is unlikely to alter the transverse physics. 
Indeed one can see that small additional linear terms do not remove the divergence in perturbation theory which was 
the hallmark of the moving glass. Let us write vi = v + w, V2 ^ aw, v:^ — v — w, V4 = bw and consider w as small 
compared to v. Also we choose for simplicity isotropic elasticity ci = C2 = C4 = C5 = C7 — C3 = c, cg = C4 0. Then 
the eigenvalues of the a; = (static) response matrix are: 

(q) = i{qxv ± w^ql + abq^) + cq^ (208) 



and eigenvectors {6ux,Suy) — (bqy,qx ± yQx + abqy. Note that one must have ab > otherwise an instability 

develops. The perturbation theory result shows that indeed a6 > at least to lowest order in the disorder. One finds 
for instance that the integral which is the key of the FRG equation for the transverse modes () reads: 

GMGyyi-q) ^ I ^— (209) 



9 (c2g4 + {yq^ -f w.^ql + abqlY){c'^q'^ + {vqx - w^g^ + abqlY) 



As is easily seen this integral is infrared divergent for d < 3, logarithmically in d = 3 and power like in d = 2 

al 



(since w is a small correction to v). The divergences occur in two hyperplanes qx = ± . ^ qy which are tilted 



symmetrically with respect of the direction of motion. 

Let us now reexamine the Moving Glass equation, i.e Model III, and ask whether cutoff effects (absence of exact 
statistical tilt symmetry) will generate new relevant terms. By definition this equation involves only Uy and thus the 
only a priori relevant terms allowed by symmetry are: 

riyydtUy + vidxUy = {cidl + C2dy + CAdl)uy + a2dyUydxUy + i^^'" + Cy (210) 



Now, we have shown in Section VI that at the Moving glass fixed point dx ^ dy and thus the KPZ term 02 is 
irrelevant by power counting. Note that a cubic KPZ terms {dyUyY dxUy is allowed by symmetry but again irrelevant 
near d = 3p-Thus the Moving Glass equation Model III is stable to cutoff effects and perfectly consistetit. This lead us 



to claim ioHjihat while previous descriptions of moving systems, such as manifolds driven in periodicii23 or disordered 
potentialsOO, focused on the generation of dissipative KPZ, such terms are much less important in the Moving glass 
equation, a problem which, because of periodicity, belongs to a new universality class. 

Finally one can also reexamine Model II, the physics of which is presumably very similar to Model III at least as 
far as the transverse degrees of freedom are concerned. This will certainly hold below a (large) length scale max(L;i„, 
Lrpz)- Above one must worry about the new terms. Power counting at d = 3 (where Uy and Ux are dimensionless) 
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in the equation for the transverse degrees of freedom Uy (using that dx ~ dy in the absence of the new terms) shows 
that the only KPZ term marginally relevant at the Model II fixed point in d = 3 (and thus the dangerous one) is the 
term a^dyUydyUx- Note also that the linear term v^ dyUr also becomes relevant there and will change the counting. 
In the end it is probable that all terms in equation (|21C| ) will have to be treated simultaneously to get the physics 
beyond max(Lii„, Lkpz)- 

Finally we note that the arguments given in the previous Section about the fact that no temperature is generated 
at T = are unspoiled by the new terms generated here in the equation of motion compared to Model II. That these 
new terms may lead to other instabilities of the periodic time ordered flow resulting in chaotic motion is clearly an 
interesting possibility deserving further investigations. 



IX. CONCLUSION 



In this paper we have studied the problem of moving structures (such as vortex lattices) in a disordered medium 
following the new physical approach developped in Ref. 68 , The main new emphasis in that approach is that because of 
degrees of freedom transverse to motion, periodic structures have a radically different physics than more conventional 
driven manifolds. The main consequence of our study is that the moving structures remain different from a perfect 
structures (e.g. a perfect crystal) at all velocities (for d < 3 for uncorrelated disorder). In particular they still 
exhibit glassy behaviour. The moving configurations can be generally described in terms of static channels which 
are the easiest paths in which particles follow each other in their motion. We have introduced here several degrees 
of approximation of this problem, embodied in several models. The simplest one. Model III, introduced in Ref. ^ 
focuses only on the transverse degrees of freedom. A more complex one Model II also contains degrees of freedom 
along the direction of motion. We have studied these models using several renormalization group techniques as well 
as physical arguments. All our calculations and results confirm that focusing on the transverse degrees of freedom 
(Model III) gives the main physics for this problem. Indeed we have shown that the more complete Model II leads to 
the same physics as Model III. 

At zero temperature we have explicitely demonstrated that the physics of the moving glass is governed by a new non 
trivial attractive disordered fixed point. Using the RG, we have explicitly demonstrated the existence of the transverse 
critical force predicted in Ref. which is related to the nonanalytic behavior of the renormalized disorder correlator 
at the fixed point. Its actual value, computed from the RG coincides with the estimate given in Ref. phased on the 
existence of a dynamical Larkin length R^. We have also found that at T = no temperature is generated because 
perfect time periodicity is maintained. 

We have also shown that a static random force is generated both along and perpendicular to the directions of motion. 
As a consequence relative displacements in both x and y directions grow logarithmically in d = 3, but algebraically 
in c? = 2. Thus in ci = 3 at weak disorder or at large velocity, the moving glass retains quasi-long range order and 
divergent Bragg peaks. Since the decay of translational order is very slow in c? = 3 we predict that a glassy moving 
structure with quasi long range order and perfect topological order in all directions exists: the Moving Bragg Glass. 
The determination of its physical properties is the main result of this paper. This phase is the natural continuation 
to non zero velocities of the static Bragg Glass. 

We have investigated the effect of a non zero initial temperature. We found that the moving Bragg Glass survives 
at finite temperature as a phase distinct from a perfect crystal and with properties continuously related to the zero 
temperature moving Bragg Glass. At low temperature the moving Bragg Glass still exhibits a highly non-linear 
transverse velocity-transfcrse force characteristics with an "effective transverse critical current" (in the same sense as 
for the longitudinal critical current). Note however that at T > the FRG calculation indicates that the asymptotic 
behaviour is linear but with a strongly suppressed transverse mobility at low temperature. 

In addition the existence of elastic channels provides a new and precise way to look at the problem of generation 
of dislocations in moving structures. The natural transition is now a decoupling of the channels with dislocations 
decoupling the adjacent layers. It is indeed easier to decouple the channels via shear deformations than to destroy 
the channel structure altogether. This leads to expect another moving glass phase which keeps a periodicity along y, 
which has been termed Moving Transverse Glass. Since it retains a periodicity along the direction perpendicular to 
motion it shares the properties of moving glasses, and in particular it exhibit a non zero transverse critical force at 
T = 0. 

We have given predictions for the phase diagram of moving systems. It shows that the existence of the Bragg glass 
phase in the statics has profound implications on the dynamical phase diagram as well. Indeed it is natural to connect 
continuously the static Bragg glass (at v = 0) to the moving Bragg glass (at v > 0). Thus there should be a wide 
range of velocities (down from the creep region to the fast moving region) where effects associated with transverse 
periodicity (such as the transverse critical force) should be observed. We have analyzed the crossovers between the 
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Bragg glass properties and the Moving Glass properties in the region where the velocity is not large. 

Further experimental consequences should be investigated in details for vortex systems in motion. A direct mea- 
surement of the transverse I — V characteristics at low temperaturje-, would be of great interest. But consequences 
for the phase diagrams should be explored too. It was predicted ir£3 that the static rBragg glass should undergo a 
transition into an amorphous glassy state upon increase of disorder. As discussed inell the field induced transition 
observed in many experiments is the likely candidate for such a transition. A similar prediction can be made in the 
dynamics, namely that the moving Bragg glass will experience a field induced transition into an amorphous moving 
phase. A detailed investigation of these transitions may help understand the nature of the high field pinned phase. 
Indeed the nature of the transition away from the Bragg glass may change once the systemj-moves if the moving 
amorphous phase is different from the static amorphous phase. If the "vortex glass" phaseEjO exists at all in the 
statics in c? = 3, one may expect that it would not survive as smoothly as the Bragg glass once the system is set in 
motion. Another consequence is that there should be a first order melting transition at weak disorder or at large v 
upon raising temperature from a Moving Bragg glass to a flux liquid. Finally, it would be interesting to investigate 
whether the anomalous response to transverse forces could have an impact on the anomalous Hall angle. As we 
have discussed, other experimental systems, such as Wigner crystals, seem to be a promising arena to investigate the 
physics presented here. Finally it would be interesting to reinvestigate more complex CDW systems such as double 
Q or triple Q. 

Another direct experimental consequence, in the case of correlated disorder is that one should observe a "Moving 
Bose glass". Static columnar disorder in vortex systems is strong but at large velocity one should expect that the 
effective disorder becomes weaker. Thus the Bose glass driven at low temperature should have interesting properties 
such as discussed here. The resulting moving Bose glass should exhibit a transverse critical force and retain a transverse 
Meissner effect in the direction perpendicular to motion. 

The novel properties of periodic driven systems discussed in this paper also suggest many other directions of 
investigation. As for the statics one outstanding problem is to treat properly dislocations in the moving glass system. 
A controlled calculation may seem out of reach, but on the other hand the existence of elastic channels suggests a new 
and wecise way to look at the problem of generation of dislocations in moving structures and may provide a starting 



pointt2£l. Solving this issue starting from large velocity is already a formidable task, but could help us to understand 
what happens close to the threshold. Indeed here again only simple cases, inspired from the manifold or CDW with 
scalar displacements and no transverse periodicity, have been considered previously. As in the statics it is possible 
that the physics is modified in a quite surprising way, and certainly all the issues about critical behavior close to 
threshold, dynamics reordering, elastic to plastic motion transitions, will have to be reconsidered. These issues are of 
major theoretical concern but also of large practical importance. Finally we note that though Model I remains to be 
tackled in order to reach a complete description of the lattice elastic flow, we have shown that the extra linear terms 
do not seem to change drastically the main features of perturbation theory. The KPZ terms remain to be treated, 
but an interesting possibility would be that again because of periodicity their effect would be weaker than expected. 

Another interesting issue is to understand to which extent a moving, or more generally a non potential system can 
be glassy. This concept may seem doomed from the start since one could conclude that the constant dissipation in 
the system would tend to kill glassy properties. However there too the situation may be more subtle and leave room 
for unexpected behavior. We have proposed the moving glass as a first physical realization of a "dissipative glass" , 
i.e. a glass with a constant dissipation Jiale in the steady state. Other realizations of non potential. glBSsy systems 
have been studied since, in spin systemsllHZI or for elastic manifolds in random flows such as polymersEEru. 

It is important to characterize these glassy effects in driven systems. Too close analogy with glassiness in the statics 
could be misleading. As we have discussed it is possible that the existence of a transverse critical force leads to 
history dependence effects. These should be checked. The role of the temperature in moving systems andJta-relation 
to entropy production remains puzzling. It is natural to expect, as in other related non potential systemstHZrEjifj that 
the absence of fluctuation dissipation theorem leads to a generation of a temperature. This is of course what happens 
in the RG approach presented here. This heating effect however is very different from the "shaking temperature" 
(since it disappears at T = in the Moving Bragg glass) and rather is likely to be related to the entropy production. 
Hopefully the methods introduced here should allow to understand this relation better. We have found within the FRG 
that at finite temperature the physics is controlled by a new non trivial finite temperature moving glass fixed point. 
This result is strengthened by the fact that another non zero T > fixed point has |alBO been obtained in the problem 
of randomly driven polymers, a problem which does have a dissipative glassy phaseEil. The problem of understanding 
dissipative glassy systems is also related to the study of general non hermitian random operators. Indeed non potential 
dynamical problems (including e.g the moving glass) are described by a Fokker Planck operator whose spectrum is not 
necessarily real (by contrast with potential problems which are purely relaxational) . These Fokker Planck problems 
with complex spectrum, ( which could be called "dynamical non hermitian quantum mechanics"!!) are, related to 
problems which have received a renewed interest recently (such as vortex lines with tilted coL urnjiar rdefectsEHSEj, spin 
relaxation in random magnetic fieldst22l diffusion of particles and polymers in random flowgli£Hii2l'Ell. Exploring this 



56 



connexion, as well as the very interesting question of the classification of these glasses and the study of their physical 
properties is still largely open. 

Finally other fascinating open questions remain. The T — Moving Bragg glass fixed point, at least within Model II, 
is a time periodic state. The general question of the stability of periodic attractors towards chaotic motion is still very 
much open. It is related to problems of time coupling and decoupling in nopJinear dynamics, such as synchronisatiop af 



oscillators in josephson arrays, which has been studied extensively recentlyE£rlii3 or synchronization by disorderti^ 
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The relation between instability to chaos and possible non perturbative generation of a temperature is also intriguing. 
Indeed one important issue is whether dislocations when present will generate an additional temperature or chaos. 
Finally, the general question of dynamical elastic instabilities is also related to recently investigated questions about 
solid friction. It would be interesting to investigate in solid friction quantities analogous to the transverse response 
and the transverse critical force once the solid is in motion. 

We are pleased to acknowledge many useful and extensive discussions with F. Williams, D. Geshkeinbeim, L. 
Gurevich F. Pardo, M. Marchevsky, P. Kes, M. Ocio, A. Kapitulnik, M. Charalambous, C. Simon. S. Bhattacharya 
E. Andrei 



APPENDIX A: DERIVATION OF THE EQUATION OF MOTION 

We derive here the contiuum version of the equation of motion for case where it does not derive from an Hamiltonian. 
It is then more simple to perform the continuum limit in the dynamical Martin-Siggia-Rose (MSR) action. One 
introduces for each particle a conjugate field (iiij) and an interpolating field similar to the one for Ui{t) such that 
u{Rj + vt,t) = Uj{t). Using these two fields the MSG action reads 

5 = ^(iM,)[77^^ + ^- [ dV{r)5{r-R'> + vt + u,{t)) + F-7^v + Q{t)] (Al) 
j ■' 

= (^5{r- R^^+vt){iu{r,t)){dt+v.dr)u{r,t)- (A2) 

f {iu{r, t))dV{r) J2 ^{r - ^° + + u^{t)) + F - r]v + C(i?j (t), t) (A3) 

Since the interpolating fields are smooth at the scale of the interparticle distance they have no Fourier components 
outside the Brillouin zone and the the kinetic term can be written in an exact manner 

{iu{-q,tMdt+v.dq)u{q,t) + c{q)u{q,t)] (A4) 

q.BZ 

Using the Fourier decomposition of the density 

J2 -R--vt~- u,{t)) = det(9„(/.^(r)) ^ ^^K(r-vt^uWr.t)+vt,t)) (^5) 

i K 

where (j) is the labelling field (^^. The pinning force is thus given by 

^^pin = - / {^uc.{<^{r. t) + vt, t))d^Vir)det{d^Mr)) E e^^(^-''*-"(*('-'*)+''*'*» (A6) 

Jr ^ 

Up to now we have made only exact transformations. Note that the disorder term depends of u{4>{r, t) + vt, t) rather 
than M(r, t). This amount in using the origin of the displacements instead of the actual position of the vortex to label 
the displacement fields. In the elastic limit u{r,t) <C r since \ui — Ui+i\ <C a where a is the lattice spacin g. One is 
therefore entitled to replace u{(j){r, t) + vt, t) by u(r, t) ■u{(p{r, t) + vt, t) by u(r, t). An integration by part of ( A6) gives 
back the action coming from the pinning force ( |3q ) and two additional terms of the form 

(A7) 

Such term coming from the fact that the continuum limit was performed after the functional derivation with respect 
to u, contains only higher gradients, and are negligeable in the elastic limit. They can in principle however generate 
relevant terms. All the allowed relevant terms will be examined in detail in section VIII B. 
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APPENDIX B: FIRST ORDER PERTURBATION THEORY 
1. general analysis 

In this Appendix we study the perturbation theory in the disorder and compute the effective action r[u, u] to lowest 
order (i.e first order) in the interacting part Sint, using the standard formula: 

T[u,u\ = So[u, u] + {Smt[u + 6u,u + Su\)su,Su (Bf) 

where the averages in ( pf| ) over 5u, Su are taken with respect to the free quadratic action 5*0 given in (p5|). We will 
remain as general as possible, in order to treat several problems and cases simultaneously, and will specify only at the 
end to particular cases. We will thus choose the following disorder term (as it appears in the MSR action ([44|)). 

S^nt^-lf {iu^t){iu^,,,,)A"^{urt-Ur't'+v{t-t'),r-r') (B2) 

^ Jrr'tV 

These allows to treat several problems. It allows to treat short range correlated disorder keeping the cutoff dependence 
which allows to generate the extra linear and KPZ terms (see section). It also allows to treat correlated disorder. The 
disorder correlator will be chosen as: 

A^'l^iurt - Ur't' + v{t - t'), r - /) = X! - r')e-*^-("'-'-"-'''+"(*-*'» (B3) 

K 

where the symbol denotes a discrete sum of lattice harmonics for a periodic problem and a continuous sum 
= J d'^k /{2ttY for a non periodic manifold. For the model ( ^7| ) one has: 

^f{r ~ r') = K^Kpg{r ~ r')e'''^---"> (B4) 

this is the hare starting correlator (it will itself be corrected and will not remain under this form, see below). In the 
case of Model II (|8|), i.e the continumm limit of the above model, one can replace: 

g(^_/)e»^('-'-') (B5) 

This is because the scale at which the displacement field varies is large compared to the correlation length of the 
disorder (see discussion of Section ()). Since we know that non potential terms may be generated under RG (from 
FDT violation) we will from the start rather study the continuous model: 

Af{r-T')^A.f5{r-r') (B6) 

We will work at finite T and also specify to T = 0. With these definitions one finds: 

Jrt ^ Jrr'W 

with: 

SwN-- / ^i?;!f^,,,(-iX^)A^'3(r-r')e-'^-("'-*-"^'*'+"(*-*'»e-5^-^-.'-'''-^ (B7) 
D'tytM = J2 ^K^^r - ^')e-»^-(«--",-'t'+-(*-*'))e-i^-^-'."'-^ (B8) 

K 

We have used that A^^'^(r' ~ r) — A^^(r — r') which comes from simple relabeling in (). Using time and space 
translational invariance of the bare action one can also write: 



KtM^ ^''(urt~Urt',t-t',r-r') (B9) 

Jr'f 

I]"(Mrt - Ur't',t- t\r- r') = ~^W^.t^, ^_,,{-iK^)Af{r - /)e-»^-(-'-*— .'t'+''(*-*'))e-i^.s.-.'.*-*' -i^ (BIO) 

K 

Drtyt' M = ■D"^(wrt - ur>t' , t - t' , r - r' ) ^ ^k^^^ " r')e-^^-("--"-'*'+"(*-*'»e-*^-^'— (Bll) 



K 
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At T = it reads simply: 

T[u,u]= S[u,u]~ I iiu'^^)R^J^^,^,A''^'^{urt-Ur't' +v{t-t'),r-r') (B12) 



Thus to this order the effect of temperature amounts to replace everywhere formally: 

A^fir - r') ^ Af{r - r')e-*^-^'--''-'' (B13) 

Temperature has thus two important effects (i) it generates a time dependence (ii) it smoothes out the disorder. 
One can already see that there will be two important different cases. Either (high enough dimension) there is a 
time persistent part to the correlator limj^oo ^r,t = -Boo < +00, in which case the disorder is smoothed out. Or 
limt^oo Br,t — 00 and the disorder gets smaller at larger scales (low dimension). Another interpretation of the above 
result is that the corrected equation of motion includes (i) a new, non random, time retarded force Sc([u] (see () ) (ii) 
a corrected pinning force which has an extra time dependence. 

{R~^)"^J.,^,u'^,^, = SqM + Fa{r,t,Urt) + fa - rja/svp (B14) 
where the new (time dependent) pinning force correlator is: 



F„(r, t, u)Ff3{r', t' , u') = ^ Af{r - r')e-'-'^-("'-'-«.'t'+^(*-*'))e-^^-B._.',._.'.-ft: (315) 

K 

We now separate the relevant contributions in this new, complicated non linear, equation of motion. We also define: 

E'^f'iUrt - Urt',t^t',r- t') = = R^J^,^,{{A'''--fPiUrt - Urt')) " Sw [ rt" {^'"''•''^ i^rt " «W"))) (B16) 

omJ;,(, Jt" 

On the functional expression () we can identify the corrections to various terms. The first thing to do is to obtain 
the corrected response and correlation functions. For that one simply has to expand () in powers of u and u up to 
quadratic order. This yields and respectively the linear and quadratic part. The linear term proportional to u 
gives the correction to the force (5/" = — = 0]. 

a. linear part of the effective action: correction to the force (or velocity) 
The linear term in the effective action in (^4|) becomes: 

{^u'^^)^{J^-^^^pvp + 5f^{v)) (B17) 

Irt 

where the correction to the force is given by: 

5Uv)= f Rf,Y.i-^K^)Af{r)e-'^-'e-'^^■^'■^^■^ (B18) 

At T = this can also be written as: 

SU{v) = I dTdrR'^^{r,T)A°'^'^^{vT,r) (B19) 



b. quadratic part of the effective action: correction to the response and correlation fonction 
The quadratic part of the effective action reads: 

Jrr'W ^ Jrr'W 
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(i) response function 

The correction to the response fonction is thus: 



{5R-^)"'\q,uj) = - [ {l-e'^'J-^+^*^)R''^{r,t)K^KsAf{r)e-'^-''*e-i^-^^-*-^ (B20) 

Jrt 



at T = it reads: 



i5R-^)"^iq,uj)= [ (l-e'(''-''+'^*))i?T*(r,t)A"*'T'5Kr) (B21) 

Jrt 

It is useful to perform the small q, co expansion to obtain the corrections to the friction coefficient, the linear terms 
and the elastic matrix. Prom the general equation of motion () one has: 



One finds: 



i6R-'f\q,io) = {iLo)6r,c,0 + {iqp)SL'^0 + QpgaSC^j; + h.o.t (B22) 
Sr]a0 = [ ii?7/i<:^i<:;3A^*(r)e-'^-''*e-5-^-^-''-^ (B23) 

Jrt 

51"^^ = [ rPR'jK^K0A'^\r)e-'''-^*e-i''-''-'^-'' (B24) 

Jrt 

= \j rPr''RlfK^K0A.f{r)e-''^-''*e-^^-^^''-^ (B25) 

(B26) 

This is a very general expression which we will particularize to special cases below. At T = one has: 

5'q'^l^{v) = - j drdrrR'^^r, t)A"^'^^{vt, r) (B27) 

Note that ddfa{v)/dvfs = —Srjapiv). In the limit v ^ Q one finds Sfa{v) ~ —SrjapVf} using an integration by part, 
provided the function A is analytic. Non analytic A yields a critical force. 

(ii) correlation functions 

The complete uu term in the quadratic part of the effective action is D"^(0, t — t',r — r'). It allows to compute the 
corrected correlation functions using the corrected response function: 

< «-g,-a,<,a, >= [R^^\uj,q).{2nT + D{Q,u:,q))M^\-u:,-q))U (B28) 

'''''^^^'^= R-^i.,q)l8R-^i.,q) ^''^ 

where . denotes the matrix multiplication of indices. Examining the large time and space behaviours one finds that 
there are two important corrections: 

(i) a correction to the temperature (from the equal time piece): 

S{vT)a0 = 11 (i?"^(0, t, r) - D"^(0, t = +0O, r)) (B30) 

^ Jrt 

(ii) a correction to a static random force. It is identified as the time persistent part of the disorder: 

D"^(0)=lim / £>"'^(0,r) (B31) 

Jrt 

It yields to a static part 6{u) in the displacement correlation. Note however that there is additional important 
corrections to the non linear part of the disorder, which we now identify. 
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c. non linear terms in the effective action: correction to the disorder and generation of KPZ terms 

It turns out that it is important to follow not just the random force but the complete non linear static part of the 
disorder term. It is identified as: 

lim D''l^{urt-Ur't',t-t',r-r') = lim A^^(r - r')e-'^-("'-*-"'-'*'+''(*-*'»e-i^-S'-'-'.*-*'-f^ (B32) 

t—t' — *oo t—t' — >oo — ^ 

K 

Fianlly, the non linear KPZ terms can be easily seen to be generated already to this order. Expanding S^M 
second order in the field u one finds: 

= \j r^r-i?^fii-^JC^(iJC,)A«^(r)e-^^-*e-^^-^'-.*-^ (B33) 

(B34) 



2. explicit evaluation of the corrections in specific models 

a. evaluation in Model II 

We first study the continuous Model II valid in the elastic limit. It is obtained by the substitution: 

Af{r-r')^Af6{r-r') (B35) 

One finds first that: 

5^0 = Q . ^^^^^ = , dK^J^^ = (B36) 

i.c that no linear, KPZ terms arc generated, and that there is no correction to the static part of the response fonction. 
These are consequences of the statistical tilt symmetry (see Section()). The only corrections are: 

6Uv) = I Rto,t E(-^^7) A^'^e-'^-^e-^^-^o.-K (B37) 
Jt ^ 

5rj^p= I ti?;!io,t^7^/3A^'e-'^-*e-^^-^°-'-^ (B38) 

Jrt 

(B39) 

For the bare problem one can further substitute A^^ = KaKjsgK- At T = it further simplifies as: 

Sfaiv) = J2i-^K,)Af [ (B40) 

K ''1 

5rj^p = J2 K.KpAl' Ai?75^|^^^^ (B41) 

(B42) 



Note the symmetries riafj{—v) = rjafj^v) and Sfa{~v) ~ —Sfa{v). 

Let us further specify to the problem to a periodic lattice. The bare perturbation theory (i.e starting from rja/S = Vo) 
gives: 

One can also look at the dressed perturbation theory (i.e adding from the start the terms which will be generated). 
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Suppose the velocity along a principal lattice direction x. Then the symmetry y — > —y ensures in (B37) that to all 



orders r]xy = 0. However rj^x and rjyy will in general be different. Thus the tensor to be inverted will be (for d = 2): 

iVvy^afS + iVxx - r]yy)e'^e'^){iuj) + ivq^Sap + CLqaqp + criSapq'^ - q^qp) (B45) 

which gives: 

P'^i.q)aP , P^{q)a0 , {Vyy - Vxx)iuje'^e^ 



criq) + iilyyUj + ivqx CLiq) + ir]yyUJ + ivqx criq) + irjxxUJ + ivqx 



(B46) 



The correlation functions can also be computed using that here D°'^{0, t — t',r — r') = D°'^{t — t')5{r — r'). Staring 
with the bare action and to lowest order in disorder () yields at T = 0: 



it is a sum of oscillating functions, plus a static one. 
At r > the expression is more complicated. It reads: 



K I=L.T.I'=L,T''1' 



Plyjq) PisiQ) 



{2TT)Siu;-K.v)—- ^ ^ 7^75M (B48) 

BZ 'ci{q)+tT]o{uj + v.q)c'j{q)-n]o{oj + v.q) ^ ' ' 

, . . , N , / N . / — r2irjT{ijj + v.q)5riys (B49) 

ci[q) + i77o(w + v.q) c'j(q) - ti]o[Lj + v.q) 



3. Corrections in Model I : linear and KPZ terms 

We now show explicitly on the first order perturbation theory that indeed linear terms of the form C"^(m"j)9cru(!( 
with: 

Cf ^- j dTdrr^W^{r, r)A"'^'T^(wT, r) (B51) 
We now compute the various terms: 

= - Y.^iK^){iK0){iK,){iKs) I dTdrr,W'{r,T)g{r)e'''-'^'—^^ (B52) 

5fo,{v)^"^{iKo,){iKs){iK^)e'^-^^'-''^^ j dTdrg{r)W\r,T)e'^-^''-'"^^ (B53) 

Sr]'^'^{v) = - "^{iK^){iKp){iK^){iKs) j dTdrTg{r)R^\r,T)e'^-^''-''^^ (B54) 

If one assumes isotropy to start with one has: 

C^^ = ^{iK^){iKp)K'^ I dTdrr„R{r,T)g{r)e'^-^''-''^^ (B55) 

Sf„{v) ^-^{iK^)K'^ I dTdrg{r)R{r,T)e'^-^^'-''^^ (B56) 
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Sr]°''^{v) =^K^{iK„){tKi3) [ dTdrTg{r)R{r,T)e'^-^''-''^^ (B57) 

K 

Let us now study the case where v is along a principal lattice direction, assuming it remains so by renormalization 
(to be checked). One finds by symmetry that only C^^ = vi, C^^ = V2, C^^ = v^, Cy^ = V4 are non zero. Note that 
by symmetry y —y Ky ~Ky one has Sr/^y — 0, Sr]yx — 0, bjy — 

One finds: 

Cl^ = Vl=^v-Y^ I xK^K^R{r, T)g{r)e'^-^''-''^^ (B58) 
= «2 = - ^ / yK^KyK,R{r, T)g{r)e'''-^-~-^^ (B59) 
Cr ^V3^VV~Y. f ^K'KlR{r, T)ff(r)e*^ ('-''^) (B60) 

J rr 



C7 = «4 



f yK^KyK,R{r,T)g{r)e'''-^--^^^ (B61) 

Note that V2 = exactly and that to lowest order in v one has v2 = w3 and rj + Srjxxiv — 0) = rj + Sriyy{v — 0). 

Srjx^iv) = - ^ K^Kl j dTdrTg{r)R{r, T)e''<-i^-^^) (B62) 

Srjyyiv) = -J2k^K^ I dTdrTg{r)R{r,T)e'^-^''-''-'^ (B63) 



Sf,{v) = - ^ K^iK, I dTdrg{r)R{r, r)e'^ ('-"") (B64) 

Note that v'q[v) has a maximum - this may be related to the two branches of solid friction (dynamical friction 
smaller than static one). 

APPENDIX C: DYNAMICAL EFFECTIVE ACTION TO SECOND ORDER AND ANALYSIS OF 

DIVERGENCES 

In this Appendix we obtain the perturbative expression of the effective dynamical action to second order in disorder. 
At each step we will remain as general as possible so that our expressions can be applied to study a large class of 
models and situations. Then we will then study particular situations and identify the terms which correct the bare 
disorder by performing a short distance or time expansion. We will focus maialy on divergences occuring near d = 4 
(for V = 0) and d = 'i (for w > 0). A similar study in d = 2 was performed irO. Note that the operators are local in 
r but non local in time, which makes the expansion more involved. 

Note that we will study here a priori both the periodic manifold case or the non periodic one. The only difference 
is that in the periodic case one has discrete to be replaced by J^^, in the continuous case. 

The effective action to second order in the interaction term isE£2l: 

- 2Ti-\W] = {SUW + 5W]%^ - {S...,[W + SW]rs^ - ( '^"^'^^ (CI) 

with W — (u, u) and SW = (Jw, Sii) and a gaussian average over SW is performed using the bare quadratic action 
So + S2 in (??). The last term merely ensures that all one particle reducible diagrams be absent. 
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One thus has to study all possible Wick contractions of the two vertex operators: 

imposing that at least two contractions joining the two vertices 1 and 2. 

We will temporarily use the shorthand notation Urw = Urt — Urf + v{t — t'). A tedious calculation gives, for the 
iuiu term in the effective action (in short notations, omitting temporarily the r variable and all integrals): 



-2r = 2(i<^)(<=)(A„,^,([/,,,,,JA„,^,;^5(f/t.,ti))ci?:,';|TO^ - 
1 

+ 2 



+(A„,^,;,(C/,,,,,JA„,^,;,([/,,,,, - R\?^^) (C3) 

where the symbol < F\u\ > means < F[u + Su] >su and < .. >c denotes a connected average between the vertices, 
i.e < F[wi]G[u2] >c=< -F[?ii]G[u2] > — < F[ui] >< G[u2] >■ Note that simplifications will occur in the particular 
case T = since the connected terms then vanish identically, and one can also drop the averages < F[u] >= F{u). 

Using the assumption of time and space translational invariance it can be put under the form: 



-4 

as a sum of four terms: ^A = ^^eff- 



rr'tit2 



5A« = 2R'\T2,O)R'^''in,r'){A0^,^s{Ur+r',t2,t2-r2) 



6A^^} = ^d{r'')R''^{T,-r')R'\T\-r'){A„0.^s{UrMM) 

[Apx{Ur+r' ,ti-T,ti-T') + Apx{Ur+r' ,t2-T,t2-T') — 
Ap\{Ur+r',ti-T,t2-T') ~ ^px{Ur+r' ,t2-T,ti-T')]) 
^^(3) ^ iJ7.(^2,r')i?'^(ri,-r')[(A,,;5([/,,t„t,_,JA;3A;7(C^r+r',t„ti-rJ) - 

(Aqp;5 ([/r,ti,ti-Ti-T2)A/3A;7(C^r+r',t2,ti-Ti)) — 
Aap;5 (C/r,t 1 ,t2 -T2 ) A/3A;7 (C^r+r',t2 ,t2 -Ti -T2 ))] 

SAI^} = 5{r")R''PiT2,O)R'\n,-r'){A^0.,s{Ur.t„t2) 

[A\p-j{Ur+r',ti-Ti,ti-Ti-T2) — AAp;7(t/r+r',t2-ri,t2-Ti-T2)])c (C5) 



Note that some terms have vanished from the Ito time discretization property Rl, Ji RJ3 = 

''2 1 12 

We have written these terms in that form for simplicity, but one must keep in mind that in addition they must be 
symmetrized under a P and r — r when necessary. 

Up to now this is very general. We will now consider several cases. 



1. static degrees of freedom at zero temperature 

In this subsection we first set T = 0, and thus 5A(i) = SA'-^'^ = and other averages can be dropped. Here we 
will also only study static disorder and we will thus drop the v{t — t') terms thus setting Urw = Urt — Urf in all 
above formulae. The bare response function R"^{t, r) remains however arbitrary. In the moving lattice problem this 
amounts to restrict ourselves to the modes Kx = 0, which is of interest for studying the transverse components u.v = 
which see only a static disorder (assuming they can be decoupled) and v still appears in the response function. Since 
we will keep u as a vector with arbitrary number of components, the equations that we will obtain can be applied to 
other problems with static disorder (e.g the usual manifold case v = 0, periodic or not, non potential problems etc.). 

The only remaining terms are: 
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^ Jr',T,T' 

[Apx{Ur+r' ,ti-T — Ur+r',ti-T') + Ap\{Ur+r' ,t2-T ~ l*r+r',t2-r' ) — 
Ap\{Ur+r',ti-T — Ur+r',t2-T') — Apx{Ur+r' ,t2-T — W^+r' ,ii -r' )] 
(JA^?-* = R'"'{T2,r')R'^^{Ti,-r')[Aap:siUr.ti - Ur,t2-r2)A/3A;7(Ur+r',t2 " Wr+r',ii -n ) - 
AQp;5(Mr,ti — Ur,ti-ri-T2)A/3A;7(Mr+r',t2 ~ '"r+r',ti-Ti ) — 
AQp;5(Mr,ti — Ur,t2-r2) A/3A;7(Ur+r',t2 ~ lir+r',i2 -ti -T2 )] 

(C6) 

It is then easy to perform a short time, short distance operator expansion in the variables r' ,ti,T2. This yields, up 
to higher order irrelevant gradient terms, the following total correction to the random force correlator: 



SAapiu) = Aai3-'ys{u){Aa'f3'{0) ~ Aa'fj'iu)) / Gja' ir)Gsf3' {r) 

J r 

-Aaa';Siu)Af3i3,.j{u) / Gja'{r)Gsp'{-r) (C7) 



where we have defined the static response G{r) = dTR{T,r). Note that this formula is valid for a large class 
of models. It does not suppose for instance that the random force correlator is the second derivative of a random 
potential. 

It is important to note that the condition that the random force is the gradient of a potential, i.e, Aap{u) — 
—dadf}R{u) where R{u) is the correlator of the random potential (not to be confused with the response function !). 
is true only when G(r) = G{—r). Indeed, in that case, assuming the symmetry that Gaf}{r) — Gpair) one finds: 

SRiu) = ^R.^s{u)R-c.'i3'iu) - R,^siu)R-a'p'{0) J^G^c.'ir)Gsf3'ir) (C8) 

If G{r) G{—r) a non potential part is generated to the disorder. If u has only n — 1 component it will remain 
a derivative. If we study models with n > 1 component fields and a non FDT response function we will generate 
nonpotential disorder. 

Higher derivatives terms have been neglected. In the periodic case, one of them is the so-called annihilation term 
(seeBa) which produces a random force term: 

r = -J/ (*<ijv2(»uf^,j (C9) 

^ Jrtit2 

This term is important in = 2 (the so called Cardy Ostlund term) but irrelevant for d > 2. 



2. full dynamical problem at zero temperature 

In this subsection we still set T — leading to the same simplifications as in the previous subsection, but we 
keep the v{t — t' terms. So we are studying the full dynamical problem of a driven lattice (i.e with transverse and 
longitudinal displacement fields). 

The effective action is the sum of the following two terms: 

Ti - /" {iu'^t){iu%)R^P{T, r')R'^{T' ,r')A^p,^s{urt - u^' + v{t ~ t')) 

^ J rr'tt' tt' 

[Apx{Ur^r\t^T - Ur-r',t-T' + v{t' - t)) + ApA (Ur-r',t'-r - Ur-r',i'-r' + v{t' - t)) - 
ApA(Ur-r',t-T - Ur^r'.t'^T' + v{t - t' + t' - t)) - Apx(Ur-rJ'-T - Ur-r'.t-r' + v{t' - t + t' - t))] 

T2 = -i(z<J(m^+,,,,,)fi''''(^, /)i?'^r', -r')[A^p,s{ur,t - Ur.v-r + vit - + r)) 

(A/3A;7('«r+r',t' " ""r+r'.t-r' + w(t' - ^ + t')) - A^A;7(Ur+r',t' - Ur+r' ,t' -t-t' + v{t + t'))) 

-Aap-5{Ur,t - Ur,t-T-T' + v{t + t' ) ) A^A;7 (Wr+r' ,t' " Ur+r' ,t~T' + v[t' - t + t' ))] (CIO) 

We can now perform a short distance and time expansion and compute the correction to the random force correlator. 
Expressed as Aa^([/), with U = u — u' + v{t — t') it reads: 
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J q,TT' 

[Ap,(«(r' - r)) ~ i(Ap,(C/ + v{t' - r)) + Ap,(-C/ + i;(r' - r)))] + 

R^p(t, q)R'\T', q)[A^p,s{U + vT)){Apx,^{^U + vt')) - Ap^,^{v{T + r'))) 

-A^p-sHt + T'))Af3x-ji-U + vt'))] (Cll) 

It is also convenient to study the Fourier transform of the correlator: 

Ac./3(f/)=5]A^''e'^^ (C12) 

K 

and to compute the correction to A'^ . We will express it using the response function in ^^{3 = iu!,q) spatial 
Fourier transform and time Laplace transform. It is the sum of two contributions. The contribution of Pi is: 



SAf = i-iK^){-tKs)Af J2 I K^'R'^'i-ivK', ^q)R'\ivK', q) (C13) 

K' -'l 

-\T.I {-^{K-K'),){~^{K-K')s)Af_^,[AP^,R■^P{~^vK',-q)R'\^vK\q) 

^ K' •'1 

+ AP_).,R''P{ivK', -q)R^^{~ivK', q)] (C14) 
The contribution of r2 is: 

5Af HK')s)H{K' - K),)A%PA^^^,_j,R^P{tvK', q)R'\iv{K' - K), q) (C15) 

K' ■'l 

SAf = {iK)sA''/ I {-~iK%Al^,R-^P{iv{K + K'), q)R'\ivK' , q) 

K' 

-{iK)^At\ I i-^K')sA"/,R''P{zvK', q)R'\iv{K' - K), q) (C16) 

K' 



3. study at finite temperature 

In this subsection we study the case of finite T > 0. We start with the static disorder case (corresponding to 
subsection () above. The study is rather tedious and we will skip some details. Since that situation has already been 
analyzed (but applied to the different case of a periodic manifold in d = 2) we refer tcl£3 for further details. We will 
concentrate mostljz. on what is needed for the analysis near d = du {du = 4 for v = and du = 3 for v > 0). 

The result (seeo) is that the short distance expansion of the effective action up to second order in disorder produces 
a iuiu term which can be written as: 

f 5A^'^e-5^-^°''i-'2-^(i<tJ(«uf^tJe'^("'-*i-"-'^) (C17) 

which is thus of the same form as the first order term and which thus corrects it. Here again, other operators (such 
as higher gradients) are produced, but they are irrelevant near du- 

We will be using extensively the assumed symmetries A'^ = A^" = A"fj^. We are not using the potential condition 
that A'^ ^ KaKp since this is wrong in general (see discussion above). 

We find that the corrections A"^ are a priori the following, starting with the terms which will not give a contribution: 

(i) the terms with connected averages (1) and (4) give: 

6Af = -2A^''^x;i^^A^^i?g^,^i?7^';^e-i^'-^«--^'(e^-(^-i-^-i--)-^' - e^-(C-i+^2-c.,.J./f') (ci8) 

K' 

+KsAfYK^K'RZ.R'Ar,e~i'''-''°'---''\e-''-^^-'-^-^^^^ 

K' 
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One can check that this term will not produce a divergence, (check). 

(ii) the last two terms of 6A^^^ give: 

= ^(if,if;(A^^A^^)U^ - i^^i4:^(A^^A^^)U^)i?7;,^i?l-^,_,^e-^^'-^°--+--^'e-^ (^----^--)-^' (C19) 

K' 

where (.•)a/3 means symmetrization over the indices a,/3. This term was unlikely to produce a divergence for the 
same reason as above, but in any case it does not since it vanishes ! Indeed one sees on this expression that this 
quantity vanishes because of the symmetry n T2 which makes the summmand over K' odd under K' — > —K'. 

(iii) finally, the terms which will produce divergences are the term SA^"^^ and the first term of SA^^h They give a 
total contribution: 

5Af = ^ (if^K5A^'='A^^,i?I';_^i?l^^^^,e^(2Co.o-c.,-.-c,,_,).i^' (C20) 

K,K'=P-K 

+jr;jr5A^^A^^i?7_^^iii^^^^,e^-(2C^°-°-^'---^'-.--')-^') (C21) 

-PjPsAf J2 ^K'R-'r,rR%r'^~'^'''-'^'''^' COShP.(a,-r - Cr,-r')-K' (C22) 
K' 

we are using extensively the assumed symmetries A*^ = A^" = A"fj^. We are not using the potential condition 

that A'^ ~ KaKjj since this is wrong in general (see discussion above). 

This term will produce a divergence at It is simply the finite temperature generalization of () above. The idea 
is that since i?(r, t) is finite (and cutoff dependent) at large r, t and since at T = the infrared divergence came 
from the large r, r values, the new IR divergence is the same as the old one, with a coefficient obtained by simply by 
taking the large r, r limit in the exponential factors. 

Near d„ the large time or space limit of -B(r, r) = 2(Co,o — C'r.r) is proportional to the temperature: 

hm B^'^ir, r) = 2Cto = B^5^0 = 2T / (^^ + (C23) 

max(r,t)^oo ' CT{q) CT(q) 



If one assumes a circular cutoff A = 2n/a (case v = 0), one finds: 
For isotropic elasticity we simply obtain 



/ q'^-^dq\ (C24) 
Jo 9 



The final divergent contribution will be: 

5A"/ = {K^KsAfA''j^MZ^rR-\r'+^KsA"/A^j,)R^P^R^^^^ (C26) 

K,K'=P~K 

-P^PsAf A''K'RZ,rR-r,r'e-^'''-''^-''' (C27) 
K' 



4. second order corrections to the mobility 

The second order renormalization of the mobility comes from the iii terms. We will only give them in the one 
component case. They are given by: 

- 25Seff = iutA-^A'iut, - Ut'^)A"{ut, - Ut0Rt'^tARt,t[ - Rt'.tOiRt^t', - Rt[tO 
-2A"{ut, - Ut,JA'{ut, - Ut0Rt,t[Rt',tA^M', - Rt[tO + 
A"{ut, - UfJA'iut, - Ut'J{Rt,t, - i?tife)(^iti - Rt[t'JiRt.t[ - Rf^tO + 
A"'{ut, - Uf;)A{ut, - ut0Rt,r^iRt,t, - Rt[tM,t', - Rt[tO) (C28) 

Note that one additional contraction vanishes since it involves Rt^t'^Rt'^t^ = 0- 

A preliminary investigation shows that no IR divergences are found either to second order (which would mean 
theta = 1 to 0(e^)), but this has to be confirmed by a more detailed analysis. 
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APPENDIX D: ANALYSIS IN D = 2 



We start from the MSR dynamical action corresponding to equation (153) is: 

S[u, u] = So [u, ii] + S2 [u, u] + Sint [u, ii] 

Sq[u,u]= / iUrt{R''^)rtr't'Ur't' - / r]TiUrtiu 

J rr'tt' J rt 



'■rt 



S2[u, u] = i / {iUq^t)[iU-q,t'){^q^ + Ao) 
^ J a.t.t' 



q,t,t' 

Sint[u,u] = - j drdtdt''^{iurt){iurt')^{urt -Urv) (Dl) 

where 

i?"' - ridt - c^dl - Cydl + vd^ {D2) 

is the two dimensional propagator defined in 

3sence of disorder the free action 6*0 y: 

R{r,t)= [ e^«''r;-ie-('=-«'+=''«'+^''«-)*/''6»(t) = / e'lr+i^^t 



In the absence of disorder the free action 5*0 yields the following correlators C in (D8) (A = in (Dl)): 

1 



CxQx + CyQy + iVQa: + irjU) 



C{r t) = I f^-(c!cQl+c^,ql)\t \ / n ^iqr-ivq^ct / V _ I ^iqr+iujt _ ^T?] 



qCxQl + Cyql \cxql + CyQl + ivQx + iriuj\ 

B{r,t)= ( ^ 2(l-e-('=-«'+'=«««)l*l/''cos(gr + i'fcV'7)) 



q,u^ ' " \Cxqi + Cyq^ + ivqx + i-quj 



(1 - cos(gr + (D3) 



Note that for ?; > C(r, t) 7^ C(-r, t). 

In the presence of disorder one studies perturbation theory expanding in the interaction term Sint using the quadratic 
part 5*0 + ^2 as the bare action. The disorder has a quadratic part S2 which is purely static and is immaterial in the 
perturbation theory. Indeed the response function of 5o + ^2 is identical to the one of Sq and the correlation function 
is changed as: Cq.t Cq.t + Cstatq,t with: 

„ Ag^ + Ao 

which is purely static and does not appear in any diagram of perturbation theory. Thus to establish the dynamical 
RG equations one can consider that Sq is used as the bare action. 

As in section ^ we compute the effective action T[u, ii] in perturbation of Sint- To lowest order one gets 

r[u,u] = So[u, ii] + S2[u, ii] + {Sint[u + Su,u + Su])su,5u (D5) 

which gives 

T[u,u\ = Sq + S2 - I Rrt,rt'{iUrt){^' {Urt - Urt')) - \ I {iUrt){iUrt'){^{Urt - Urt')) (D6) 
Jrtt' ^ Jrtt' 

where as usual in {F[u]) one has split the fields and use So to make the averages {F[u + 5u])su- The average in ( p^ ) 
is easy to perform and gives 

5T ^ I {iUrt) gSin{Urt- Urt')Rrtrt'e~^^'-<^'-<^' -\-{iUrt){iUrt') gCO^ (D7) 
Jrtt' ^ 

where: 

Brtr't' = Crtrt + Cr't'r't' — Crtr't' — Cr't'rt (D8) 



68 



The expression (D7) immediately yields the corrections to first order in g for the friction coefficient, the temperature, 
and the disorder given in the text. The correction to 77 comes from a a gradient expansion in time which yields a 
correction to the term iiirtvdtUrt- The correction to 77T from a correction to the iiirtiurt and the correction to disorder 
from the long time limit of the exponen tial. 

To compute the RG equations from ( |156| ), we have to decide on a regularization scheme. Here we choose to take 
an infrared regulator by defining a large time tmax but no infrared regulator in momentum q. The ultraviolet cutoff 
is enforced via a gaussian cutoff in momentum, i.e we define: 

(fq 2T 



Bir,t,a) 



(27r)2 cq- 



■(1 



-c? n\t\ iqr+ivqa:Ht 



(D9) 



where the mobility /i = has been introduced. (D9) can be readily evaluated as: 



2T 

B{r,t,a) = — 
c 



+00 



ds 



<fq 



■sq 



-{s+ctlt)q 



T 

cos{vqx^t)) ^ - — 
lire 



1 



y-' + Cx + i^fit)^ 



u '1 



/, (27r)2 

where in the intermediate stage we have integrated over q and performed an intermediate change of variable u 
c^t/{c^t + s). Using 

r dt 

■™) = C + ln(rz)-Si[-rz]) Ei[x\ 



(DIO) 



du 
u 



e*- 

, t 



(Dll) 



Ei[-x\ -:r«n C + log(-a;) -x+ — + 0{x^) 



one obtains for (DIO) 



B{r,t,a) = —{Log[^^ ]+C + Log[^ ' ' 



27rc 



4(c/i|i| +a2) 



Ei\ 



-{y'^ + {x + v^itf) 
4(cAi|t| + a2) 



]) 



(D12) 



APPENDIX E: STABILITY OF ERG EQUATION 

One can show that the fixed point solution is locally attractive within the space of periodic functions on [0,1]. One 
can in fact diagonalize it completely. This allows to study crossovers in precise way. 

The function A(m) — A(0) — A(u) is positive on the interval [0,1] and satisfies the RG equation: 

^A(u) = (1 + A"H)A(«) (El) 

with the conditions A(0) = A(l) = 0. 

We have checked numerically that analytic initial functions (i.e with zero odd derivatives at 0) converge towards 
the non analytic fixed point A (u) — u{l — u)/2. The stability analysis is performed by writing A(it) = A [u) + 5{u). 
One then has to solve the eigenvalue problem: 

iu(l - u)5"{u) = ~\5{u) (E2) 

The eigenfunctions arc such that: 

]^u{l-u)5':,{u)^-]^n{n^l)8^{u) (E3) 

One can also define the variable u = (1 + v)/2. Then the eigenfunctions are the Jacobi polynomials 6n{u) — 
Pn^'~^{v) (see Abramowitz and Stegun p 779) and form an orthonormal complete set. They can be written as: 

Sniu) = ^.(1 - u) — [uil - .)]"- ='-^(1- - vr-^ (E4) 

Because of the u — > —u symmetry, which due to periodicity becomes u — > (1 — m) symmetry, we can restrict 
ourselves to n an even and non zero positive integer. The lowest eigenmodes are thus (u) = (u^ — l)/4 (eigenvalue 
-1), 64{v) = 3(1 - 6v^ + 5v^)/16 (eigenvalue -6), 6e{v) = 5(^2 - l)(l - Uv"^ + 21u4)/32 (eigenvalue -15), etc.. 
Note that they satisfy J„(w = — 1) = (5„(w = 1) = as requested. Note that all these eigenfunctions are non analytic 
(though by combining several one may get analytic ones). Rendering the initial function non analytic is presumably 
the role of the non linearity. This equation is interesting since it is the simplest case on which one can work out the 
full stability spectrum and it may enlight us about the generation of non analyticity in these type of RG equations. 
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APPENDIX F: PARAMETRIZATION OF MOVING STRUCTURES 



Moving structures can be generally parametrized by their internal space D, the number of components n of the 
displacements fields (characterizing its deformations - or the number of components of the order parameter) and the 
embedding space d. We denote for convenience by the same symbol the space itself and its dimension. 

In the statics one can distinguish several cases. The problem of manifolds in random potentials has been studied 
for e.g (i) fully oriented manifolds where D and n are orthogonal [d = D + n) such as directed polymers or interfaces 
(ii) isotropic manifolds n ~ d such as self avoiding chains in random potentials (iii) problem of lattices where usually 
d = D and n < d. Lattices with D < d are possible in principle, such as flat but fluctuating tethered mcbranes 
D — n < d OT isotropic tethered membranes D < d = n or any intermediate case (the so-called tubules). 

In the driven dynamics, let us call x the direction in the embedding space along which the system is driven. 

One can distinguish the following cases: 

(A) the structure is elastic (i.e not liquid) in the direction where it is driven. Then there is a displacement Ux along 
X and X also belongs to the n-space. There are two subcases: 

(Al) X also belongs to the internal space D. This is the problem of manifolds driven along one of their internal 
dimension, to which the moving glass studied here belongs. 
A general parametrization in that case would be 

D = {x,yi,z) n= [ux,Uy = [uy^,Uy^)) d= {x,yi,y2,z) (Fl) 

it allows for manifolds with D < d which do not entirely fill space (i.e with "height" degrees of freedom Uy2)- Then 
a parametrization of the dimensions (and the subspaces) is: 

D = l + di+dz n^l + di+d2 d = 1 + rfi + ^2 + 4 (F2) 

where di and d2 are the number of components of Uy-^ and Uy^ respectively, and dz the number of components of 
z. In this paper we have mainly studied the case d — D {d2 = 0) but with di > 0. Note that a single Q CDW has 
di = d2 = [uy = 0) and d = D. 

(A2) X does not belong to the internal space D. This is the problem of manifolds driven perpendicular to their 
internal dimension. The general parametrization in that case is: 

D = {yi,,z) n ^ {ux,Uy = {uy,,Uy^)) ( 2; = u j. , J , ^2 , ^ ) (F3) 

and thus D = di + dz, n = 1 + di + d2 and d = l + di + d2 + dz. It also contains the case of a driven order parameter 
u which does not couple at all to internal space (such as a vector spin order parameter). Indeed in that particular 
case one can define the "embedding space" as the sum d = D + n (and thus Uy^ = ) . 

(B) The structure is a liquid in the direction where it is driven. Then x belongs to Z?-space but not to n space. 
Then one sets u^; = in the case (Al) above, i.e n ~ di + c?2. The parametrization is thus: 

D={x,yi,z) n^{uy^,Uy^) d = {x,yi,y2, z) (F4) 

with D — 1 + di + dz, n ^ di + d2, d — 1 + di + d2 + dz- The Transverse Moving Glass studied here is one example 
with d2 — and (di = I, dz — 0) in d = 2 and {di = 1, = 1) in d = 3 (a moving line lattice giving a smectic 
sheets structure of channels) and {di = 2, = 0) in d = 3 (a moving point lattice giving a line crystal structure 
of channels). Note that as for any liquid scalar density fluctuations should in principle be also incorporated in a 
complete description. 



APPENDIX G: HARTREE METHOD 



For completeness we give here the Hartree equations exact in_the large N limit, for Model III generalized to N 
components. The equations at u = were derived and analyzed inEj (see alsoE£3E£il. These equations will be analyzed 
further in a future publication. 

The Hartree equations are: 



dtRkw 
dtCkw 



(fc2 + ivkx)Rktt' + 4 / dsV^'{Bts)Rts{Rku' - Rkst') 



-(P + ikxv)Cktt' + 2 / ds Vi{Bts)R-kt's 
Jo 

+4 / dsV^'{Bts)Rts{Cktt' - Cks.t') + 2r 
Jo 



(Gl) 
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where Ru' = Rkw , Btf = Jj. But' and But' = Cktt + Ckt't' - Cktt' - C-ktt': noting tlmt Cktt' = C-kt't- 

where V2 contains only the potential part of disorder while Vi contain all disorder (seeE2l for definitions). One can 
look for a time-translational invariant solution of these equations: Rktt' = rk{t — t'), Cktt' = Ck{t — t']Ei (in the statics 
this is the equivalent of a replica symmetric solution) . It can be written in Fourier transform version: 

. ^ 1 . ^ 2T + Diu;) 

Note that Cfc(tj) = c_fe(— tj). We have defined: 

Y.{uj) = -4 / dte"^%' ib{t))r{t) D{uj) = 4 / dte*"V/(6(t)) (G3) 



— 00 



where r{t) = /^^ r{t) and b{t) = /^^ bkit) - /fc(2cfc(t 0) - Ck{t) - c_fc(i)). Note that 6fc(t) = fjl - e''^*)(cfe(t^) + 
c-k{uj) and b{t) — /^(l — cos {tut)) c{uj) where c{lu) — J^. Ckitt') is an even function of uj. 

A very superficial analysis of the above equation would indicate that this non periodic problem has asymptotically 
linear rosoonse and is noit-glassy for w > for iV = 00 (while it has non linear asymptotic response for v = both at 
iV finitely and N infiniteH) Indeed the response function seems to be massive since integrating over k^ one has: 

r(uj) = - f 5-— (G4) 

24 (fc2 + ,^ + S(0)-E(c.) + ^f)i/2 

Thus the response to an applied force being F/V = (1 — S'(0)) would be linear at least at_the most naive level (for 



a more detailed behaviour one must add a small transverse force and follow the methods olE£3) . This is related to the 



absence of divergence for 77 noticed in the FRG Section VI . Further investigations would be necessary however before 



reaching a conclusion. One should makes sure that no transition occur in the above equations (such can happen in 



the case v = 0). Also, it is possible that the glassy physics found in Section VI which comes from a renormalization 



of the disorder may not be fully captured here by the most naive large N limit. 



APPENDIX H: DIRECT METHOD OF RENORMALIZATION 



We present here a simple method to obtain the FRG equations for Model III at T = 0, based on a mode elimination 
by hand. This approach is less rigorous than the full calculation of the effective action performed in (^). However it 
is quite direct and illustrate a simple way how non potential disorder is generated in this problem. 

We start from: 

iG-Xr,4' = Fo.ir,Ur) (HI) 
and separate u ^ u + 5u. The equation for the fast modes 5u is: 

{,G-Xr'54' = Fo.{r, ur + Sur) « F„(r, u) + F^.pir, u)Su^ (H2) 
to lowest order. Solving to the same order one has: 

<5< = C^^^^Fpiruu) + G?,^^Fp,^{n,u)Gj'^,^Fs{r2,u) (H3) 
The equation for the slow modes is now: 

(G'~^)"^u(', = Fa(r,Ur + Sur) « Fa{r, u) + Fa^pir, u)Su'^ + ^Fa;i3^{r,u)6u'^6uJ 

« Fa + Fi^^ + Fi'^ + Fi^) (H4) 

One finds: 

Fi'^ = F^.,^ir,u)G^,lF^ 

fi^^ = Fa,.Ar,u)G^,lF^.sin,u)Gi^^.^Fx 

Fi^^ = \Fc.-13^ ir, u)GP4^ Fs in , u)G1^^ Fx (r2 , u) (H5) 
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It is easy to see that if Gr.r' 7^ Gr'^r the above forces are non potential even if the original one is the gradient of a 
potential. One can now compute the new disorder correlator Acq/(m — u' ,r — r') which has several contributions: 



^ia' = Fo.' (r', (r, u) - -A„,a(7. - u')A^^MO)G?:^Gtr' 

^al = F^,{r',u')Fi^\r,u) - ^(6rr'A^^,,p^{u - u')Asx{Q)G^^^^G^^^^ + 

Ao,,siu - u')A^x.0^iO)G'^J,GJ^ + Aa,,xiu - u')A^s-p^iO)G^^^Gj^,) (H6) 
Putting everything together, in the long wavelength limit: 

SA^^, [u - u')6rr. - A^Z + Ai^'^, + Ait + ^^1' + Ait (H7) 
one notes that A^^^^ cancels the last two terms of A^''^ and thus finds: 

5Aai3{u-u') = Aai3-^s(u){Aa^l3^{0)- Aa2P^{u)) I 5G~^a2ir)5G sp^{r) 

J r 

-A^o..,s{u)App,,^{u) / 5G^o..{r)5Gsp'{-r) (H8) 



One thus recovers the result of (^. The time-dependent problem at T = can also be analyzed in t he s ame way. 
All the above equations up to (H5) are identical with r replaced by rt and Grr' replaced by Rrtr't'- In (H6) one gets 
integrals over times and again the zero- frequency Grr' ■ 
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